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Introduction 


Background 

Prostate  Cancer.  Every  year  over  180,000  American  men  are  diagnosed  with  prostate  cancer  [1], 
Although  most  will  receive  some  type  of  treatment,  such  as  chemotherapy  and/or  radiotherapy,  mortality 
for  the  10%  of  patients  with  recurrences  and  or  metastases  is  nearly  100%. 

MicroRNAs  and  Cancer.  MicroRNAs  are  regulatory,  non-protein-coding,  endogenous  RNAs  that  have 
recently  gained  considerable  attention  in  the  scientific  community.  They  are  18-24  nucleotides  in  length, 
and  are  thought  to  regulate  gene  expression  through  translational  repression  by  binding  to  a  target 
mRNA  [2-4],  Many  microRNAs  are  defined  to  be  either  tumor  suppressors  or  oncogenes  [5-7]  and  play 
a  crucial  role  in  variety  of  cellular  processes  such  as  cell  cycle  control,  apoptosis,  haematopoiesis  and 
tumorgenesis  [8,9]. 

Recent  studies  indicate  that  selected  miRNAs  may  play  a  role  in  human  cancer  pathogenesis  [6,10]. 
The  results  of  large-scale  miRNA  profiling  studies,  using  normal  and  cancer  tissues,  suggest  that  a 
number  of  microRNAs  are  either  overexpressed  or  downregulated  in  tumors  [11-18].  The  genes 
encoding  mir-15  and  mir-16  are  located  on  chromosome  13ql4,  a  region  that  is  deleted  in  a  majority  of 
B-cell  chronic  lymphocytic  leukemias  (B-CLL),  suggesting  that  mir-15  and  mir-16  may  function  as 
tumor  suppressors.  Let-7  miRNA  family  members  are  known  to  downregulate  the  oncogene  RAS  [19]. 
In  contrast,  several  microRNAs  are  upregulated  in  cancer  and  may  function  as  oncogenes.  Members  of 
the  miR-17  cluster  are  overexpressed  in  cancer  via  c-Myc  and  function  as  oncogenes,  leading  to 
dysregulation  of  downstream  targets  that  mediate  cell  cycle  and  apoptosis  [20].  Additionally,  miR-21, 
overexpressed  in  prostate,  breast,  lung  and  colon  cancers,  promotes  cell  transformation  by  targeting 
PDCD4  [21], 


Expression  of  specific 
microRNA  may  be  of  prognostic 
value  as  determinants  of  clinical 
aggressiveness  [12,15,22],  Thus, 
microRNA  expression  profiles  may 
serve  as  a  promising  new  class  of 
cancer  biomarkers.  Prostate  cancer 
microRNA  profiling  studies  by 
several  groups  [12,23,24]  revealed 
significant  differential  expression 
patterns  of  several  microRNAs 
between  normal  and  neoplastic 
tissues. 

EZH2  and  Cancer.  EZH2,  a 
Polycomb  Group  (PcG)  protein, 


EZH2 

-H3K27me3 
«-Ace-H3K27 
«-H3K27me1 
«-  total  H3 
-GAPDH 


TargetScan 


Fig  1.  miR-101  regulates  EZH2  expression.  (A)  Venn  diagram 
displaying  miRNAs  computationally  predicted  to  target  EZH2  by 
use  of  bioinformatic  software.  (B)  miR-101  downregulates  EZH2 
expression  and  decreases  H3K27me3  levels  in  DU145  cells,  while 
Ace-H3K27  and  H3K27mel  levels  are  not  altered  by  miR-101. 
GAPDH  and  total  H3  were  used  as  loading  controls. 


4 


plays  a  master  regulatory  role  in  controlling  important  cellular  processes  such  as  maintaining  stem  cell 
pluripotency  [25-27],  cell  proliferation  [28,29],  early  embryogenesis  [30],  and  X  chromosome 
inactivation  [31].  EZH2  functions  in  a  multi-protein  complex  called  Polycomb  Repressive  Complex  2 
(PRC2),  which  includes  SUZ12  (Suppressor  of  Zeste  12)  and  EED  (Embryonic  Ectoderm  Development) 
[32,33].  The  primary  function  of  the  PRC2  protein  complex  is  to  tri-methylate  histone  H3  lysine  27 
(H3K27)  at  target  gene  promoters  to  enact  epigenetic  silencing  [34,35].  Mounting  evidence  suggests 
that  EZH2  has  properties  that  are  consistent  with  those  of  an  oncogene,  as  overexpression  promotes  cell 
proliferation,  colony  formation,  increased  invasion  of  benign  cells  in  vitro  [28,29,36],  and  induction  of 
xenograft  tumor  growth  in  vivo  [37].  Likewise,  knock-down  of  EZH2  in  cancer  cells  results  in  growth 
arrest  [28,37]  as  well  as  diminished  tumor  growth  [34]  and  metastasis  in  vivo  [38], 


Fig  2.  The  role  of  miR-101  in  regulating  cell  proliferation,  invasion  and  tumor  growth.  (A)  miR-101 
overexpression  reduces  cell  proliferation.  Cell  growth  assay  of  DU145  cells  treated  with  either  precursor  miR- 
101  or  siRNA  targeting  EZH2.  (B)  miR-101  expression  decreases  cell  invasion  of  DU145  cells.  Cells  were 
transfected  with  miR-101,  EZH2-specific  siRNA,  control  miR,  and  non-targeting  siRNA.  All  cells  were 
subjected  to  a  matrigel  invasion  assay.  (C)  Over-expression  of  miR-101  attenuates  prostate  tumor  growth. 
Over-expression  of  miR-101  reduces  DU  145  tumor  growth  in  a  mouse  xenograft  model.  Plot  of  mean  tumor 
volume  trajectories  over  time  for  the  mice  inoculated  with  miR-101  (red)  and  vector  (green)  expressing  DU145 
cells. 

miR-101  represses  EZH2  expression.  In  order  to  examine  whether  microRNAs  regulate  EZH2 
expression  during  cancer  progression,  we  used  prediction  software  programs  to  nominate  candidate 
microRNAs  and  found  that  only  miR-101  and  miR-217  are  predicted  to  bind  to  EZH2  3’UTR  by  all 
four  of  the  programs  we  used  (PicTar  [39],  TargetScan  [40],  miRanda  [41],  and  miRInspector  [42]) 
(Fig.lA),  with  miR-101  as  the  top  ranked  candidate  by  the  four  programs,  despite  using  different 
algorithms.  Upon  ectopic  overexpression  of  miR-101,  miR-217,  and  control  microRNAs  into  the 
prostate  cancer  cell  line  DU  145,  we  found  that  EZH2  protein  expression  was  only  repressed  by  miR-101, 
equal  to  EZH2  siRNA  duplex,  but  not  by  miR-217  compared  to  control  microRNA,  or  siRNA  duplex  as 
measured  by  immunoblot  analysis  (Fig  IB).  Further,  H3K27  tri-methylation  levels  were  decreased  by 
miR-101  and  EZH2  siRNA  duplex,  while  Ace-H3K27  and  H3K27mel  levels  were  not  altered. 
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miR-101  inhibits  oncogenic  potential.  To  determine  whether  miR-101  inhibits  the  function  of  EZH2 
and  PRC2,  we  evaluated  cell  proliferation  and  invasion  following  transfection  of  miRNA  precursors.  As 
predicted,  miR-101  overexpression  significantly  attenuated  DU145  prostate  cancer  cells 
proliferation  and  invasion.  Co-expression  of  mutant  EZH2  (without  endogenous  3’UTR)  rescued 
the  miR-101  mediated  inhibition  of  cell  growth  and  invasion  (Fig  2A,  B).  Furthermore,  in  vivo, 

DU  145  cells  stably  overexpressing  miR-101  grew  markedly  slower  than  the  vector  control  xenograft 
(P=0.0001),  (Fig  2C).  All  of  these  findings  demonstrated  that  miR-101  was  a  tumor  suppressor  and 
inhibited  these  cancer  cell  specific  properties. 


miR-101  expression  is  inversely  correlated 
with  EZH2  expression  in  cancer  progression. 

We  examined  miR-101  expression  in  human 
tissues,  and  observed  that  miR-101  and  EZH2 
expression  were  inversely  correlated  during 
prostate  cancer  progression  (Fig  3A).  miR- 
101  expression  was  significantly  decreased  in 
metastatic  prostate  cancer  relative  to  clinically 
localized  disease  or  benign  adjacent  prostate 
tissue,  while  EZH2  was  significantly 
upregulated  in  metastatic  prostate  cancer 
relative  to  clinically  localized  disease  or  benign 
adjacent  prostate  tissue  (P<0.0001).  This  is 
consistent  with  the  functional  observation  based 
on  genomic  PCR,  2  of  16  clinically  localized 
prostate  cancers  and  17  of  33  metastatic  prostate 
cancers  exhibited  loss  of  the  miR- 101-1  locus 
(Fig.  3B).  To  formally  demonstrate  that 
genomic  loss  of  miR-101  loci  was  somatic  in 
nature,  we  identified  9  metastatic  prostate 
cancers  that  exhibited  loss  of  miR- 101-1  and 
obtained  DNA  from  matched  normal  tissue.  As 
expected,  8  of  9  cases  exhibited  a  marked 
decrease  in  relative  levels  of  miR- 101-1  copy 
number  in  cancer  compared  to  matched  normal 
tissue  (Fig.  3C).  In  addition,  our  meta-analysis 
of  the  majority  of  publicly  available  microRNA 
expression  datasets  suggests  that  miR-101  is 
significantly  under-expressed  in  prostate, 
breast,  ovarian,  lung  and  colon  cancers 
[14,15,22]. 


-4  _ 

miR-217 


miR-101 


-iJ _ : _ 

EZH2 


B 


miR-101-1 


Fig.  3.  Genomic  loss  of  the  miR-101  locus  may  explain 
overexpression  of  EZH2  in  solid  tumors.  (A)  miR-101 
transcript  levels  are  inversely  correlated  with  EZH2 
expression  in  prostate  cancer  progression.  (B)  Genomic  PCR 
of  miR-101-1.  Vertical  axes  represent  log  (base  2)  relative 
quantification  values;  dashed  lines  are  shown  at  the  deletion 
threshold  of  log2  (0.7)~-0.51.  For  clarity,  points  have  been 
horizontally  displaced  within  each  sample  class.  (C) 
Evidence  that  the  miR-101-1  locus  is  somatically  lost  in 
tumor  samples  relative  to  matched  normal  samples.  Nine 
metastatic  prostate  cancers  were  chosen  that  have  copy 
number  loss  in  the  miR-101-1  locus,  and  matched  normal 
tissues  were  analyzed  for  comparison.  Bar  heights  represent 
differences  in  log2  (RQ)  values  between  metastatic  and 
matched  normal  tissues. 
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Research  progress: 


Specific  Aim  1:  To  characterize  miR-101  as  a  biomarker  of  metastatic  prostate  cancer,  and  the 
role  of  miR-101  in  prostate  cancer  progression. 


1)  miR-101  regulates  PRC1  proteins  through  PRC2  and  other  miRNAs. 


Through  microRNA  profiling,  we  identified  that  the  miR-101  target,  EZH2,  could  repress  over  100 
miRNAs.  Among  these  EZH2 -regulated  miRNAs,  we  found  that  miR-181b,c,  miR-200b,c  and  miR-203 
could  repress  PRC1  proteins  BMI1  and  RING2,  as  well  as  the  levels  of  their  substrate  ubiquityl-H2A- 
K1 19  (Fig.  4). 
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Fig.  4.  PRC2-regulated  miRNAs  repress  PRC1 
proteins  BMI1  and  RING2. 

(A)  miRNA  profiling  of  DU145  prostate  cancer 
cells  in  which  EZH2  was  knocked-down 
compared  to  DU  145  cancer  cells  relative  to 
benign  cells  HME,  PrEC,  RWPE  and  H16N2. 
Shades  of  red  represent  increased  gene  expression 
while  shades  of  green  represent  decreased 
expression.  (B)  Overexpression  of  indicated  miRs 
in  DU145  and  BT-549  cells,  and  expression  of 
PRC  components  and  PRC1  target  histone  mark 
and  genes  by  immunoblot  analysis.  P-actin  was 
used  as  loading  control. 


Furthermore,  introduction  of  miR-101  or  EZH2  siRNA  restored  the  repression  of  the  miRNAs  by  EZH2 
(Fig.  5). 
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Fig.  5.  miR-101  and  EZH2  siRNA  increase  miR-181b,c  miR-200a,b,c  and  miR-203  expression  in  BT-549 
and  DU145  cells,  but  not  control  miRNAs  miR-217  or  miR-219. 


Most  importantly,  miR-101  target  EZH2  is  negatively  correlated  with  miR-101,  miR-181b,c,  miR- 
200b, c  and  miR-203  during  prostate  cancer  progression  (Fig.  6A).  A  proposed  model  of  coordinated 
regulation  of  PRC1  and  PRC2  complexes  in  cancer  is  presented  in  Fig.  6B. 
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Fig.  6.  Coordinated  expression  of  PcG  proteins  and  PRC  regulatory  miRNAs  during  prostate  cancer 
progression.  (A)  Expression  of  indicated  miRs  as  assessed  by  q-PCR  in  benign  prostate,  clinically  localized 
prostate  cancer,  and  metastatic  prostate  cancer  tissues.  (B)  A  proposed  model  role  for  microRNAs  in 
regulating  PRCs.  Specifically,  PRC2  is  molecularly  linked  to  PRC1  via  a  set  of  regulatory  miRs. 
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Fig.  7.  EZH2  and  STMN1  are  positively  correlated  in  prostate  cell  lines  (A),  and  prostate  tissue  specimens  (B). 
RNA-seq  was  performed  using  prostate  cell  lines  and  tissue  specimens. 


2)  miR-101  targets  STMN1  (Stathmin)  in  prostate  cancer 


From  our  RNA-seq  data  analysis,  we  found  that  STMN1  and  EZH2  are  positively  correlated  in  prostate 
cell  lines  (Fig.  7A)  and  cancer  patient  tissues  (Fig.  7B).  Interestingly,  we  found  that  STMN1  3’UTR  has 
a  conserved  binding  site  of  miR-101  by  TargetScan  (targetscan.org)  (Fig.  8A).  Next  we  cloned  the 
predicted  miR-101  binding  site  of  STMN1  3’UTR  as  well  as  a  mutant  versions  into  a  luciferase  reporter 
construct,  and  examined  whether  miR-101  regulates  STMN1  3’UTR.  As  predicted,  miR-101  decreased 
wild-type  STMN1  3’UTR  luciferase  activities,  but  not  the  mutant  (Fig.  8B),  confirming  that  miR-101 
directly  binds  to  the  STMN1  3’UTR  and  regulates  STMN1  expression  levels  during  prostate  cancer 
progression. 
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Fig.  8.  miR-101  binds  to  STMN1  3’UTR.  (A)  Schematic  of  miR-101  binding  site  in  the  STMN1  3’UTR  (left 
panel),  and  derived  mutant  which  miR-101  may  not  bind  to.  (B)  DU145  prostate  carcinoma  cells  were  co¬ 
transfected  with  miR-101  or  control  miR  together  with  a  reporter  vector  containing  firefly  luciferase  linked  to 
the  3’UTR  of  EZH2.  *p<0.001.  S.E.M.  of  3  replicates.  The  p-value  was  calculated  between  control  miR  and 
miR-101. 
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Specific  Aim  2:  To  identify  the  mechanisms  by  which  miR-101  is  dysregulated  in  prostate 
cancer. _ 

1)  miR-101  binding  sites  in  EZH2  3’UTR  were  not  mutated  in  prostate  cancer. 

Our  preliminary  data  show  that,  in  some  of  our  prostate  cancer  samples,  miR-101  was  underexpressed, 
but  miR-101  genomic  loci  were  not  deleted.  To  investigate  whether  miR-101  target  EZH2  3’UTR 
regions  are  mutated  in  metastatic  prostate  cancers,  we  employ  the  next  generation  sequencing  powered 
by  Solexa  technology  to  sequence  the  cohort  of  prostate  cancer  tissue  specimens.  As  shown  in  Table  1, 
we  found  that  10  cancer  specimens  have  point  mutations  at  3  positions  but  none  were  within  the  miR-101 
binding  sites. 


Table  1  EZH2  3’UTR  point  mutations 

Position 

Variant 

Cancer  Samples 

Chr  7:  148135656 

A— >G 

2 

Chr  7: 148135655 

G— >A 

6 

Chr  7:  148135650 

G— >A 

2 
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Fig.  9  EZH2  represses  miR-101  expression.  (A)  Depletion  of  EZH2  in  DU145  cells  restores  miR-101  levels. 

(B)  EZH2  inhibitor  DZNep  increased  miR-101  levels  in  DU145  cells  at  two  different  doses  and  time  points. 

(C)  ChIP-qPCR  analysis  of  H3K27me3  at  indicated  genes  and  microRNAs  in  BT-549  cells  treated  with 
SAHA  and/or  5-aza-dC.  All  bar  graphs  are  shown  with  +SEM. 


2)  EZH2  represses  miR-101  through  its  promoter  regions. 

Through  miRNA  profiling,  we  found  that  miR-101  is  a  target  of  EZH2  (Fig.  4A).  To  validate  this 
finding,  we  performed  miRNA  qPCR,  and  found  that  miR-101  levels  were  increased  in  EZH2  stable 
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knockdown  DU145  cells  compared  to  control  cells  (Fig.  9A).  Furthermore,  when  we  treated  DU145 
cells  with  EZH2  inhibitor  DZNep,  miR-101  levels  were  increased  in  a  concentration  and  incubation  time 
dependent  manner  (Fig.  9B).  To  confirm  that  EZH2  regulates  miR-101  by  epigenetic  repression,  we 
performed  chromatin  immunoprecipitation  (ChIP)  assays  with  anti-H3K27me3  antibodies  in  BT-549 
cells.  As  expected,  H3K27me3  occupies  miR-101  promoter  regions  as  well  as  known  PRC2  target 
promoter  regions  (Fig.  9C).  Treatment  of  5-aza-dC  and  SAHA,  either  alone  or  in  combination,  markedly 
decreased  the  H3K27me3  occupancy  in  these  regions  (Fig.  9C). 


Key  Research  and  Training  Accomplishments: 

❖  We  identified  miRNAs  that  are  regulated  by  miR-101  and  EZH2. 

♦♦♦  miR-101  and  EZH2-regulated  miRNAs  in  turn  regulate  PRC1  proteins,  BMI1  and  RING2. 
♦♦♦  miRNAs  and  PRC  protein  levels  are  inversely  correlated  in  prostate  cancer. 

❖  PRC1  and  PRC2  activities  are  coordinately  regulated  via  miRNAs. 

♦♦♦  EZH2  can  in  turn  repress  miR-101. 

❖  miR-101  binding  sites  in  EZH2  3’UTR  are  mutated  in  prostate  cancer 

❖  STMN1  is  a  miR-101  target  in  prostate  cancer. 


REPORTABLE  OUTCOMES: 
Manuscript  Published: 


1)  Cao  O,  Dahiya  A,  Chinnaiyan  AM.  Histone  Methyltransferase  EZH2,  A  Master  Regulator  in 
Cancer.  Review.  Cell  Research.  2012  In  press 

2)  Cao  Q,  Mani  RS,  Ateeq  B,  Dhanasekaran  SM,  Asangani  IA,  Prensner  JR,  Kim  JH,  Brenner  JC, 
Jing  X,  Cao  X,  Wang  R,  Li  Y,  Dahiya  A,  Wang  L,  Pandhi  M,  Lonigro  RJ,  Wu  Y-M,  Tomlins  SA, 
Palanisamy  N,  Qin  Z,  Yu  J,  Maher  CA,  Varambally  S,  Chinnaiyan  AM.  Coordinated  regulation 
of  Polycomb  group  complexes  through  microRNAs  in  cancer.  Cancer  Cell.  2011  Aug 
16;20(2):187-99. 

3)  Prensner  JR,  Iyer  MK,  Balbin  OA,  Dhanasekaran  SM,  Cao  O,  Brenner  JC,  Laxman  B,  Asangani 
IA,  Grasso  C,  Kominsky  HD,  Cao  X,  Jing  X,  Wang  X,  Siddiqui  J,  Wei  JT,  Robinson  D,  Iyer  HK, 
Palanisamy  N,  Maher  CA,  Chinnaiyan  AM.  Transcriptome  Sequencing  Identifies  PCAT-1,  a 
Novel  lincRNA  Implicated  in  Prostate  Cancer  Progression.  Nat.  Biotechnol.  201 1  Jul 
31;29(8):742-9.  doi:  10.1038/nbt.l914. 
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4)  Ateeq  B,  Tomlins,  SA,  Laxman  B,  Asangani  IA,  Cao  O,  Cao  X,  Yong  L,  Wang  X,  Feng  FY, 
Pienta  KJ,  Varambally  S,  Chinnaiyan  AM.  Therapeutic  targeting  of  SPINK  1 -positive  prostate 
cancer.  Sci.  Transl.  Med.  2011  3:  72ral7  (2011).  PMID:  21368222. 


Seminars  and  Extramural  Invited  Presentations 

1)  American  Association  for  Cancer  Research  (AACR)  102nd  Annual  Meeting,  “An  onco-protein 
axis  linking  polycomb  repressive  complex  2  and  polycomb  repressive  complex  1  through 
miRNAs  in  cancer”,  Orlando,  FL,  April  2-6,  2011. 

2)  Department  of  Molecular  and  Cellular  Oncology,  The  University  of  Texas  MD  Anderson 
Cancer  Center,  “Coordinated  Regulation  of  Polycomb  Repressive  Complexes  Through  miRNAs 
in  Cancer”,  Houston,  TX,  September  1,  2011. 

3)  Division  of  Life  Sciences,  Hong  Kong  University  of  Science  and  Technology,  “Coordinated 
Regulation  of  Polycomb  Repressive  Complexes  Through  miRNAs  in  Cancer”,  Hong  Kong, 
China,  September  27,  2011. 

4)  Center  for  Inflammation  and  Epigenetics,  The  Methodist  Hospital  Research  Institute, 
“Coordinated  Regulation  of  Polycomb  Repressive  Complexes  Through  miRNAs  in  Cancer”, 
Houston,  TX,  December  14,  2011. 

5)  Institute  for  Integrative  Genome  Biology,  University  of  California  Riverside,  “Coordinated 
Regulation  of  Polycomb  Repressive  Complexes  Through  miRNAs  in  Cancer”,  Riverside,  CA, 
January  30,  2012. 

6)  Department  of  Pathology,  Mount  Sinai  School  of  Medicine,  “Coordinated  Regulation  of 
Polycomb  Repressive  Complexes  Through  miRNAs  in  Cancer”,  New  York,  NY,  February  29, 
2012. 

7)  Fels  Institute  for  Cancer  Research  &  Molecular  Biology,  Temple  University,  “Coordinated 
Regulation  of  Polycomb  Repressive  Complexes  Through  miRNAs  in  Cancer”,  Philadelphia,  PA, 
March  21,  2012. 

8)  School  of  Biomedical  Sciences,  Chinese  University  of  Hong  Kong,  “Coordinated  Regulation  of 
Polycomb  Repressive  Complexes  Through  miRNAs  in  Cancer”,  Hong  Kong,  China,  March  13, 
2012. 

Abstracts  and  Posters 

1)  Cao  Q,  Mani  RS,  Ateeq  B,  Dhanasekaran  SM,  Asangani  IA,  Prensner  JR,  Kim  JH,  Brenner  JC, 
Jing  X,  Cao  X,  Wang  R,  Li  Y,  Dahiya  A,  Wang  L,  Pandhi  M,  Lonigro  RJ,  Wu  Y-M,  Tomlins  SA, 
Palanisamy  N,  Qin  Z,  Yu  J,  Maher  CA,  Varambally  S,  Chinnaiyan  AM.  An  onco-protein  axis 
linking  polycomb  repressive  complex  2  and  polycomb  repressive  complex  1  through  miRNAs  in 
cancer.  2011  Innovative  Minds  in  Prostate  Cancer  Today  (IMPaCT)  Conference  (DOD  Prostate 
Cancer  Research  Program),  Orlando,  FL.  Mar.  9-12,  2011 

2)  Cao  Q,  Mani  RS,  Ateeq  B,  Dhanasekaran  SM,  Asangani  IA,  Prensner  JR,  Kim  JH,  Brenner  JC, 
Jing  X,  Cao  X,  Wang  R,  Li  Y,  Dahiya  A,  Wang  L,  Pandhi  M,  Lonigro  RJ,  Wu  Y-M,  Tomlins  SA, 
Palanisamy  N,  Qin  Z,  Yu  J,  Maher  CA,  Varambally  S,  Chinnaiyan  AM.  Coordinated  regulation 
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of  Polycomb  group  complexes  through  microRNAs  in  cancer.  Multi-Institutional  SPORE 
Prostate  Cancer  Program  Retreat,  Fort  Lauderdale,  FL.  Mar.  21-23,  2011 

3)  Cao  Q,  Mani  RS,  Ateeq  B,  Dhanasekaran  SM,  Asangani  IA,  Prensner  JR,  Kim  JH,  Brenner  JC, 
Jing  X,  Cao  X,  Wang  R,  Li  Y,  Dahiya  A,  Wang  L,  Pandhi  M,  Lonigro  RJ,  Wu  Y-M,  Tomlins  SA, 
Palanisamy  N,  Qin  Z,  Yu  J,  Maher  CA,  Varambally  S,  Chinnaiyan  AM.  An  onco-protein  axis 
linking  polycomb  repressive  complex  2  and  polycomb  repressive  complex  1  through  miRNAs  in 
cancer.  American  Association  for  Cancer  Research  (AACR)  102nd  Annual  Meeting,  Orlando, 

FL.  Apr.  2-6,  2011 

4)  Mani  RS,  Tomlins  SA,  Brenner  JC,  Wang  L,  Iyer  M,  Yu  J,  Cao  Q,  Callahan  K,  Ghosh  A,  Nyati 
MK,  Varambally  S,  Palanisamy  N,  Chinnaiyan  AM.  The  origins  and  functional  consequences  of 
ETS  gene  fusions  in  prostate  cancer.  American  Association  for  Cancer  Research  (AACR)  102nd 
Annual  Meeting,  Orlando,  FL.  Apr.  2-6,  2011 

5)  Prensner  JR,  Iyer  MK,  Balbin  OA,  Dhanasekaran  SM,  Cao  Q,  Brenner  JC,  Asangani  IA,  Grasso 
C,  Kominsky  HD,  Cao  X,  Jing  X,  Siddiqui  J,  Wei  JT,  Palanisamy  N,  Robinson  D,  Iyer  H,  Maher 
CA,  Chinnaiyan  AM.  Discovery  and  characterization  of  PCAT-1,  a  novel  lincRNA  implicated  in 
prostate  cancer  tumorigenesis.  American  Association  for  Cancer  Research  (AACR)  102nd 
Annual  Meeting,  Orlando,  FL.  Apr.  2-6,  2011 

6)  Ateeq  V,  Tomlins  SA,  Laxman  B,  Asangani  IA,  Cao  Q,  Cao  X,  Li  Y,  Feng  FY,  Pienta  KJ, 
Varambally  S,  Chinnaiyan  AM.  Therapeutic  targeting  of  S PINK  1 -positive  prostate  cancer. 
American  Association  for  Cancer  Research  (AACR)  102nd  Annual  Meeting,  Orlando,  FL.  Apr. 
2-6,  2011 

7)  Iyer  MK,  Prensner  JR,  Balbin  OA,  Dhanasekaran  SM,  Cao  Q,  Brenner  JC,  Asangani  IA,  Grasso 
C,  Kominsky  HD,  Cao  X,  Jing  X,  Siddiqui  J,  Wei  JT,  Palanisamy  N,  Robinson  D,  Iyer  H,  Maher 
CA,  Chinnaiyan  AM.  Transcriptome  sequencing  identifies  novel  non-coding  RNAs  associated 
with  prostate  cancer  progression.  American  Association  for  Cancer  Research  (AACR)  102nd 
Annual  Meeting,  Orlando,  FL.  Apr.  2-6,  2011 

8)  Wang  XS,  Shankar  S,  Dhanasekaran  SM,  Ateeq  B,  Sasaki  AT,  Jing  X,  Robinson  D,  Cao  Q, 
Gopalan  A,  Reuter  VE,  Cantley  LC,  Rubin  MA,  Palanisamy  N,  Varambally  S,  Chinnaiyan  AM. 
Characterization  of  KRAS  rearrangements  in  metastatic  prostate  cancer.  American  Association 
for  Cancer  Research  (AACR)  102nd  Annual  Meeting,  Orlando,  FL.  Apr.  2-6,  2011 

9)  Brenner  JB,  Ateeq  B,  Li  Y,  Yocum  A,  Cao  Q,  Asangani  IA,  Patel  S,  Liang  H,  Yu  J,  Palanisamy 
N,  Siddiqui  J,  Yan  W,  Wang  X,  Cao  X,  Mehra  R,  Basrur  V,  Lonigro  R,  Yang  J,  Tomlins  SA, 
Maher  CA,  Elenitoba- Johnson  K,  Hussain  M,  Navone  NM,  Pienta  KJ,  Varambally  S,  Feng  FY, 
Chinnaiyan  AM.  Mechanistic  rationale  for  inhibition  of  Poly(ADP-Ribose)  Polymerase  in  ETS 
gene  fusion  positive  prostate  cancer.  American  Association  for  Cancer  Research  (AACR)  102nd 
Annual  Meeting,  Orlando,  FL.  Apr.  2-6,  2011 

10)  Mani  RS,  Tomlins  SA,  Brenner  JC,  Wang  L,  Iyer  M,  Yu  J,  Cao  Q,  Callahan  K,  Ghosh  A,  Nyati 
MK,  Varambally  S,  Palanisamy  N,  Chinnaiyan  AM.  The  origins  and  functional  consequences  of 
ETS  gene  fusions  in  prostate  cancer.  Eighteenth  Annual  Prostate  Cancer  Foundation  Scientific 
Retreat,  Lake  Tahoe,  NV.  Sep  201 1 
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1 1)  Mani  RS,  Iyer  M,  Cao  Q,  Brenner  CJ,  Wang  L,  Ghosh  A,  Nyati  MK,  Wu  YM,  Cao  X,  Lonigro 
RJ,  Qin  Z,  Tomlins  SA,  Palanisamy  N,  Varambally  S,  Chinnaiyan  AM.  Functional  consequences 
of  ETS  gene  fusions  in  prostate  cancer.  5th  Annual  Multi-Institutional  SPORE  Prostate  Cancer 
Program  Retreat,  Fort  Lauderdale,  FL.  Mar  19-21,  2012 

Award: 

2011  American  Association  for  Cancer  Research  102nd  Annual  Meeting  AACR- Aflac  Scholar-in- 
Training  Award 

Dataset 

Gene  expression  data  from  DU145  and  BT-549  overexpressing  miRNAs  were  deposited  into  GEO 
(GSE26996). 


Conclusion: 

Our  previous  studies  showed  that  microRNA  miR-101  represses  histone  methyltransferase  EZH2  at 
transcript  and  protein  levels.  miR-101  expression  levels  were  decreased  during  prostate  cancer 
progression  and  negatively  correlated  with  EZH2.  In  addition,  we  demonstrated  that  genomic  loci 
encoding  miR-101  were  deleted  in  metastatic  prostate  cancer,  and  the  genomic  loss  of  miR-101  leads  to 
EZH2  overexpression  in  cancer. 

Based  on  the  aims  of  this  proposal,  we  performed  miRNA  profiling  and  identified  additional  miR-101 
downstream  targets,  BMI1  and  RING2,  through  EZH2-regulated  miRNAs  miR-181a,b  miR-200b,c  and 
miR-203.  Importantly,  we  determined  these  miRNAs,  including  miR-101,  are  inversely  correlated  with 
PRC  proteins  EZH2,  BMI1  and  RING2  during  prostate  cancer  progression.  Additionally  we  revealed 
that  the  coordinated  regulation  between  PRC2  and  PRC1  is  mediated  through  miRNAs.  Furthermore,  we 
found  that  EZH2  represses  miR-101  by  binding  to  its  promoter  regions.  By  RNA-seq,  we  did  not  find 
miR-101  binding  site  mutations  in  EZH2  3’UTR  regions. 

Finally,  we  identified  that  another  potential  prostate  cancer  marker  STMN 1  is  positively  correlated  with 
EZH2  during  prostate  cancer  progression.  Interestingly,  miR-101  binds  to  the  STMN1  3’UTR  and 
represses  it  at  protein  and  transcript  levels. 
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ABBREVIATIONS 


PC  A:  prostate  cancer 

miRNA:  microRNA 

EZH2:  enhancer  of  zeste  homolog  2 

BMI1:  B  lymphoma  Mo-MLV  insertion  region  1  homolog 

PRC:  Polycomb  Repressive  Complex 

STMN1:  Stathmin  1 
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EZH2  is  a  member  of  Polycomb  Repressive  Complex  2  (PRC2)  that 


contains  histone  methyl  transferase  activity,  and  is  involved  in  epigenetic  regulation 
of  target  genes  in  normal  cells  as  well  as  cancer.  Along  with  other  components  of 
PRC2,  EZH2  trimethylates  histone  H3K27  causing  chromatin  condensation  and 
gene  silencing.  In  this  review,  we  discuss  the  role  of  Polycomb  Repressive  Complexes 
in  cellular  processes  and  cancer,  and  the  potential  for  EZH2  to  serve  as  a  cancer 
biomarker  and  a  therapeutic  target. 


Polycomb  Group  Protein  Complexes 

During  development,  cell  fate  is  determined  by  the  initial  gene  expression 
pattern  in  a  cell,  after  which  cellular  identity  must  be  maintained  epigenetically  by 
establishing  regions  of  active  and  silenced  chromatin  The  Polycomb  Group  (PcG)  and 
Trithorax  Group  (TrxG)  proteins  are  critical  for  maintaining  the  cellular  memory  and  cell 
type-specific  transcription  programs2,  3  by  enacting  histone  modification,  DNA 
methylation  and  chromatin  remodeling  4~10.  In  mammals,  PcG  gene  expression  exhibits  a 
tissue-specific,  temporal  pattern  that  regulates  cell  fate11  . 

The  genes  encoding  the  Polycomb  group  and  Trithorax  group  are  highly 
conserved  between  Drosophila  and  mammals  and  regulate  other  genes  involved  in 
development  and  cell  cycle  control.  TrxG  and  PcG  proteins  were  initially  identified  in 
Drosophila  by  genetic  screening  and  shown  to  establish  and  maintain  homeobox  (Hox) 
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gene  expression  patterns11'13.  However,  TrxG  and  PcG  targets  are  not  limited  to  Hox 
genes,  as  they  are  reported  to  bind  to  thousands  of  other  chromosomal  sites  u’ 13, 14 .  Both 
TrxG  and  PcG  form  multi-protein  complexes  that  bind  to  specialized  DNA  elements 
called  Polycomb/Trithorax  Response  Elements  [PREs/TREs]  present  in  the  regulatory 
regions  of  many  genes15’ 16 . 

Two  distinct  Polycomb  complexes  have  been  characterized  by 
immunoprecipitation,  yeast  two-hybrid,  and  size-fractionation  experiments  in  mammalian 
systems  (Figure  1  and  Table  1).  Polycomb  Repressive  Complex  1  (PRC1)  includes  BMI1, 
RING1,  PHI/2/3,  and  PCI/2/3  (Psc,  dRING,  Ph  and  Pc  respectively  in  Drosophila )  6’7’ 17; 
Polycomb  Repressive  Complex  2  (PRC2)  includes  EZH2,  EED,  SUZ12,  RbAp46/48  and 
AEBP2  (E(z),  Esc,  Su(z)12,  RbAp48  and  Jing  respectively  in  Drosophila )  6’7’  n’ 17.  EZH2 
is  histone  methyltransferase  (HMT),  the  core  enzyme  of  PRC2.  Several  studies  have 
suggested  that  HDACs  are  associated  with  PRC2  through  interaction  with  EED18’  19. 
Supporting  this  notion  is  the  fact  that  HD  AC  inhibitors  such  as  trichostatin  A  and  SAHA 
can  inhibit  EZH2/PRC2  function20. 

Based  on  recent  discoveries,  a  model  of  the  molecular  mechanism  of  Polycomb 
group  complex  function  has  been  proposed.  First,  histone  acetyltransferases  mark 
transcriptionally  active  DNA  regions  by  acetylating  histone  tails.  Specific  molecular  cues 
are  then  transmitted  and  PRC2  binds  to  the  PRE.  HDACs  that  deacetylate  the  histone 
tails  are  then  recruited  followed  by  the  methylation  of  the  histone  tails  by  PRC2.  The 
resulting  remodeled  chromatin  renders  the  promoter  regions  of  target  genes  accessible  for 
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further  histone  tail  deacetylation  followed  by  methylation.  The  methylation  marks 


function  as  a  binding  site  for  the  N-terminal  chromodomain  of  Pc  proteins  that  then 
recruits  PRC1  to  the  promoters  of  target  genes  resulting  in  the  inhibition  of  transcription 
initiation  4’ 21  ’ 22  .  Collectively  these  experiments  provide  the  hierarchical  model  of  how 
PRC2  and  PRC1  complexes  are  assembled  and  the  molecular  events  in  the  repression  of 
gene  expression  (Figure  1).  Accordingly,  PRC2  is  also  referred  to  as  PRCi  (initiation), 
and  PRCI  as  PRCm  (maintenance). 

In  mammals  there  are  two  homologues  of  E(z),  EZH1  and  EZH2.  Sequence 
analysis  shows  that  EZH1/2  proteins  have  four  conserved  domains  aligning  to  E(z)  23’24: 
homologue  domain  I  (HI  domain),  homologue  domain  II  (H2  domain),  cysteine-rich 
domain  and  C-terminal  SET  (SU(var)3-9,  E(z)  and  Trithorax)  domain.  The  histone 
methyltransferase  activity  of  this  family  of  enzymes  is  contained  within  the  SET  domain 
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Using  yeast  two-hybrid  screen,  human  Esc  homolog  EED  was  found  to  interact 
with  EZH2  in  vitro ;  using  co-immunoprecipitation,  EZH2  and  EED  were  shown  to  also 
form  a  complex  in  vivo.  Several  studies  demonstrated  that  the  WD40  domains  of  EED  are 
essential  for  the  EZH2-EED  interaction,  which  can  be  disrupted  by  specific  point 
mutations  in  the  WD40  domains7’ 23-25 . 

Another  PcG  protein,  SUZ12,  is  an  essential  component  of  PRC2  and  is 
required  for  its  histone  methyltransferase  (HMTase)  activity.  Isolation  and 
characterization  of  the  enzymatic  complex  revealed  that  SUZ12,  RbAp48,  and  AEBP2 
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interact  with  the  EZH2-EED  complex7  and  the  fully  reconstituted  complex  can 
specifically  methylate  Histone  H3  in  vitro.  In  this  complex,  EZH2,  EED,  and  SUZ12  are 
all  required  for  the  HMTase  activity.  RbAp48  facilitates  SUZ12’s  interaction  with  EZH2 
and  EED.  AEBP2  significantly  enhances  the  EZH2-EED-SUZ12  HMTase  activity7 . 

Another  PcG  protein,  polycomblike  (Pel)  or  PHF1,  was  identified  and  found  to 
interact  with  PRC2  in  Drosophila  and  mammals26’ 21 .  Using  a  tandem  affinity  purification 
strategy,  Nekrasov  et  al.  characterized  a  novel  PRC2  complex  containing  Pel  from  a 
transgenic  Drosophila  strain  expressing  the  tagged  Pel  fusion  protein  .  Like  the 
previously  characterized  PRC2  complex,  the  Pcl-PRC2  complex  is  an  H3K27-specific 
HMTase.  In  the  Pel'7"  strains,  H3K27me3  levels  were  much  lower  compared  to  wild-type 
strains;  however,  H3K27mel  and  H3K27me2  levels  were  higher,  suggesting  that 
Pcl/PHFl  is  required  for  generating  high  levels  of  tri-methylated  H3K27  to  maintain  the 
Polycomb-repressed  chromatin  structure. 

JARID2  (Jumonji,  AT  Rich  Interactive  Domain  2)  was  reported  to  interact  with 
PRC2  in  embryonic  stem  cells  and  cancer  cells  by  five  independent  groups28"32  . 
Although  it  is  not  clear  whether  JARID2  has  histone  demethylase  activity,  JARID2  can 
bind  to  GC  and  GA  rich  motifs  in  DNA  directly  through  its  C-terminal  domain. 
Interestingly,  JARID2  was  shown  to  modulate  PRC2  activity  by  promoting  its 
recruitment  to  target  genes  while  simultaneously  attenuating  its  methyltransferase  activity. 

While  the  PRC2  complex  possesses  HMTase  activity,  the  PRC1  complex 
harbors  H2AK119  ubiquitin  E3  ligase  activity4,21’22  .  The  mammalian  PRC1  H2AK119 
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ubiquitin  E3  ligase  complex  consists  of  several  PcG  proteins.  RING2/RING1B  is  the 


catalytic  subunit,  but  the  presence  of  BMI1,  RING1/RING1A,  and  PC3  significantly 
augments  the  ubiquitin  E3  ligase  activity.  Notably,  BMI1  interacts  with  RING1/RING1A, 
RING2/RING1B/RNF2,  PH2,  and  PC3,  indicating  that  BMI1  is  important  for  the 
integrity  of  this  complex.  Although  another  PcG  protein,  Mel-18,  can  substitute  for  BMI1 
in  maintaining  the  PRC1  complex,  Mel- 18  cannot  stimulate  the  ubiquitin  E3  ligase 
activity  of  PRC1. 

Another  PRC1  protein  Pc2  is  reported  to  have  SUMO  E3  ligase  activity33  . 
Kagey  et  al.  demonstrated  that  Pc2  can  catalyze  the  sumoylation  of  the  transcriptional 
repressors  CtBPl  and  CtBP2.  Pc2  can  interact  with  CtBPl  and  CtBP2  via  a  PLDLS-like 
motif  and  recruit  them  to  PRC1  34 .  Pc2  does  not  exhibit  obvious  homology  to  any  other 
known  E3s.  In  2005,  the  same  group  found  that  the  C-terminal  part  of  Pc2  can  act  as  a 
docking  site  for  UBC9  35 ,  a  well-known  E2  for  sumoylation  and  their  substrates,  CtBPl 
and  CtBP2. 

EZH2  and  PcG  in  cellular  processes 

In  the  Drosophila  embryo  the  segmental  body  plan  is  established  within  the  first 
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three  and  a  half  hours  of  development  .  Each  segment’s  developmental  fate  is 
determined  by  a  specific  combination  of  homeotic  gene  products,  their  concentrations, 
and  temporal  expression.  Furthermore,  maintenance  of  segment  fate  decision  is  carried 
out  by  another  unique  set  of  homeotic  proteins.  The  initial  expression  pattern  in  each 
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segment  is  established  by  transcription  factors  encoded  by  segmentation  genes  that  decay 


within  five  to  seven  hours  of  development,  after  which  the  Polycomb  (PcG)  and 
Trithorax  (TrxG)  groups  assume  responsibility  for  maintaining  the  expression  pattern. 
PcG  and  TrxG  proteins  have  been  demonstrated  to  maintain  long-term  gene  silencing 
during  development  in  Drosophila.  However,  E(z.)  and  Esc  differ  from  other  PcG  proteins 
in  that  E(z)  and  Esc  are  required  early  on  in  development  whereas  other  PcG  proteins 
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appear  to  be  required  relatively  later  .  Furthermore,  E(z)  and  Esc  are  the  most  highly 
conserved  PcG  genes  throughout  evolution  as  they  are  the  only  two  PcG  genes  found  in 
C.  elegans  genome. 

Several  groups  utilized  a  global  approach  to  identify  target  loci  of  PcG  proteins 
by  genome-wide  mapping  of  PcG  binding  in  human  fibroblasts,  human  ES  cells,  mouse 
ES  cells,  and  Drosophila  37-41  .  All  of  the  studies  demonstrated  that  PcG  proteins  directly 
repress  many  genes  involved  in  early  development,  a  wide  variety  of  developmental 
processes  and  cell  differentiation.  In  hematopoietic  stem  cells  (HSCs),  PRC2  and  PRC1 
regulate  the  self-renewal  and  differentiation  of  HSCs  42  . 

PcG  proteins  also  regulate  genes  in  many  signaling  pathways  including  WNT  43, 
FGF  (Fibroblast  Growth  Factor)  44,  BMP  (Bone  Morphogenic  Protein)  45  and  TGFp 
(Transforming  Growth  Factor  P)  46  that  are  required  for  linage  differentiation  and 
associated  with  cancer  progression.  Therefore  PcG  studies  in  ES  cells  can  be  informative 
in  understanding  how  dysregulation  of  these  genes  can  promote  tumorigenesis. 
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PRC2  and  trxG  group  in  cancer 


Recently,  Wilson  et  al.  reported  that  the  loss  of  the  tumor  suppressor  SNF5,  a 
TrxG  protein,  led  to  the  overexpression  of  EZH247.  They  also  found  that  Polycomb  target 
genes  are  marked  by  H3K27me3  and  transcriptionally  repressed  in  5A^F5-negative 
cancers.  SNF5  is  a  component  of  the  TrxG  SWI/SNF  complex  thought  to  regulate  gene 
expression  through  perturbations  in  nucleosome  positioning.  Tumor  formation  via  loss  of 
SNF5  can  be  inhibited  by  EZH2  inactivation  in  vivo.  SNF5  binds  directly  to  the  EZH2 
promoter  and  the  SWI/SWF  complex  negatively  regulates  EZH2  expression. 

Disruption  of  EZH2  in  cancer 

Every  year,  over  230,000  American  women  are  diagnosed  with  breast  cancer, 
and  over  240,000  American  men  are  diagnosed  with  prostate  cancer  48.  Of  these,  15%  of 
breast  cancer  and  11%  of  prostate  cancer  cases  will  result  in  death  48.  Given  the  high 
mortality  rate,  there  is  an  urgent  need  to  identify  robust,  specific  prognostic  markers  of 
aggressive  prostate  and  breast  cancer  so  early  therapeutic  interventions  can  be  pursued. 

In  epithelial  derived  tumors,  EZH2  was  first  observed  to  be  significantly 
associated  with  metastatic  prostate  cancer.  Along  with  MTA1,  HPN,  PIM1  and  several 
other  genes,  EZH2  is  overexpressed  in  castration-resistant  metastatic  prostate  cancer  20, 49  . 
These  data  imply  that  EZH2  may  be  a  potent  biomarker  of  prostate  metastases  20 .  Higher 
EZH2  protein  level  is  associated  with  a  shorter  disease-free  interval  after  initial  surgical 
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treatment,  lower  overall  survival,  and  a  high  probability  of  cancer-specific  death20, 50. 
Since  the  initial  discovery  of  EZH2  overexpression,  EZH2  has  also  been  shown  to  be 
involved  in  a  number  of  other  cancers  such  as  bladder  51-54  ,  gastric  55  ,  lung  56  ,  liver 
cancer 57  and  leukemia  58, 59  . 

EZH2  is  aberrantly  expressed  due  to  a  variety  of  mechanisms  that  all  lead  to 
increased  PRC2  activity.  For  example,  point  mutations  within  the  EZH2  SET  domain 
(Tyr641)  have  been  identified  in  diffuse  large  cell  B-cell  lymphomas  of  germinal  center 
origin  60  and  myeloid  neoplasms  61 ' 62 .  Intriguingly,  compared  to  the  wild-type  enzyme, 
disease-associated  Tyr641  mutations  could  not  mono-methylate  H3K27  but  enhanced 
catalytic  efficiency  for  the  di-  and  trimethylation  activities,  implying  that  both  WT  and 
variants  of  EZH2  may  be  involved  in  neoplastic  transformation  following  this  specific 

/TO 

EZH2  mutation  .  Recently,  in  addition  to  EZH2  mutations,  high  frequency  mutations 
have  been  discovered  in  the  PRC2  proteins  EED  and  SUZ12  in  T-cell  acute 
lymphoblastic  leukemia  (T-ALL)  64-66  . 

Li  et  al.  identified  a  series  of  DNA  rearrangements  within  PRC2  components  in 
endometrial  stromal  tumors  .  The  authors  discovered  a  (7;  17)  translocation  resulting  in 
the  fusion  between  C-terminal  sequences  of  SUZ12  to  N-terminal  sequences  of  JAZF1  in 
three  primary  endometrial  stromal  sarcomas  (ESSs)  cell  lines.  Importantly,  this  fusion 
protein  can  augment  EZH2  activity  and  increase  H3K27  methylation.  Hence,  gene 
fusions  represent  another  mechanism  by  which  PRC2  can  promote  carcinogenesis. 

While  EZH2  and  PRC2  components  can  be  genetically  disrupted  in  cancer, 


9 


EZH2  overexpression  can  also  be  caused  by  cancer-associated  signaling  pathways.  In 


2003,  Bracken  et  al.  demonstrated  that  the  E2F-RB  axis  can  regulate  EZH2  and  PRC2  in 

/TO 

cancer  cells  .  Here,  it  was  shown  that  activation  of  E2F1,  E2F2  and  E2F3  can  markedly 
increase  EZH2,  EED  and  SUZ12  expression  levels  by  directly  binding  to  their  promoter 
regions.  In  contrast,  overexpression  of  the  upstream  tumor  suppressors  of  this  pathway, 
pRB  and  p!6,  repress  the  activity  of  E2Fs  and  decrease  expression  of  EZH2  and  EED. 
Consistent  with  the  findings  above,  EZH2  and  EED  are  increased  in  pRb_/"  mouse 

/TO 

embryonic  fibrobroblasts  (MEFs)  ;  loss  of  pl6  leads  to  overexpression  of  PRC2 
expression  at  both  transcript  and  protein  levels  and  induces  the  DNA  hypermethylation  of 
PRC2  targets  69 . 

In  addition  to  direct  regulation  by  the  E2F-RB  signaling  axis,  EZH2  was 
reported  to  be  downregulated  in  senescent  cells  as  a  result  of  p53  activation70.  This 
downregulation  is  dependent  on  an  inhibitor  of  the  E2F  pathway,  p21/wafl.  Interestingly, 
BMI1  is  reportedly  a  target  of  E2F1-3  but  not  of  E2F4  or  E2F5,  suggesting  that  PRC1 
activity  may  also  be  regulated  by  this  pathway71.  Figure  2  proposes  a  model  of  how  PcG 
proteins  are  coordinately  regulated  in  cancer. 

The  functional  role  of  EZH2  in  cancer 

To  specifically  investigate  the  role  of  PRC2  in  metastatic  prostate  cancer,  Yu  et 
al.  performed  ChIP-on-chip  with  SUZ12  and  H3K27me3  in  late-stage,  aggressive 
prostate  cancer  tissues72.  Genome-wide  mapping  showed  that  there  was  a  strong  overlap 
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between  the  genomic  sites  occupied  by  SUZ12  and  H3K27me3  in  metastatic  prostate 


cancer  and  prostate  cancer  cell  lines.  Further,  Oncomine  Molecular  Concept  Map  (MCM) 
analysis  -  a  tool  for  comparing  gene  lists  from  high  throughput  profiling  studies  - 
showed  that  H3K27me3-occupied  genes  are  consistently  downregulated  in  prostate, 
breast  and  lung  cancers.  Importantly,  there  was  a  strong  link  between  H3K27me3- 
occupied  genes  in  prostate  metastases  and  H3K27me3-,  SUZ12-,  or  EED-occupied  genes 
in  embryonic  stem  cells,  supporting  the  overlapping  role  of  PRC2  in  regulating  stem  cell 
pluripotency  and  differentiation  and  prostate  cancer  progression.  Importantly,  the 
Polycomb  repression  signature  successfully  predicted  the  clinical  outcome  of  prostate  and 
breast  cancer  patients. 

Many  of  the  gene  repressed  by  EZH2  are  considered  tumor  suppressors 
including  ADRB2  73  ,  CDH1  74  ,  CIITA  75  ,  RUNX3  76  ,  SLIT2  11  ,  raplGAP  78  among 
others  in  various  cancers.  Besides  tumor  suppressor  genes  EZH2  also  represses  genes 
associated  with  DNA  repair  79, 80  .  In  mammals,  hampered  double-strand  break  (DSB) 
repair  processes  can  lead  to  chromosomal  instability  resulting  in  cell  death  or  cancer. 
Recently  Chang  et  al.  found  that  RAD51  downregulation  by  EZH2  leads  to  impaired 
DNA  damage  repair  and  increase  in  oncogenic  hits  that  promote  breast  tumor  initiating 
cells  (BTICs)79.  Likewise,  high  levels  of  EZH2  are  also  associated  with  decreased  nuclear 
levels  of  BRCA1,  a  tumor  suppressor  protein  present  in  invasive  breast  carcinomas81. 
BRCA1  regulates  DNA  repair,  activates  cell-cycle  checkpoints  and  has  a  vital  role  in  the 
maintenance  of  chromosomal  stability.  Gonzalez  et  al.  reported  that  EZH2  alters  the 
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intracellular  location  of  BRCA1  protein  in  benign  and  ER-negative  breast  cancer  cells81. 

Consistent  with  a  role  in  regulating  the  expression  of  many  cancer-related  genes, 
EZH2  knockdown  by  RNA  interference  (RNAi)  in  prostate  cell  lines  PC3  and  RWPE  as 
well  as  the  osteosarcoma  cell  line  U20S  resulted  in  attenuated  cell  proliferation,  but  did 
not  induce  apoptosis  20’ 68.  Additionally,  EZH2-depleted  cells  displayed  cell-cycle  arrest 
in  the  G2/M  phase  indicating  that  EZH2  inhibits  cell  proliferation  by  mitigating  G2/M 
transition82.  Min  et  al.  discovered  that  EZH2  could  promote  cancer  metastasis  by 
inducing  K-ras  and  NF-kB  pathways  by  inhibiting  DAB2IP,  a  Ras  GTPase-activating 
protein  (GAP)  in  prostate  cancer83.  EZH2  overexpression  has  also  been  shown  to  trigger 
the  activation  of  AKT-1  cascade,  suggesting  that  EZH2  may  promote  genomic  stability 
through  the  PI3K/AKT-1  pathway79.  For  example,  EZH2-mediated  RAF1-  P-catenin 
activation  may  play  a  role  in  aggressive  BTIC  proliferation.  Targeting  the 
EZH2/Polycomb  complex  and  its  downstream  activation  pathways  such  as  RAF-ERK 
signaling  by  AZD6244  (also  known  as  Selumetinib),  a  potent,  selective,  and  ATP- 
uncompetitive  inhibitor  of  MAPK/ERK  kinase  1/2,  may  potentially  block  BTICs  and 
prevent  breast  cancer  recurrence. 

Beyond  its  role  in  epigenetic  regulation,  EZH2  protein  may  have  cytoplasmic 
functions84'86  .  In  late  stage  cancer  cell  lines,  for  example,  mutant  EZH2  protein  lacking 
the  NLS  continues  to  form  a  complex  with  other  PRC2  components.  This  cytosolic 
complex  retains  HMTase  activity  and  promotes  TCR-  and  PDGF-induced  actin 
polymerization.  It  is  still  unclear  if  the  aforementioned  signaling  pathways  are  dependent 
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on  EZH2  nuclear  and/or  cytoplasmic  pathways,  and  understanding  these  differences  may 


be  critical  for  therapeutic  design  in  the  future. 

Phosphorylation  of  EZH2 

EZH2  is  regulated  through  direct  phosphorylation.  In  2005  Cha  et  al.  identified 
EZH2  S21  as  a  downstream  target  of  PI3K-Akt  signaling  pathway87.  In  this  study, 
phosphorylation  of  EZH2  at  S21  did  not  affect  the  PRC2  complex  formation,  but  reduced 
the  affinity  of  EZH2  for  H3K27  and  caused  diminished  HMTase  activity87.  Recent 
studies  have  also  demonstrated  that  EZH2  threonine  350/345  and  487  can  be 
phosphorylated  by  CDK1  and  CDK288’ 89.  These  studies  demonstrated  that  Thr350/345- 
phosphorylated  EZH2  is  responsible  for  marking  target  genes  and  to  prime  them  for 
recruitment  of  unphosphorylated  PRC2,  leading  to  an  increase  in  H3K27me3  and 
subsequent  repression  of  target  genes  89.  In  contrast,  Thr487  phosphorylation  by  CDK1 
indirectly  hinders  EZH2  methyltransferase  activity  and  disrupts  its  interactions  with 
SUZ12  and  EED90.  It  has  been  suggested  that  a  combination  of  EZH2  phosphorylation 
events  determines  the  distribution  of  PRC2  along  chromatin  during  G2/M.  The  sequential 
phosphorylation  of  Thr350/345  and  Thr487  also  regulates  enzymatic  activity  and  PRC2 
association/dissociation  from  target  genes88’ 90. 

EZH2  and  microRNAs 

microRNAs  are  a  class  of  negative  regulatory  small  non-coding  RNAs  that  play 
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a  role  in  stem  cell  maintenance  and  cancer  progression.  In  2008,  Varambally  et  al. 


discovered  that  microRNA-101  (miR-101)  repressed  EZH2  at  transcript  and  protein 
levels  by  targeting  its  3’UTR  91.  Consistent  with  this  observation,  miR-101  and  EZH2  are 
negatively  correlated  during  prostate  cancer  progression.  In  contrast  to  EZH2,  miR-101  is 
highest  in  benign  tissue,  but  lowest  in  metastatic  prostate  cancer.  By  genomic  DNA 
qPCR  and  aCGH,  the  authors  found  that  the  30-50%  of  metastatic  prostate  cancer 
exhibited  genomic  loss  of  miR-101  loci  leading  to  the  overexpression  of  EZH2  in 
prostate  cancer91.  Additionally,  the  genomic  loss  of  miR-101  loci  was  observed  in  breast 
cancer,  gastric  cancer,  gliobastoma  multiforme,  lung  adenocarcinoma,  and  acute 
lymphocytic  leukemia91.  Subsequently,  miR-101  was  found  to  represses  EZH2  in  bladder 
cancer  92 ,  squamous  cell  carcinoma 78 ,  and  nasopharyngeal  carcinoma  93. 

Other  microRNAs  have  also  been  reported  to  regulate  EZH2  and  other  PcG 
proteins  in  cancer  (Table  2).  For  example,  miR-26a  was  first  reported  to  repress  EZH2 
during  myogenesis94  .  Later,  several  groups  discovered  that  c-Myc  repressed  miR-26a 
expression  resulting  in  EZH2  overexpression  in  lymphoma  95,  nasopharyngeal  carcinoma 
93,96  and  rhabdomyosarcoma  97 .  miR-98  and  miR-214  are  also  linked  to  EZH2  regulation 
in  cancer  93  and  embryonic  stem  cells  98.  It  has  also  been  reported  that  miR-183,  miR- 
200b,  miR-200c  and  miR-203  repress  BMI1  in  breast,  pancreatic  and  colorectal  cancer99, 
100  and  are  linked  to  epithelial  to  mesenchymal  transition  10°. 

Coordinated  regulation  of  PRC2  and  PRC1  in  cancer  by  microRNAs 


14 


EZH2  is  highly  overexpressed  in  prostate  carcinoma  metastasis  precursor  cells 
along  with  PRC1  protein  BMI1101.  Utilizing  quantitative  immunofluorescence  co¬ 
localization  analysis,  Berezovska  et  al.  demonstrated  a  marked  enrichment  of  a 
population  of  circulating  human  prostate  carcinoma  metastasis  precursor  cells  with  dual¬ 
positive  high-BMIl/EZH2-expressing  cells.  High  levels  of  BMI1  and  EZH2  allow 
prostate  cancer  cells  to  resist  apoptosis  induction  in  cells  of  epithelial  origin  in  response 
to  attachment  deprivation.  Depletion  of  BMI1  or  EZH2  in  prostate  carcinoma  metastasis 
precursor  cells  results  in  attenuated  tumorigenic,  metastatic  and  proliferative  potential 
when  injected  into  mice. 

Our  group  recently  reported  that  coordinated  regulation  of  PRC2  and  PRC1 
expression  is  mediated  by  specific  microRNAs  including  miR-203,  miR-181a,b  and  miR- 
200b, c  102  .  The  miRNA  levels  were  decreased  and  PRC1  proteins,  BMI1,  RING2,  and 
ubiquityl-H2A  were  elevated  in  metastatic  prostate  cancer;  these  PRC2-regulated 
miRNAs  repressed  PRC1  proteins,  BMI1  and  RING2,  as  well  as  H2AK119 
ubiquitination.  PRC2  epigenetically  represses  miR-181a,  b,  miR-200b,c  and  miR-203 
expression  by  enacting  H3K27me3  trimethylation  at  these  loci.  Exogenous 
overexpression  of  miR-181a,  b,  miR-200b,  c  and  miR-203  inhibit  tumorigenic  phenotype 
in  vitro.  We  also  showed  that  miR-203,  miR-200b  and  miR-181b  overexpression 
suppressed  prostate  tumor  formation  and  growth  in  mouse  xenografts.  Taken  together, 
these  data  reveal  that  PRC2,  PRC1  and  microRNAs  form  an  integral  regulatory  axis  in 
the  epigenetic  regulation  of  gene  expression. 
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EZH2  and  non-coding  RNAs 


Non-coding  RNAs  -  RNAs  without  an  extended  open  reading  frame  -  were 
recently  shown  to  display  abnormal  expression  patterns  in  many  cancers  103’  104  and 
expression  of  antisense  ncRNAs  may  lead  to  silencing  of  respective  tumor  suppressor 
genes  triggering  carcinogenesis.  ncRNAs  may  also  play  a  role  in  recruiting  factors 
mediating  chromatin  repression.  HOTAIR  is  known  to  repress  transcription  across  40  kb 
of  the  HOXD  locus  by  altering  the  chromatin  trimethylation  state105.  The  genes  in  the 
human  HOXD  locus  are  regulated  in  trans  by  HOTAIR  and  produced  by  the  HOXC  locus. 
Intriguingly,  EZH2,  SUZ12,  and  EED  proteins  possess  potential  binding  sites  for 
HOTAIR105.  Although  HOTAIR  works  in  trans,  evidence  from  antisense  transcript 
studies  indicate  that  many  ncRNAs  can  function  by  directing  PRC2  in  cis  106  .  Taken 
together,  these  results  show  that  long  ncRNAs  are  ideal  for  targeting  chromatin 
remodeling  complexes  to  specific  regions  of  the  genome  106  108  ,  utilizing  both  cis  and 
trans  mechanisms. 

Most  recently,  Prensner  et  al.  reported  that  there  is  a  significant  group  of 
unannotated  ncRNAs  with  cancer-specific  expression  patterns  104.  They  nominated  121 
such  Prostate  Cancer  Associated  Transcripts  (PCATs)  that  collectively  functioned  as  well 
as  known  biomarkers  in  differentiating  benign,  cancerous,  and  metastatic  prostate  cancer. 
In  particular,  they  identified  an  uncharacterized  prostate-specific  ncRNA  called  PCAT-1. 
PCAT-1  is  negatively  regulated  by  the  PRC2  in  vitro  and  is  overexpressed  in  a  subset  of 
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prostate  cancers,  particularly  in  metastatic  tissue.  Further,  patterns  of  PCAT-1  and  PRC2 


expression  stratified  patient  tissues  into  molecular  subtypes  distinguished  by  expression 
signatures  of  PCAT-1 -repressed  target  genes  including  the  putative  target  genes  BRCA2 
and  CENPF.  Consistent  with  the  role  of  PCAT-1  and  the  PRC2  in  advanced  prostate 
cancer,  PCAT-1  was  shown  to  promote  cell  proliferation  in  prostate  cancer  cell  lines. 
Taken  together,  these  data  add  to  the  growing  importance  of  unannotated  ncRNA  species 
in  oncogenesis  and  suggest  additional  mechanisms  by  which  PRC2  can  regulate  cancer 
progression. 

EZH2  as  a  target  for  cancer  therapy 

Given  the  role  of  EZH2  in  a  wide  range  of  tumorigenic  processes  and  pathways, 
EZH2  is  an  attractive  therapeutic  target  for  targeted  therapy  in  cancer.  As  discussed  above, 
many  groups  have  now  demonstrated  that  knockdown  of  EZH2  inhibits  cancer  cell 
growth,  motility,  invasion  potential  and  tumorigenesis.  Accordingly,  several  different 
approaches  to  inhibit  EZH2  function  have  been  explored. 

One  approach  involves  the  use  of  HDAC  inhibitors  such  as  TSA  and  SAHA 
that  are  known  to  attenuate  EZH2-induced  invasion.  However,  HDAC  inhibitors  have 
pleiotropic  effects  and  regulate  EZH2  indirectly. 

In  order  to  directly  target  EZH2,  many  groups  are  currently  screening  small 
molecule  inhibitors  of  EZH2  using  various  approaches109, 110.  Tan  et  al.  discovered  that  S- 
adenosylhomocysteine  hydrolase  inhibitor  3-Deazaneplanocin  A  (DZNep)  can  deplete 
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EZH2,  EED  and  SUZ12  protein  levels  and  inhibit  H3K27  methylation110.  PRC2 
depletion  by  DZNep  allows  RNA  polymerase  II  complex  to  occupy  the  promoters  of 
PRC2  targets  and  initiate  transcription.  Interestingly,  DZNep  can  re-activate  some 
hypermethylated  genes  while  5’-Aza-dC  and  TSA  fail  to,  suggesting  that  DZNep  could 
trigger  DNA  demethylation.  DZNep  was  also  shown  to  induce  apoptosis  in  cancer  cells 
with  high  EZH2  expression  but  not  in  normal  cells  1U.  Thus  DZNep  may  have  promising 
clinical  potential  for  cancer  therapy. 

Besides  small  molecule  inhibitors  of  EZH2,  additional  therapeutic  strategies  are 
being  developed.  Using  synthesized  fragment  peptides  of  EZH2,  peripheral  blood 
mononuclear  cells  (PBMCs)  were  stimulated  to  produce  EZH2-specific  cytotoxic  T 
lymphocytes  (CTLs).  The  EZH2-specific  CTLs  can  generate  IgG  against  EZH2  and  are 
toxic  to  HLA-A24  positive  cells  that  have  high  EZH2  expression112.  This  finding 
provides  another  possible  mechanism  by  which  EZH2  function  can  be  inhibited  during 
cancer  progression. 

Taken  together,  our  understanding  of  EZH2  disruption,  regulation  and  function 
has  increased  immensely  in  the  last  few  years.  However,  many  questions  remain  about 
the  mechanisms  by  which  EZH2  regulates  cancer  progression,  especially  in  different 
cellular  contexts.  In  the  future,  targeted  therapies  can  hopefully  be  developed  to  both 
indirectly  and  directly  inhibit  EZH2  and  thus  PRC2  activity  in  cancer.  In  this  light, 
biomarker-driven  clinical  trials  can  be  designed  to  specifically  target  the  large  subset  of 
cancers  with  aberrant  EZH2  activity. 
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Table  1  Main  components  of  Polycomb  group  proteins  and  their  functions  and  roles  in  cancer 


Drosophila 

proteins 

Human 

homologues 

Complex 

Domains 

Biochemical  activities 

Role  in  cancer 

Pc 

HPC1 

HPC2 

HPC3 

PRC1 

Chromodomain 

Recognition  of 
methylated  histones; 
SUMO  E3  ligase 
(HPC2) 33’35 

Ph 

HPH1 

HPH2 

HPH3 

PRC1 

Zinc-finger 

SPM  domain 

Stabilization  of  PRC1 
complex  113 

Psc 

BMI1 

MEL  18 

PRC1 

Zinc-finger 

domain 

Integrity  of  PRC1 
complex,  stimulation  of 
ubiquitin  E3  ligase 
activity  4’ 113 

Target  of  E2F1-3  71 ; 

Immortalize  normal  human  oral 
keratinocytes  114; 

Promote  the  generation  of  B-  and  T-cell 
lymphomas  with  c-Myc  115; 

Overexpressed  in  B  cell  non-Hodgkin 
lymphoma116,  breast  cancer  m,  colorectal 
carcinoma  118,  Hodgkin’s  lymphoma  119, 
liver  carcinoma  ,  medulloblastoma 
non-small  cell  lung  cancer  122  and  prostate 

123 

cancer 

Ring 

RING1/RING1 A 
RING2/RING1B 

PRC1 

Zinc-finger 

domain 

Ubiquitin  E3 

ligase(RINGlB) 4’ 21  ’ 22 

Overexpressed  in  Hodgkin’s  lymphoma 
and  prostate  cancer  123; 

Increase  c-Jun  and  c-fos  ; 

Induce  anchorage-independent  growth  and 
enhance  tumor  formation  in  nude  mice  124 

Pho 

YY1 

PRC1 

PRC2 

Zinc-finger 

domain 

DNA  binding; 

Interact  with  E2Fs  and 
Mdm2 

Cooperate  with  c-myc  and  E2Fs  ; 

Inhibit  the  activation  of  p53  126  127; 

Enhance  cyclooxygenase-2  gene 

20 


expression  in  macrophages  ; 
Overexpressed  in  osteosarcoma  129  and 
prostate  cancer  130 

E(z) 

EZH1 

EZH2 

PRC2 

SET  domain 

HI,  H2 

protein-protein 

interaction 

domain; 

Cysteine-rich 

domain 

Histone 

methyl  transferase ; 

Target  ofE2Fl-3  68 ; 

Promote  anchorage-independent  growth 
and  invasion  ; 

Essential  for  cell  proliferation  20’ 68; 

oz: 

Regulation  of  actin  polymerization  ; 
Prognostic  marker  of  aggressive  breast 
cancer 50  and  prostate  cancer  20; 

Impair  DNA  repair  in  breast  cancer  cells 79, 

80. 

Overexpressed  in  B  cell  non-Hodgkin 
lymphoma  116,  bladder  cancer  51'54,  breast 
cancer  50,  gastric  cancer  55 ,  liver  cancer  57 , 

1  Q  1 

lung  cancer,  Hodgkin’s  lymphoma 
mantle  cell  lymphoma  132,  melanoma  58 
and  prostate  cancer  20’ 49’ 58 

Esc 

EED 

PRC2 

WD40  domain 

Interact  with  HDACs 

Substrate  preference  of  PRC2  '' 133 

Su(z)12 

SUZ12 

PRC2 

Zinc-finger 

domain 

Integrity  of  PRC2 
complex 

Fusion  protein  of  JAZF1-SUZ12  in  breast 
cancer  67 . 

Jing 

AEBP2 

PRC2 

three  Gli- 

Kriippel  (Cys2- 
His2)-type  zinc 
fingers 

Interact  with  and 

stimulate  the  activity  of 
the  PRC2  complex 

Substrate  preference  of  PRC2  '' 133 

Jarid2 

JARID2 

PRC2 

AT-rich 

interaction 

Interact  with  and 

Substrate  preference  of  PRC2  28-32 

domain 

(ARID) 

modulate  the  activity  of 
the  PRC2  complex 

Pci 

PHF1 

PRC2 

PHD  finger 

Interact  with  EZH2 

Enhance  EZH2  Histone  H3K27  tri- 

21 
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Table  2  PcG  proteins  and  associated  miRNAs 


Human 

proteins 

Complex 

miRNAs 

Associated  cancers 

EZH2 

PRC2 

miR-101  78’  91_ 

93 

miR-26a  93’ 95-97 
miR-214  98 

Non-Hodgkin  lymphoma,  breast,  bladder, 
liver,  lung  and  prostate  cancers 

EED 

PRC2 

miR-3()b  174 
miR-101  91 

Breast,  prostate  and  liver  cancers 

SUZ12 

PRC2 

miR-101  91 
miR-200b  135 

Metastatic  breast  cancer 

BMI1 

PRC1 

miR-15  136 
miR-16  136 
miR-128  137 

miR-200b,c 

100,  102 

Medulloblastoma,  breast  and  ovarian 

cancers 

miR-203  100’ 

102 

miR-429  100 

RING2 

PRC1 

miR-181a,b  1U2 
miR-200b,c  102 

Prostate  and  breast  cancers 

YY1 

PRC1 

PRC2 

miR-29 
miR-34a  139 

Cervical  cancer,  Rhabdomyosarcoma 
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Legends 

Figure  1  Polycomb  Repressive  Complexes.  (A)  Schematic  depicting  Polycomb 
Repressive  Complex  1  (PRC1).  RING1,  RING2  and  BMI1  form  the  core  of  Polycomb 
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Repressive  Complex  1  (PRC1).  PC3,  PH2,  RYBP,  YY1,  CtBP  and  E2F6  directly  or 
indirectly  bind  to  RING1,  RING2  or  BMI1.  (B)  Schematic  depicting  Polycomb 
Repressive  Complex  2  (PRC2).  EZH2,  EED,  SUZ12,  RbAp48,  and  AEBP2  consist  of  the 
core  of  Polycomb  Repressive  Complex  2  (PRC2).  Jarid2,  HDACs,  YY1,  VAV,  DNMT1, 
DNMT3A,  and  DNMT3B  interact  with  PRC2  and  regulate  its  activities. 

Figure  2  The  regulation  of  Polycomb  Repressive  Complexes.  (A)  PRC1  and  PRC2 
Regulatory  Axis.  PRC1  and  PRC2  are  known  to  be  regulated  by  E2F/Rb  pathway.  AR  has 
been  shown  to  regulate  ERG,  which  in  turn  regulates  PRC2.  miR-101,  miR-26a  and  miR- 
214  repress  PRC2  by  binding  to  their  3’UTRs.  PRC2  epigenetically  represses  miR-181a, 
b,  miR-200b,c,  and  miR-203  by  enacting  H3K27  trimethylation  at  these  particular  loci. 
These  miRNAs  in  turn  repress  PRC1  proteins,  BMI1,  RING1  and  RING2,  showing  that 
PRC2,  PRC1,  and  miRNAs  form  a  vital  regulatory  axis  in  the  epigenetic  regulation  of 
gene  expression.  (B)  PRC1  and  PRC2  recruitment  leads  to  gene  repression.  EZH2 
phosphorylation,  via  CDK1,  CDK2,  and  AKT,  has  been  shown  to  be  vital  in  terms  of 
regulating  EZH2  activity.  Acetylation  of  the  histone  tails  by  histone  acetyltransferases  on 
the  chromatin  marks  transcriptionally  active  DNA  regions.  Once  specific  molecular 
signals  are  released,  PRC2  binds  to  PRE  and  recruits  HDACs  that  deacetylate  the  histone 
tails,  followed  by  the  methylation  of  the  histone  tails  by  PRC2.  The  resulting  remodeled 
chromatin  mediates  access  to  the  promoters  of  the  target  genes  so  that  the  histone  tails  of 
the  promoters  can  in  turn  become  deacetylated  and  then  methylated.  Methylation  marks 
its  function  as  a  binding  site  for  the  N-terminal  chromodomain  of  PcG  proteins,  which 
recruits  PRC1  to  the  promoters  of  target  genes.  This  in  turn  results  in  the  inhibition  of 
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transcription  initiation.  The  ubiquitination  of  H2AK119  by  PRC1  leads  to  inactivation. 
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SUMMARY 

Polycomb  Repressive  Complexes  (PRC1  and  PRC2)-mediated  epigenetic  regulation  is  critical  for  maintain¬ 
ing  cellular  homeostasis.  Members  of  Polycomb  Group  (PcG)  proteins  including  EZH2,  a  PRC2  component, 
are  upregulated  in  various  cancer  types,  implicating  their  role  in  tumorigenesis.  Here,  we  have  identified 
several  microRNAs  (miRNAs)  that  are  repressed  by  EZH2.  These  miRNAs,  in  turn,  regulate  the  expression 
of  PRC1  proteins  BMI1  and  RING2.  We  found  that  ectopic  overexpression  of  EZH2-regulated  miRNAs  atten¬ 
uated  cancer  cell  growth  and  invasiveness,  and  abrogated  cancer  stem  cell  properties.  Importantly,  expres¬ 
sion  analysis  revealed  an  inverse  correlation  between  miRNA  and  PRC  protein  levels  in  cell  culture  and 
prostate  cancer  tissues.  Taken  together,  our  data  have  uncovered  a  coordinate  regulation  of  PRC1  and 
PRC2  activities  that  is  mediated  by  miRNAs. 


INTRODUCTION 

Polycomb  group  (PcG)  proteins  are  evolutionarily  conserved 
regulators  of  gene  silencing  important  in  metazoan  development 
(Surface  et  al.,  201 0),  stem  cell  pluripotency  (Pereira  et  al.,  201 0), 
and  X  chromosome  inactivation  (Cao  et  al.,  2002;  Margueron  and 
Reinberg,  2011).  PcG  proteins  form  multiprotein  repressive 
complexes  called  PRCs.  Both  PRC1  and  PRC2  play  a  critical 
role  in  the  maintenance  of  normal  and  cancer  stem  cell  pop¬ 
ulations  (Ezhkova  et  al.,  2009;  Lukacs  et  al.,  2010;  Pietersen 
et  al.,  2008).  Dysregulation  of  PcG  proteins  can  contribute  to 


a  number  of  human  diseases,  most  notably,  cancer  (Bracken 
and  Helin,  2009;  Margueron  and  Reinberg,  2011). 

Key  components  of  the  human  PRC2  include  the  histone 
methyltransferase  Enhancer  of  Zeste  Homolog  2  (EZH2),  and 
its  binding  partners,  Embryonic  Ectoderm  Development  (EED) 
and  Suppressor  of  Zeste  12  (SUZ12),  which  function  as  a  multi¬ 
subunit  complex  that  trimethylates  histone  H3K27.  PRC2  is 
thought  to  be  recruited  to  target  genomic  loci  by  long  noncoding 
RNAs  (ncRNAs)  such  as  HOTAIR  (Gupta  et  al.,  2010;  Kaneko 
et  al.,  2010;  Rinn  et  al.,  2007).  EZH2,  which  is  the  enzymatic 
component  of  PRC2,  is  elevated  in  aggressive  forms  of  prostate 


Significance 

Polycomb  group  (PcG)  proteins  are  chromatin-modifying  complexes  that  regulate  epigenetic  silencing  and  play  an  impor¬ 
tant  role  in  determining  cell  fate.  PcG  proteins  form  two  major  complexes,  Polycomb  Repressive  Complex  1  (PRC1)  and 
Polycomb  Repressive  Complex  2  (PRC2).  PRC2  methylates  histone  H3  on  Iysine27  (H3K27),  a  chromatin  mark  that  stimu¬ 
lates  PRC1  to  enact  gene  silencing  at  target  genes.  Employing  in  vitro  and  in  vivo  cancer  models  and  human  tumor  studies, 
we  demonstrate  that  PRC2  and  PRC1  coordinate  their  functions  through  regulation  of  specific  microRNAs.  Increased  PRC2 
activity  in  cancer  leads  to  repression  of  these  microRNAs,  and  subsequent  increase  of  PRC1  components.  Thus,  we  pro¬ 
pose  that  key  microRNAs  link  PRC2  to  PRC1  forming  an  integral  regulatory  axis  of  the  epigenetic  silencing  machinery. 
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and  breast  cancer  (Kleer  et  al.,  2003;  Varambally  et  al.,  2002),  as 
well  as  multiple  other  solid  tumors  (Matsukawa  et  al.,  2006;  Sudo 
et  al.,  2005).  Loss  of  microRNA  (miRNA)-1 01 ,  has  been  shown  to 
be  one  mechanism  that  leads  to  elevated  EZH2  and  PRC2 
activity  in  tumors  (Cao  et  al.,  201 0;  Chiang  et  al.,  201 0;  Friedman 
et  al.,  2009;  Varambally  et  al.,  2008;  Wang  et  al.,  2010).  Also, 
miR-26a  was  reported  to  target  EZH2  in  cancer  and  myogenesis 
(Lu  et  al.,  201 1 ;  Wong  and  Tellam,  2008).  Accumulating  evidence 
suggests  that  increased  activity  of  PRC2  is  oncogenic  as 
measured  by  cell  proliferation  (Bracken  et  al.,  2003;  Varambally 
et  al.,  2002),  cell  invasion  (Cao  et  al.,  2008;  Kleer  et  al.,  2003), 
anchorage-independent  growth  (Bracken  et  al.,  2003;  Kleer 
et  al.,  2003),  maintenance  of  tumor-initiating  cells,  tumor  xeno¬ 
graft  growth  (Yu  et  al.,  2007b),  and  metastasis  in  vivo  (Min 
et  al.,  2010)  . 

A  key  collaborator  of  PRC2  in  epigenetic  silencing  is  human 
PRC1,  which  comprises  B  lymphoma  Mo-MLV  insertion  region 
1  (BMI1),  RING1  (also  known  as  RING1A  or  RNF1)  and  RING2 
(also  known  as  RING1 B  or  RNF2),  and  functions  as  a  multiprotein 
complex  to  ubiquitinate  histone  H2A  at  lysine  119  (uH2A)  (Cao 
et  al.,  2005;  Wang  et  al.,  2004).  The  prevailing  hypothesis  is 
that  PRC2-mediated  trimethylation  of  FI3K27  recruits  PRC1  to 
gene  loci,  which  enacts  chromatin  condensation  and  epigenetic 
silencing  of  target  genes  (Bracken  and  Helin,  2009).  Like  PRC2 
component  EZH2,  BMI1  and  RING2  have  been  shown  to  be 
elevated  in  a  number  of  tumor  types  (Glinsky  et  al.,  2005;  San- 
chez-Beato  et  al.,  2006)  and  regulate  self-renewal  of  embryonic 
stem  cells  and  cancer  stem  cells  (Galmozzi  et  al.,  2006;  Valk- 
Lingbeek  et  al.,  2004).  The  mechanism  of  how  PRC2  and 
PRC1  coordinate  their  functions  is  still  unclear.  In  this  study, 
we  sought  to  explore  the  regulatory  axis  between  PRCs  and 
whether  miRNAs  mediate  the  synergy  between  the  two 
complexes. 

RESULTS 

PcG  Proteins  Are  Regulated  by  miRNAs 

Previously,  it  has  been  reported  that  EZH2,  the  methyltransfer- 
ase  subunit  of  the  PRC2  complex,  is  repressed  by  miR-101 
(Friedman  et  al.,  2009;  Varambally  et  al.,  2008)  and  miR-26a 
(Lu  et  al.,  2011;  Wong  and  Tellam,  2008).  We  hypothesized 
that  PcG  proteins  (comprising  the  mammalian  PRC  complexes) 
may  in  general  be  regulated  by  miRNAs.  To  test  this  hypothesis, 
we  knocked  down  Dicer,  a  key  protein  required  for  miRNA  pro¬ 
cessing,  by  employing  Dicer-specific  siRNA  duplexes.  By  immu- 
noblot  analysis,  we  found  that  PRC2  proteins  EZH2,  EED,  and 
SUZ12,  and  PRC1  proteins  BMI1  and  RING2  were  increased 
significantly  by  three  different  Dicer  siRNA  duplexes  (Figure  1  A; 
see  Figure  SI  A  available  online).  These  experiments  support  the 
general  notion  that  miRNAs  function  to  repress  PcG  expression. 

Identification  of  EZH2-Regulated  miRNAs 

To  explore  miRNAs  regulated  by  PRC2  globally,  we  knocked 
down  EZH2  in  DU145  prostate  cancer  cells  with  a  validated 
siRNA  targeting  EZH2  and  monitored  miRNA  expression  with 
lllumina  BeadChips.  In  parallel,  we  compared  these  miRNA 
profiles  with  DU145  cells  relative  to  four  benign  epithelial  cell 
lines  of  either  prostate  (PrEC  and  RWPE)  or  breast  (H16N2 
and  HME)  origin.  We  primarily  observed  miRNAs  that  were 


decreased  in  cancer  cells  relative  to  benign  that  are  targets  of 
repression  by  EZH2,  and  thus  PRC2.  We  found  63  miRNAs 
that  were  downregulated  in  DU145  cells  compared  with  the 
normal  cell  lines,  and  inhibition  of  EZH2  by  knockdown  restored 
expression  of  these  miRNAs  (Figure  1 B;  Table  SI).  Similarly,  the 
expression  levels  of  these  63  miRNAs  were  downregulated  in 
breast  cancer  cells  BT-549  and  SKBr3  compared  with  breast 
benign  epithelial  cells  H16N2  and  HME  (Figure  IB;  Table  SI) 
Using  miRNA  target  analysis  (www.targetscan.org),  we  identified 
14  miRNAs  as  top  candidates  with  the  following  properties:  (1) 
upregulated  by  EZH2  knockdown  in  DU145  cancer  cells  which 
express  high  levels  of  PRC2;  (2)  higher  in  benign  cell  lines 
compared  with  DU145  cells,  and  (3)  predicted  to  bind  to  the  3' 
untranslated  region  (UTR)  of  target  PRC1  components  based 
on  TargetScan  (Figure  1C).  Thirteen  of  the  14  miRNAs  meeting 
these  criterion  fell  into  several  known  miRNAs  clusters  and  fami¬ 
lies,  including  miR-200b  and  miR-200c  in  the  miR-200  family, 
which  has  previously  been  reported  to  repress  BMI1  (Shimono 
et  al.,  2009;  Wellner  et  al.,  2009).  Of  the  14  miRNAs,  only  miR- 
203,  which  is  also  known  to  target  BMI1  (Wellner  et  al.,  2009), 
does  not  belong  to  any  known  cluster  or  family  (Figure  SI  B). 

EZH2-Regulated  microRNAs  Inhibit  Expression  of  PRC1 
Proteins  BMI1  and  RING2 

To  pinpoint  the  specific  miRNAs  that  target  PRC1  (out  of  the  14 
that  were  nominated  by  computational  approaches)  (Figure  1C), 
we  overexpressed  each  of  them  in  BT-549  and  DU145  cancer 
cell  lines  and  monitored  EZH2,  BMI1 ,  and  RING2  protein  expres¬ 
sion  (Figure  2A;  Figure  S2A).  Of  these,  miR-181a,  b  decreased 
RING2  protein  levels,  miR-203  decreased  BMI1  protein  levels 
while  miR-200b,c  decreased  both  BMI1  and  RING2  (Figure  2A). 
Attenuation  of  these  PRC1  members  resulted  in  decreased  glo¬ 
bal  ubiquityl-H2A,  a  known  PRC1  substrate  and  mark  of  gene 
repression.  Furthermore,  PRC1  targets  including  p16INK4A 
(Jacobs  et  al.,  1999a)  and  p21  (Wafl/Cip)  (Fasano  et  al.,  2007) 
were  derepressed  (Figure  2A).  Several  of  the  miRNAs  com¬ 
putationally  predicted  to  inhibit  PRC1  failed  to  do  so  by  overex¬ 
pression  including  miR-17,  miR-19b,  and  others  (Figure  S2A). 
Similar  to  protein  levels,  real-time  qPCR  showed  miR-181a,b 
and  miR-200b,c  decreased  RING2  transcript  levels  and  miR- 
200b, c  and  miR-203  decreased  BMI1  transcript  levels  in 
BT-549  cells  (Figure  2B).  As  expected,  overexpressing  miR- 
200b  or  miR-203  decreased  BMI1  occupancy  on  known  PRC1 
target  gene  pi 6,  pi 9  (Jacobs  et  al.,  1999b),  p21,  and  HoxC13 
(Cao  et  al.,  2005)  regions  (Figure  S2B). 

To  further  corroborate  our  miRNA  overexpression  studies,  we 
also  extinguished  expression  of  miRNAs  using  antagomiRs 
(Krutzfeldt  et  al.,  2005).  Consistent  with  our  predictions,  anta- 
gomiR-200b,  antagomiR-200c,  and  antagomiR-203  increased 
BMI1  protein  levels,  while  antagomiR-181a,  antagomiR-181  b, 
antagomiR-200b,  and  antagomiR-200c  increased  RING2  pro¬ 
tein  levels  in  H16N2  cells  (Figure  2C). 

To  evaluate  whether  these  miRNAs  directly  bind  to  the  3'  UTR 
of  BMI1  or  RING2,  we  cloned  the  predicted  binding  sites  of  the 
wild-type  or  mutant  3'  UTR  into  a  luciferase  reporter  system 
and  cotransfected  them  with  miRNA  expression  vectors  into 
BT-549  cells  (Figure  2D;  Figures  S2C-S2F).  As  expected,  inhibi¬ 
tion  of  luciferase  activity  was  observed  in  cells  transfected  with 
constructs  containing  wild-type  binding  sites  but  not  the  mutant 
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Figure  1.  PcG  Proteins  Are  Regulated  by 
miRNAs 

(A)  Knockdown  of  Dicer  in  DU145  and  BT-549  cells 
by  three  different  Dicer-specific  duplexes  and  PcG 
protein  expression  was  assessed. 

(B)  miRNA  profiling  of  DU1 45  prostate  cancer  cells 
in  which  EZH2  was  knocked  down  compared  with 
DU145  cancer  cells  relative  to  benign  cells  HME, 
PrEC,  RWPE,  and  HI  6N2.  Shades  of  red  represent 
increased  gene  expression  while  shades  of  green 
represent  decreased  expression. 

(C)  A  Venn  diagram  depicting  1 4  miRNAs  that  were 
upregulated  by  EZH2  knockdown,  had  high 
endogenous  levels  in  normal  cells,  and  were  pre¬ 
dicted  to  target  PRC1  proteins. 

See  also  Figure  SI  and  Table  SI . 


constructs.  The  RING2  3'  UTR  reporters  were  down  regulated  by 
miR-181a,  miR-181b,  miR-200b,  and  miR-200c  while  the  BMI1 
3'  UTR  reporters  were  downregulated  by  miR-200b,  miR-200c, 
and  miR-203  (Figure  2D). 

We  next  determined  whether  the  miRNAs  that  regulate  PRC1 
were  directly  regulated  by  PRC2  in  BT-549  and  DU145  cells. 
Cells  were  transfected  with  either  a  validated  EZH2  siRNA  or 
miR-101  (both  of  which  target  and  downregulate  the  PRC2), 
and  expression  levels  of  target  miRNAs  were  measured  by 
real-time  PCR.  miR-181a,  miR-181b,  miR-200a,  miR-200b, 
miR-200c,  and  miR-203  expression  levels  were  increased  in 
EZH2  siRNA  or  miR-101 -transfected  cells.  Expression  of 
miRNAs  miR-217  and  miR-219,  two  control  microRNAs  not  pre¬ 
dicted  to  be  regulated  by  EZH2,  were  not  altered  (Figure  3A). 


Further,  we  observed  similar  expression 
changes  in  these  microRNAs  upon  sta¬ 
ble  overexpression  of  miR-101  or  EZH2 
shRNA  in  DU145  and  SKBr3  cells  (Fig¬ 
ure  S3A).  Also  we  observed  that  miR- 
101  was  increased  in  DU145  cells  in 
which  EZH2  was  stably  knocked  down, 
suggesting  the  existence  of  feedback 
regulation  between  EZH2  and  miR-101. 
In  contrast,  overexpression  of  EZH2,  but 
not  EZH2ASET  (which  is  missing  its  cata¬ 
lytic  SET  domain),  decreased  miR-181a, 
miR-181  b,  miR-200a,  miR-200b,  miR- 
200c,  and  miR-203  levels  in  H16N2  cells 
(Figure  S3B). 

Next,  we  treated  DU145  cells  with  the 
global  histone  methylation  inhibitor,  dea- 
zaneplanocin  A  (DZNep),  that  depletes 
PRC2  and  thus  attenuates  H3K27me3 
(Tan  et  al.,  2007).  Interestingly,  DZNep 
treatment  led  to  derepression  of  the  puta¬ 
tive  PRC2-targeted  miRNAs  including  miR 
-181a,b,  miR-200a,b,c,  and  miR-203 
(Figure  3B).  This  effect  was  both  concen¬ 
tration  and  incubation  time  dependent. 
Control  microRNAs,  miR-217,  miR-219, 
and  miR-21  were  not  affected  by  DZNep 
treatment. 

In  addition  to  DZNep,  we  evaluated  other  chemical  inhibitors  of 
epigenetic  pathways.  As  HDAC  activity  is  essential  for  EZH2 
function  (Cao  et  al.,  2008;  Kleer  et  al.,  2003),  and  EZH2  directly 
or  indirectly  facilitates  DNA  methylation  (Vire  et  al.,  2006),  we  pre¬ 
dicted  that  treatment  with  the  HDAC  inhibitor  suberoylanilide 
hydroxamic  acid  (SAHA)  and/or  the  DNA  methylation  inhibitor 
5-aza-2'-deoxycytidine  (5-aza-dC)  would  inhibit  EZH2-mediated 
epigenetic  modifications,  leading  to  an  increase  in  miRNA 
expression.  Treatment  of  BT-549  and  DU145  cells  with  5-aza- 
dC  or  SAHA  alone  or  in  combination,  resulted  in  a  marked  incre¬ 
ase  in  miR-1 81  a,b,  miR-200a,b,c,  and  miR-203  expression,  sug¬ 
gesting  epigenetic  regulation  of  these  microRNAs  (Figure  3C). 

Importantly,  when  we  overexpressed  EZH2  by  adenovirus  in 
DZNep  or  SAHA  and  5-aza-dC-treated  DU1 45  cells,  EZH2  could 
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completely  abolish  DZNep-mediated  miRNA  upregulation  (Fig¬ 
ure  S3C),  and  partially  decreased  SAHA  and  5-aza-dC-mediated 
miRNA  upregulation  (Figure  S3D)  presumably  because  SAHA 
and  5-aza-dC  also  inhibited  HDAC  and  DNMT  activities. 

To  confirm  that  EZH2  regulates  these  microRNAs  by  epige¬ 
netic  repression,  we  performed  chromatin  immunoprecipitation 
(ChIP)  assays  with  anti-H3K27me3,  EZH2  and  BMI1  antibodies 
in  BT-549  cells.  Interestingly,  H3K27me3  and  EZH2  occupied 
the  PRC2-regulated  miRNAs  regions  as  expected.  In  addition, 
BMI1  also  occupied  these  regions  (Figure  S3E),  suggesting 


Figure  2.  PRC2-Regulated  miRNAs  Repress 
PRC1  Proteins  BMI1  and  RING2 

(A)  Overexpression  of  indicated  miRs  in  DU145 
and  BT-549  cells  and  expression  of  PRC  compo¬ 
nents,  PRC2  histone  mark  H3K27me3,  PRC1 
target  histone  mark  ubiquityl-H2A  and  indicated 
genes  by  immunoblot  analysis.  (3-actin  and  total 
H3  were  used  as  loading  controls. 

(B)  As  in  (A),  except  transcript  level  was  assessed 
in  BT-549  by  qPCR. 

(C)  Transfection  of  indicated  antagomiRs  (anti- 
miR)  in  H16N2  cells  and  immunoblot  analysis  for 
BMI1  and  RING2.  |3-actin  was  used  as  a  loading 
control. 

(D)  TargetScan  analysis  depicting  potential 
binding  sites  for  EZH2-regulated  miRNAs  in  the  3' 
UTR  of  BMI1  and  RING2.  Luciferase  reporter 
assays  with  wild-type  or  mutant  3'  UTR  constructs 
of  BMI1  or  RING2  demonstrate  that  miR-1 81  a, 
miR-1 81b,  miR-200b,  miR-200c,  and  miR-203 
repress  BMI1  and/or  RING2  activity. 

All  bar  graphs  are  shown  with  ±  SEM.  See  also 
Figure  S2. 


that  a  negative  feedback  system  between 
PRC2-regulated  miRNAs  and  PRC1  may 
exist.  Furthermore,  an  EZH2-specific 
siRNA  (Figure  S3F)  or  treatment  with 
5-aza-dC  and  SAHA,  either  alone  or  in 
combination  (Figure  3D),  markedly  decre¬ 
ased  the  H3K27me3  occupancy  in  these 
regions. 


EZH2-Regulated  miRNAs  Attenuate 
Growth,  Invasiveness,  and  Self- 
Renewal  of  Cancer  Cells 

Because  EZH2  has  been  shown  to 
repress  several  tumor  suppressor  genes 
(Cao  et  al.,  2008;  Chen  et  al.,  2005;  Fujii 
et  al.,  2008;  Min  et  al.,  2010;  Yu  et  al., 
2007b,  2010),  we  postulated  that  the 
EZH2-regulated  microRNAs  also  func¬ 
tioned  as  tumor  suppressors.  Consistent 
with  this  hypothesis,  overexpression  of 
either  miR-1 81  a,  miR-1 81b,  miR-200a, 
miR-200b,  miR-200c,  or  miR-203  mark¬ 
edly  attenuated  BT-549  and  DU145  cell 
proliferation  to  levels  similar  to  that  of 
cells  transfected  with  EZH2  siRNA,  or 
cells  overexpressing  miR-101  (Figure  4A 
*p  <  0.001,  **p  <  0.01;  Figure  S4A).  Likewise,  overexpression 
of  either  miR-1 81  a,  miR-1 81b,  miR-200a,  miR-200b,  miR- 
200c,  or  miR-203  inhibited  the  in  vitro  invasive  potential  of 
BT-549  and  DU145  cells  through  modified  Boyden  chambers 
coated  with  Matrigel  (Figure  4B,  *p  <  0.005,  **p  <  0.02).  However, 
overexpressing  EZH2-repressed  miRNAs  had  no  effect  on  the 
invasiveness  of  RWPE-UBE2L3-KRAS  and  RWPE-SLC45A3- 
BRAF  stable  cells,  in  which  fusion  proteins  UBE2L3-KRAS 
(Wang  et  al.,  2011)  and  SLC45A3-BRAF  (Bond  et  al.,  2008;  Pal- 
anisamy  et  al.,  2010)  confer  neoplastic  properties  to  RWPE  cells 
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Figure  3.  PRC2  Silences  Multiple  miRNAs  by  Epigenetic  Mechanisms 

(A)  Taqman  miRNA  qPCR  analysis  of  indicated  miRs  in  BT-549  and  DU1 45  cells  in  which  EZH2  was  knocked  down  using  siRNA  or  miR-1 01  (a  microRNA  which 
targets  EZH2).  Quantitative  microRNA  levels  were  normalized  against  U6. 

(B)  As  in  (A),  except  DZNep  at  two  different  doses  and  time  points  was  incubated  with  DU145  cells. 

(C)  As  in  (A),  except  SAHA  and/or  5-aza-dC  was  used  in  BT-549  and  DU145  cells. 

(D)  ChIP-qPCR  analysis  of  H3K27me3  at  indicated  genes  and  microRNAs  in  BT-549  cells  treated  with  SAHA  and/or  5-aza-dC. 

All  bar  graphs  are  shown  with  ±  SEM.  See  also  Figure  S3. 


(Figure  S4B),  suggesting  that  EZH2-repressed  miRNAs  miR- 
181a,b,  miR-200b,c,  and  miR-203  may  inhibit  cell  invasion 
through  acting  on  PRC1  proteins.  However,  EZH2-repressed 
miRNAs  still  decreased  RWPE-UBE2L3-KRAS  and  RWPE- 
SLC45A3-BRAF  proliferation  (Figure  S4C),  consistent  with  a  crit¬ 
ical  role  of  PcG  proteins  in  cell  growth. 

To  investigate  whether  miR-1 81  a,  miR-1 81b,  miR-200a,  miR- 
200b,  miR-200c,  or  miR-203  inhibit  anchorage-independent 


growth,  we  performed  soft  agar  colony  formation  assays.  Similar 
to  miR-101  and  EZH2  knockdown  controls,  overexpression  of 
miR-1 81  a,  miR-1 81b,  miR-200a,  miR-200b,  miR-200c,  and 
miR-203  markedly  suppressed  DU145  colony  formation  (Fig¬ 
ure  4C,  *p  <  0.001 ,  **p  <  0.01).  Next,  we  evaluated  the  ability  of 
DU1 45  to  form  prostatospheres  in  sphere-promoting  cell  media. 
This  assay  serves  as  a  surrogate  measure  of  stem  cell-like 
phenotypes,  and  cells  that  are  able  to  form  spheres  have 
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Figure  4.  PRC2-Mediated  Regulation  of  microRNAs  Potentiates  the  Cancer  Cell  Phenotype 

(A)  Overexpression  of  PRC2- regulated  miRNAs,  but  not  control  miR-217  or  miR-219,  inhibited  BT-549cell  proliferation.  EZH2siRNAand  miR-101  overexpression 
were  positive  controls  and  miR-217  and  miR-219  overexpression  were  negative  controls.  *p  <  0.001 ,  **p  <  0.01 .  (Student’s  t  test). 

(B)  Overexpression  of  PRC2-regulated  miRNAs  decreased  BT-549  and  DU145  cell  invasion  in  vitro.  *p  <  0.01.  (Student’s  ttest). 

(C)  Overexpression  of  PRC2-regulated  miRNAs  suppressed  DU145  anchorage-independent  growth  in  soft  agar.  *p  <  0.01.  (Student’s  t  test). 

(D)  Overexpression  of  EZH2-regulated  miRNAs  decreased  prostatosphere  formation  by  DU145  cells.  *p  <  0.01.  (Student’s  t  test).  Representative  images  of 
prostatospheres  (scale  bar:  100  |im)  were  shown  in  the  inset. 

(E)  qPCR  analysis  demonstrating  EZH2,  BMI1  and  RING2  transcript  levels  were  higher  in  spheres  compared  with  monolayer  culture,  while  miR-101 ,  miR-181a,b, 
miR-200a,b,c,  and  miR-203,  but  not  miR-217  or  miR-219,  were  lower  in  spheres  compared  with  monolayers.  Expression  level  of  each  gene  was  normalized  to 
GAPDH  or  U6  and  normalized  to  corresponding  monolayer  cultured  cell  line. 

(F)  qPCR  analysis  showing  EZH2,  BMI1  and  RING2  levels  were  higher  in  sorted  CD24-/CD44+  DU145  and  RWPE  cells  compared  with  the  unsorted  population, 
while  miR-101,  miR-181a,b,  miR-200a,b,c,  and  miR-203,  but  not  miR-217  or  miR-219,  were  lower  in  CD24-/CD44+  DU145  and  RWPE  cells  compared  with  an 
unsorted  population. 

(G)  Genes  regulated  by  EZH2-repressed  miRNAs  cluster  into  multiple  functional  concepts.  BT-549  and  DU145  cells  were  transfected  with  EZH2-repressed 
miRNAs  followed  by  gene  expression  profiling  and  Molecular  Concepts  analysis.  Each  node  represents  a  molecular  concept  or  set  of  biologically  related 
genes.  miR-101,  miR-181a,  miR-181b,  miR-200a,  miR-200b,  miR-200c,  and  miR-203  (miRNA  signatures,  purple  for  BT-549,  orange  for  DU145)  were  enriched 
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enhanced  stem  cell  characteristics  (Lawson  et  al.,  2007).  We 
found  that  miR-181a,  miR-181b,  miR-200a,  miR-200b,  miR- 
200c,  and  miR-203  overexpression,  as  well  as  miR-101  overex¬ 
pression  and  EZH2  siRNA  controls,  significantly  inhibited  the 
ability  of  DU145  cells  to  form  spheres  in  this  assay  (Figure  4D, 
*p  <  0.001,  **p  <  0.01).  Intriguingly,  several  genes  implicated  in 
pluripotency  and  cellular  reprogramming  by  induced  pluripo- 
tency,  such  as  Klf4,  Sox2,  and  c-Myc,  were  markedly  downregu- 
lated  by  miR-200b,  miR-200c,  and  miR-203,  and  marginally 
decreased  by  miR-101,  miR-181a,  miR-181b,  and  miR-200a 
expression,  but  not  by  miR-21 7  or  miR-21 9  controls  (Figure  S4D). 
Relative  to  the  human  embryonic  stem  cell  H7,  BT-549  and 
DU145  cancer  cells  have  comparable  expression  levels  of  iPS 
factors  and  PcG  proteins  (Figure  S4E). 

Next,  we  measured  expression  levels  of  EZH2,  BMI1,  RING2, 
and  key  microRNAs  relevant  to  this  study  in  spheres  and  mono- 
layers.  In  BT-549,  SKBr3,  DU145,  and  PC3  cells,  we  observed 
that  EZH2,  BMI1,  and  RING2  levels  were  higher  in  spheres 
than  in  monolayers;  conversely  miR-101,  miR-181a,  miR-181b, 
miR-200a,  miR-200b,  miR-200c,  and  miR-203  levels  were  lower 
in  spheres  than  in  monolayers  (Figure  4E).  Using  DU145  and 
RWPE  parental  cell  lines,  we  employed  flow  cytometry  to  isolate 
cells  with  high  expression  of  the  CD44  surface  antigen  and  low 
expression  of  the  CD24  surface  antigen  (CD24-/CD44+),  a  cell 
population  enriched  for  stem  cell-like  phenotypes  (Hurt  et  al., 
2008).  We  measured  EZH2,  BMI1,  RING2,  and  miRNA  levels  in 
CD24-/CD44+  cells  compared  with  total,  unsorted  cells.  We 
observed  that  EZH2,  BMI1,  and  RING2  levels  were  increased 
in  CD24-/CD44+  cells,  but  miR-101,  miR-181a,  miR-181b, 
miR-200a,  miR-200b,  miR-200c,  and  miR-203  expression  were 
decreased  in  this  cell  population  (Figure  4F).  Taken  together, 
the  data  provide  compelling  evidence  for  the  coordinated  regu¬ 
lation  of  PRC2,  PRC1 ,  and  miRNAs  in  the  maintenance  of  a  differ¬ 
entiated  cellular  state  and  inhibition  of  stem  cell-like  phenotypes. 

In  order  to  understand  the  functional  biology  of  the  miRNAs 
identified  in  this  study,  we  sought  to  identify  global  gene  expres¬ 
sion  patterns  and  molecular  pathways  to  which  they  might 
contribute.  We  conducted  gene  expression  microarray  analyses 
of  DU145  and  BT-549  cells  transfected  with  control  miRNA,  miR- 
101,  miR-181a,  miR-181b,  miR-200a,  miR-200b,  miR-200c,  or 
miR-203.  As  shown  in  Table  S2  and  Table  S3,  EZH2-repressed 
miRNAs  targeted  many  predicted  genes.  When  we  analyzed 
the  miRNA-regulated  genes  using  Molecular  Concepts  Maps 
(MCM)  (Tomlins  et  al.,  2007b),  as  expected,  molecular  concepts 
associated  with  these  miRNAs  were  highly  overlapping,  showing 
a  high  correlation  to  gene  sets  representing  multiple  cancers, 
metastatic  cancer  processes,  cancer  survival,  Polycomb  Group 
targets,  and  stem  cell-related  genes  (Figure  4G;  Table  S4). 

In  order  to  further  examine  the  molecular  link  between  PRC1 
and  PRC2  activities,  we  generated  DU145  cells  stably  overex¬ 
pressing  miR-200b  and  miR-203  (Figure  S5A)  and  monitored 
levels  of  BMI1  and  RING2.  BMI1  and  RING2  were  decreased  in 
miR-200b  stable  cells  while  only  BMI1  was  decreased  in  miR- 
203  stable  cells.  In  addition,  uH2A,  the  histone  modification 
mediated  by  PRC1 ,  was  similarly  decreased  in  both  miR-200b- 


and  miR-203-expressing  cells.  Interestingly,  BMI1,  RING2,  and 
uH2A,  as  well  as  EZH2  and  H3K27me3,  were  decreased  in 
miR-101  stable  expressing  DU145  cells  (Figure  5A)  suggesting 
that  prolonged  knockdown  of  PRC2  components  leads  to 
suppression  of  PRC1.  Using  cell  count  and  Boyden  chamber 
invasion  assays,  we  found  that  similar  to  miR-101,  miR-200b 
and  miR-203  stably  expressing  cells  grew  more  slowly  and 
were  less  invasive  than  vector-transfected  cells  (Figures  5B 
and  5C).  Intriguingly,  coexpression  of  BMI1  or  EZH2  (control) 
without  the  3'  UTR  both  restored  the  proliferation  and  invasion 
properties  of  DU145  cells  despite  the  presence  of  miR-101, 
miR-200b,  or  miR-203  (Figures  5B  and  5C).  Importantly,  murine 
xenograft  experiments  demonstrated  that  DU145  cells  with 
stable  knockdown  of  PRC1  proteins  BMI1  or  RING2  (Figure  S5B), 
or  expressing  miR-181b  (Figure  S5C),  miR-200b,  or  miR-203 
grew  more  slowly  than  the  vector  control  in  vivo  (p  =  0.0001 , 
Figures  5D  and  5E). 

EZH2-Regulated  miRNAs  Inversely  Correlate  with  PRC 
Protein  Levels  in  Prostate  Cancer 

Since  miR-101,  miR-181a,  miR-181b,  miR-200a,  miR-200b, 
miR-200c,  and  miR-203  appear  to  play  an  important  role  in 
cancer  progression,  we  next  measured  the  endogenous  expres¬ 
sion  levels  of  these  miRNAs  by  qPCR  analysis  of  a  cohort  of 
benign  prostate,  localized,  and  metastatic  prostate  cancers  in 
which  we  had  measured  miR-101,  miR-21 7,  and  EZH2  levels 
previously  (Varambally  et  al.,  2008).  As  expected,  miR-181a, 
miR-181b,  miR-200a,  miR-200b,  miR-200c,  and  miR-203  levels 
were  lowest  in  metastatic  prostate  cancer  tissues,  and  highest  in 
benign  prostate  tissues  (Figure  6A).  In  addition,  immunoblot 
analyses  showed  that  BMI1,  RING2,  and  uH2A,  as  well  as 
EZH2,  but  not  RING1,  were  increased  in  metastatic  prostate 
cancer  compared  with  benign  tissues  and  localized  cancer 
samples  (Figure  6B;  Figure  S6A).  EZH2  levels  were  highly  corre¬ 
lated  with  BMI1 ,  RING2,  and  H2A  protein  levels  (Figure  S6B), 
further  supporting  a  molecular  link  between  PRC1  and  PRC2 
expression  and  activities  during  cancer  progression.  As  ex¬ 
pected,  ChIP  assays  showed  that  H3K27me3-marked  chromatin 
occupied  the  miR-203  upstream  region  in  metastatic  prostate 
cancer,  but  not  in  localized  prostate  cancer  (PCA)  (Figure  S6C). 
Similarly,  DNA  methylation  of  the  miR-203  genomic  region  was 
observed  in  localized  and  metastatic  prostate  cancer  but  not 
benign  prostate  tissue  (Figure  6C).  Taken  together,  these  data 
suggest  that  EZH2-mediated  epigenetic  repression  of  miR- 
181a,  miR-181b,  miR-200b,  miR-200c,  and  miR-203  results  in 
an  upregulation  of  PRC1  proteins  BMI1  and  RING2  and  histone 
code  ubiquityl-H2A  in  advanced  prostate  cancer. 

DISCUSSION 

This  study  unravels  the  intricacies  in  the  regulation  of  the  poly¬ 
comb  protein  complexes  mediated  by  various  miRNAs,  and 
substantiates  the  essential  role  played  by  PRC  in  cancer.  We 
demonstrated  that  increased  PRC2  activity  results  in  repression 
of  numerous  miRNAs  that  are  known  to  be  important  in  the 


for  concepts  related  to  cancer  (yellow),  cancer  survival  (red),  stem  cell  likeness  (blue),  and  function  of  polycomb  group  (green).  All  bar  graphs  are  shown 
with  ±SEM. 

See  also  Rgure  S4,  and  Table  S2,  Table  S3,  and  Table  S4. 
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Figure  5.  PRC2-Repressed  miRNAs  Inhibit 
Tumor  Growth 

(A)  DU145  cells  stably  overexpressing  miR-101, 
miR-200b,  and  miR-203  demonstrated  repression 
of  EZH2,  BMI1 ,  or  RING2,  as  well  as  decreased 
H3K27me3  and  ubiquityl-H2A  (uH2A)  levels. 

(B  and  C)  Coexpression  of  EZH2A3/UTR  or 
BMIIAS'UTR  rescued  cell  proliferation  (B),  and 
invasiveness  (C)  of  DU145  cells  stably  overex¬ 
pressing  miR-101,  miR-203,  or  miR-200b. 

(D)  Stably  knocking  down  BMI1  or  RING2  by 
BMI1 -specific  shRNA  (BMI1-sh3)  or  RING2-spe- 
cific  shRNA  (RING2-sh1)  decreased  Dili 45  tumor 
growth  in  mice.  N  =  8  for  DU145  control  (scra¬ 
mble),  BMI1-sh3,  and  RING2-sh1,  respectively, 
were  used  for  the  xenograft. 

(E)  Stable  overexpression  of  miR-1 81  b,  miR-200b, 
or  miR-203  decreased  DU145  tumor  growth  in 
mice.  DU145  miR-vector  (N  =  9),  miR-NT  (non¬ 
targeting)  (N  =  8),  miR-1 81b  (N  =  8),  miR-200b 
(N  =  8),  or  miR-203  (N  =  7)  were  used  for  the 
xenograft  experiment.  DU145  stable  pools  were 
injected  subcutaneously. 

All  bar  graphs  are  shown  with  ±  SEM.  See  also 
Figure  S5. 
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maintenance  of  stem  cell-like  phenotypes  in  cancer  cells.  We 
show  that  PRC2  epigenetically  represses  miR-1 81  a,  miR-1 81b, 
miR-200b,  miR-200c,  and  miR-203  expression  by  facilitating 
H3K27me3  trimethylation  at  these  loci,  and  that  exogenous 
overexpression  of  miR-1 81  a,  miR-1 81b,  miR-200b,  miR-200c, 
and  miR-203  inhibits  a  cancer  phenotype  in  vitro.  Furthermore, 
miR-1 81b,  miR-200b,  and  miR-203  overexpression  suppressed 
prostate  tumor  formation  and  growth  in  mouse  xenografts. 
Recently,  several  groups  have  also  reported  roles  for  miR- 
200b,  miR-200c,  and  miR-203  in  controlling  stem  cell  differenti¬ 
ation  (Yi  et  al.,  2008),  epithelial-to-mesenchymal  transition  (EMT) 
(Park  et  al.,  2008;  Wellner  et  al.,  2009),  and  cancer  progression 
(Faber  et  al.,  2008;  Shimono  et  al.,  2009). 

Here,  we  demonstrated  that  PRC1  proteins  BMI1  and  RING2 
are  direct  targets  of  miR-1 81  a,  miR-1 81b,  miR-200b,  miR- 
200c,  and  miR-203  in  breast  and  prostate  cancer.  Furthermore, 
we  observed  a  significant  negative  correlation  between  PRC2 
expression  and  miR-1 81  a,  miR-1 81b,  miR-200b,  miR-200c, 
and  miR-203,  as  well  as  a  strong  positive  correlation  between 
EZH2,  BMI1 ,  and  RING2  protein  levels.  Intriguingly,  earlier 
studies  suggested  a  discrepancy  between  BMI1  protein  and 


RNA  levels  in  prostate  tissues,  as  protein 
levels  were  increased  while  RNA  levels 
were  decreased  during  prostate  cancer 
progression  (Varambally  et  al.,  2005).  It 
is  possible  that  regulation  of  PRC 
proteins  occurs  at  both  transcriptional 
and  posttranscriptional  levels  by  sepa¬ 
rate  mechanisms.  We  provide  evidence 
that  EZH2-regulated  microRNAs  con¬ 
tribute  to  the  maintenance  of  a  differenti¬ 
ated  cellular  state,  and  that  miR-1 81  a, 
miR-1 81b,  miR-200b,  miR-200c,  and 
miR-203  function  as  tumor  suppressors  during  prostate  cancer 
progression. 

Interestingly,  several  recent  studies  have  reported  similar 
microRNA-protein  regulatory  networks  that  play  critical  roles  in 
cancer.  In  one  study,  the  RAS  proto-oncogene  was  shown  to 
be  coordinately  regulated  by  the  let-7  family  of  miRs  (Johnson 
et  al.,  2005).  Likewise,  the  miR-1 5a-miR-1 6-1  cluster,  located 
on  chr13q14,  was  proposed  to  serve  as  a  tumor  suppressor  in 
prostate  tissue  by  regulating  levels  of  cancer-related  genes 
such  as  BCL2,  CCND1,  and  WNT3A  (Bond  et  al.,  2008). 
Recently,  Poliseno  et  al.  (2010)  reported  a  proto-oncogenic 
miRNA-dependent  network  in  prostate  cancer  progression  in 
which  the  miR-1 06b~25  cluster  regulates  PTEN  expression 
and  cooperates  with  MCM7  in  cellular  transformation.  These 
studies,  along  with  our  present  study,  strongly  suggest  that  dys- 
regulation  of  miRNA  and  target  protein  networks  may  contribute 
to  cancer  development. 

Here,  we  propose  a  model  for  a  coordinated  PRC2-PRC1 
oncoprotein  axis,  and  epigenetic  link  between  H3K27me3 
and  ubiquityl-H2A,  mediated  by  PRC2-regulated  miRNAs  (Fig¬ 
ure  7).  Recently,  lliopoulos  et  al.  (2010)  reported  that  miR-200b 
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Figure  6.  Coordinated  Expression  of  PcG  Proteins  and  PRC  Regulatory  miRNAs  in  Prostate  Cancer  Progression 

(A)  Expression  of  indicated  miRs  as  assessed  by  q-PCR  in  benign  prostate,  clinically  localized  prostate  cancer  and  metastatic  prostate  cancer  tissues.  Data  for 
EZH2,  miR-217,  and  miR-101  were  reported  previously  (Varambally  et  al.,  2008)  and  displayed  here  for  comparison  (Student’s  t  test). 

(B)  Immunoblot  analysis  of  EZH2,  BMI-1 ,  RING2,  RING1 ,  and  ubiquityl-H2Ain  benign  prostate,  clinically  localized  prostate  cancer,  and  metastatic  prostate  cancer. 

(C)  Bisulfite  sequencing  analysis  of  the  miR-203  genomic  region  revealed  cancer-specific  DNA  methylation  in  a  region  proximal  to  miR-203  in  prostate  cancer 
tissues. 

All  bar  graphs  are  shown  with  ±  SEM.  See  also  Figure  S6. 


regulates  PRC2  protein  SUZ1 2  in  a  manner  similar  to  that  of  miR- 
101,  lending  further  support  for  microRNA-mediated  PRC 
activity  during  cancer  progression.  These  findings  offer  multiple 
targets  for  therapeutic  interventions  in  the  treatment  of  aggres¬ 
sive  cancers  (Garzon  et  al.,  2010). 

EXPERIMENTAL  PROCEDURES 

Cell  Lines 

Breast  cancer  cell  line  BT-549  was  grown  in  RPM1 1640  (Invitrogen,  Carlsbad, 
CA)  with  0.023  lU/ml  insulin  and  10%  FBS  (Invitrogen)  in  5%  C02  cell  culture 


incubator;  breast  cancer  cell  line  SKBr3  was  grown  in  RPMI  1640  (Invitrogen) 
with  10%  FBS  (Invitrogen)  in  5%  C02  cell  culture  incubator;  and  prostate 
cancer  cell  line  DU145  was  grown  in  MEM  with  10%  FBS  in  5%  C02  cell 
culture  incubator.  Immortalized  breast  cell  lines  HME  and  H16N2  were  grown 
in  F-12  Nutrient  Mixture  with  5^ig/ml  Insulin  (Sigma,  St.  Louis,  MO),  1  ^ig/ml 
Flydrocortisone  (Sigma),  10  ng/ml  EGF  (Invitrogen),  5  mM  Ethanolamine 
(Sigma),  5  ^ig/ml  Transferrin  (Sigma),  10  nM  Triiodo  Thyronine  (Sigma), 
50  nM  Sodium  Selenite  (Sigma),  10  mM  FIEPES  (Invitrogen)  and  50  unit/ml 
Penstrep  (Invitrogen),  10%  C02.  The  PrEC  (Lonza,  Conshohocken,  PA)  and 
RWPE  (ATCC,  Manassas,  VA)  cells  were  grown  in  their  respective  medium 
as  specified  by  the  suppliers.  miR-181b,  miR-200b,  and  miR-203  overex¬ 
pression  constructs  were  obtained  from  Openbiosystems  (Huntsville,  AL). 


Cancer  Cell  20,  1 87-1 99,  August  1 6,  201 1  ©201 1  Elsevier  Inc.  1 95 


Cancer  Cell 

microRNAs  Link  PRC2  to  PRC1 


Cel 


PRC2 


miR-101  - 

miR-26a  I 


PRC1 


miR-181a 
miR-181b 
miR-200b  — | 
miR-200c 
miR-203 


H3K27 

Tri-methylation 


H2A 

ubiquitination 


Normal 

Cells 

Cancer/ 

pluripotent 

cells 


T 

i 


i 

T 


T 

i 


i 

T 


Lentiviruses  were  generated  by  the  University  of  Michigan  Vector  Core.  BMI1 , 
RING2  and  control  shRNA  lentivirus  were  obtain  from  Sigma.  Prostate  cancer 
cell  line  DU1 45  was  infected  with  lentiviruses  expressing  BMI1  shRNA,  RING2 
shRNA,  miR-181b,  miR-200b,  and  miR-203  or  controls  only,  and  stable  cell 
lines  were  generated  by  selection  with  300  [xg/ml  puromycin  (Invitrogen). 

Benign  and  Tumor  Tissues 

In  this  study,  we  utilized  tissues  from  clinically  localized  prostate  cancer 
patients  who  underwent  radical  prostatectomy  as  a  primary  therapy  between 
2004  and  2006  at  the  University  of  Michigan  Hospital.  Samples  were  also  used 
from  androgen-independent  metastatic  prostate  cancer  patients  from  a  rapid 
autopsy  program  described  previously  (Tomlins  et  al.,  2005,  2007a).  The 
detailed  clinical  and  pathological  data  are  maintained  in  a  secure  relational 
database.  This  study  was  approved  by  the  Institutional  Review  Board  at  the 
University  of  Michigan  Medical  School.  Informed  consent  was  also  obtained 
from  all  subjects  through  the  Institutional  Review  Board  at  the  University  of 
Michigan  Medical  School.  Both  radical  prostatectomy  series  and  the  rapid 
autopsy  program  are  part  of  the  University  of  Michigan  Prostate  Cancer 
Specialized  Program  of  Research  Excellence  Tissue  Core. 

Illumina  microRNA  Profiling 

Total  RNA  (500  ng)  from  each  sample  was  labeled  and  hybridized  on  the 
Human  v2  microRNA  Expression  BeadChips  (Illumina,  San  Diego,  CA)  accord¬ 
ing  to  the  manufacturers  recommendations.  BeadChips  were  scanned  with 
the  Illumina  iScan  Reader.  Data  were  then  average  median  normalized  before 
generating  differential  expression  values  between  treated  and  control 
samples. 

microRNA  Transfection,  AntagomiR  Transfection,  and  Small  RNA 
Interference 

Knockdown  of  EZH2  or  Dicer  was  accomplished  by  RNA  interference  using 
siRNA  duplexes  (Dharmacon,  Lafayette,  CO)  as  previously  described  (Varam- 
bally  et  al.,  2002).  Precursors  of  respective  microRNAs,  antagomiRs  and  nega¬ 
tive  controls  were  purchased  from  Ambion  (Austin,  TX).  Transfections  were 
performed  with  oligofectamine  (Invitrogen).  EZH2  siRNA  duplexes  sequences, 
(duplex  1:  GAG  GTT  C  AG  AC  GAG  CT  GAT ;  duplex  2:  AGACTCT  GAATGCA 
GTTGC). 


Figure  7.  A  Proposed  Model  Role  for  micro¬ 
RNAs  in  Regulating  PRCs 

Specifically,  PRC2  is  molecularly  linked  to  PRC1 
via  a  set  of  regulatory  miRs. 


every  other  day)  and/or  1[xM  suberoylanilide 
hydroxamic  acid  (SAHA)  for  2  days.  DU145  cells 
were  treated  with  2.5  or  5  ^M  deazaneplanocin  A 
(DZNep)  for  2  or  3  days  followed  by  RNA  extraction 
or  chromatin  immunoprecipitation. 

Cell  Proliferation  Assay  and  Basement 
Membrane  Matrix  Invasion  Assays 

Invasive  breast  cancer  cell  BT-549  and  prostate 
cancer  cell  DU145  were  transfected  with  miRNAs 
or  controls.  The  cell  proliferation  and  invasion 
assays  were  performed  as  described  (Cao  et  al.,  2008;  Kleer  et  al.,  2003;  Var- 
ambally  et  al.,  2008;  Yu  et  al.,  2007b). 

Soft  Agar  Colony  Formation  Assays 

A  50  |xl  base  layer  of  agar  (0.6%  Agar  in  DM  EM  with  1 0%  FBS)  was  allowed  to 
solidify  in  a  96-well  flat-bottom  plate  prior  to  the  addition  of  a  75  ^1  miRNAs  or 
control-transfected  or  stable  DU145  cell  suspension  containing  4000  cells  in 
0.4%  Agar  in  DMEM  with  10%  FBS.  The  cell  containing  layer  was  then  solid¬ 
ified  at  4C  for  1 5  min  prior  to  the  addition  of  1 00  pil  of  MEM  with  5%  FBS.  Colo¬ 
nies  were  allowed  to  grow  for  21  days  followed  by  counting  and  imaging  under 
a  light  microscope. 

Spheres  Culture 

Spheres  culture  was  performed  as  described  (Dontu  et  al.,  2003;  Yu  et  al., 
2007a).  Briefly,  cells  (1000  cells/ml)  were  cultured  in  suspension  in  serum- 
free  DMEM-F12  (Invitrogen),  supplemented  with  B27  (1:50,  Invitrogen), 
20  ng/ml  EGF  (BD  Biosciences),  0.4%  bovine  serum  albumin  (Sigma),  and 
4  [xg/ml  insulin  (Sigma).  To  propagate  spheres  in  vitro,  spheres  were  collected 
by  gentle  centrifugation,  dissociated  to  single  cells  as  described  (Dontu  et  al., 
2003;  Yu  et  al.,  2007a),  and  then  cultured  to  generate  prostatospheres  of  the 
next  generation.  Spheres  larger  than  50  [xm  were  counted. 

Gene  Expression  Profiling 

Expression  profiling  was  performed  using  the  Agilent  Whole  Human  Genome 
Oligo  Microarray  (Santa  Clara,  CA)  according  to  the  manufacturer’s  protocol. 
BT-549  and  DU145  cells  were  transfected  with  miRNAs  or  negative  control  for 
precursor  microRNA.  Over-  and  underexpressed  signatures  were  generated 
by  filtering  to  include  only  features  with  significant  differential  expression 
(Log  ratio,  p  <  .01)  in  all  hybridizations  and  2-fold  average  over-  or  under¬ 
expression  (Log  ratio).  Gene  expression  data  are  deposited  into  GEO 
(GSE26996). 

Gene  Set  Enrichment  Analysis 

Molecular  Concept  Map  (MCM)  analysis  was  performed  using  gene  list  of 
putative  targets  to  search  for  all  concepts  available  in  the  Oncomine  database 
as  previously  described  (Yu  et  al.,  2007c).  Representative  concepts  with  signif¬ 
icant  enrichment  (p  <  0.001)  were  displayed  as  a  network  (Figure  4G;  Table  S4). 


Transcriptional 
Repression  of 
Target  Genes 


miR  Reporter  Luciferase  Assays 

The  50  bp  of  wild-type  or  mutant  3'  UTR  of  BMI1  and  RING2  containing  the 
predicted  miR-181a,b,  miR-200b,c  or  miR-203  binding  sites  (as  described  in 
Figures  S2C-S2F)  were  cloned  into  the  pMIR-REPORT  miRNA  Expression 
Reporter  Vector  (Ambion).  BT-549  cells  were  transfected  with  miRNAs  or 
controls  and  then  cotransfected  with  wild-type  3'  UTR-luc  or  mutant  3'  UTR- 
luc,  as  well  as  pRL-TK  vector  as  internal  control  for  luciferase  activity.  After 
48  hours  of  transfection,  the  cells  were  lysed  and  luciferase  assays  were  con¬ 
ducted  using  the  dual  luciferase  assay  system  (Promega,  Madison,  Wl).  Each 
experiment  was  performed  in  triplicate.  Drug  Treatment. 

BT-549  and  DU145  cells  were  treated  with  5  [xM  5-aza-2'-deoxycytidine 
(5-aza-dC)  for  6  days  (fresh  media  change  containing  the  drug  was  performed 


Prostate  Tumor  Xenograft  Model 

All  procedures  involving  mice  were  approved  by  the  University  Committee  on 
Use  and  Care  of  Animals  (UCUCA)  at  the  University  of  Michigan  and  conform 
to  their  relevant  regulatory  standards.  Five-week-old  male  nude  athymic 
BALB/c  nu/nu  mice  (Charles  River  Laboratory,  Wilmington,  MA)  were  used 
for  examining  tumorigenicity.  To  evaluate  the  role  of  BMI1  and  RING2  knock¬ 
down,  or  miR-181b,  miR-200b,  and  miR-203  overexpression  in  tumor  forma¬ 
tion,  the  DU145  stably  overexpressing  BMI1  shRNA,  RING2  shRNA,  scramble 
shRNA,  miR-181b,  miR-200b,  miR-203,  nontargeting  miR  or  vector  control 
cells  were  propagated  and  5  x  106  cells  were  inoculated  subcutaneously 
into  the  dorsal  flank  of  mice  (n  =  7  for  miR-203,  n  =  9  for  vector  control,  and 
n  =  8  for  Scramble,  BMI1  -sh3,  RING2-sh1 ,  miR-1 81  b,  miR-200b,  and  miR-NT, 
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respectively).  Tumor  size  was  measured  every  week,  and  tumor  volumes  were 
estimated  using  the  formula  (tc/6)  {L  x  W2),  where  L  =  length  of  tumor  and  1/1/  = 
width. 

Bisulfite  Modification  and  Methylation-Specific  PCR  of  miR-203 
in  Prostate  Tissues 

Bisulfite  conversion  was  carried  out  using  EZ  DNA  methylation  gold  kit  (Zymo 
Research  Corporation,  Orange,  CA)  according  to  manufacturer’s  instructions. 
Purified  DNA  (2  pil)  was  used  as  template  for  PCRs  with  primers  (Integrated 
DNA  Technologies  Inc.,  San  Diego,  CA)  and  synthesized  according  to  bisul¬ 
fite  converted  DNA  sequences  for  the  regions  of  interest  using  the  Meth- 
primer  software  (Li  and  Dahiya,  2002).  The  PCR  product  was  gel  purified 
and  cloned  into  pCR4  TOPO  TA  sequencing  vector  (Invitrogen,  Carlsbad, 
CA).  Plasmid  DNA  isolated  from  ten  colonies  from  each  sample  was 
sequenced  by  conventional  Sanger  Sequencing  (University  of  Michigan 
DNA  Sequencing  Core).  The  “BIQ  Analyzer”  (Bock  et  al.,  2005)  online  tool 
was  used  to  calculate  the  methylation  percentage  and  to  generate  the  bar 
graphs. 

ACCESSION  NUMBERS 

Coordinates  have  been  deposited  in  Gene  Expression  Omnibus  database  with 
accession  code  GSE26996. 

SUPPLEMENTAL  INFORMATION 

Supplemental  Information  includes  six  figures,  four  tables,  and  Supplemental 
Experimental  Procedures  and  can  be  found  with  this  article  online  at  doi:10. 
1 01 6/j.ccr.201 1.06.01 6. 
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Transcriptome  sequencing  across  a  prostate  cancer 
cohort  identifies  PCAT-1,  an  unannotated  lincRNA 
implicated  in  disease  progression 

John  R  Prensner1,8,  Matthew  K  Iyer1,8,  O  Alejandro  Balbin1,  Saravana  M  Dhanasekaran1,2,  Qi  Cao1, 

J  Chad  Brenner1,  Bharathi  Laxman3,  Irfan  A  Asangani1,  Catherine  S  Grasso1,  Hal  D  Kominsky1,  Xuhong  Cao1, 
Xiaojun  Jing1,  Xiaoju  Wang1,  Javed  Siddiqui1,  John  T  Wei4,  Daniel  Robinson1,  Hari  K  Iyer5, 

Nallasivam  Palanisamy1,2,6,  Christopher  A  Maher1,2  &  Arul  M  Chinnaiyan1,2, 4,6,7 


Noncoding  RNAs  (ncRNAs)  are  emerging  as  key  molecules  in  human  cancer,  with  the  potential  to  serve  as  novel  markers  of  disease 
and  to  reveal  uncharacterized  aspects  of  tumor  biology.  Here  we  discover  121  unannotated  prostate  cancer-associated  ncRNA 
transcripts  (PCATs)  by  ab  initio  assembly  of  high-throughput  sequencing  of  polyA+  RNA  (RNA-Seq)  from  a  cohort  of  102  prostate 
tissues  and  cells  lines.  We  characterized  one  ncRNA,  PCAT-1 ,  as  a  prostate-specific  regulator  of  cell  proliferation  and  show  that  it 
is  a  target  of  the  Polycomb  Repressive  Complex  2  (PRC2).  We  further  found  that  patterns  of  PCAT-1  and  PRC2  expression  stratified 
patient  tissues  into  molecular  subtypes  distinguished  by  expression  signatures  of  PCA T- J-repressed  target  genes.  Taken  together, 
our  findings  suggest  that  PCAT-1  is  a  transcriptional  repressor  implicated  in  a  subset  of  prostate  cancer  patients.  These  findings 
establish  the  utility  of  RNA-Seq  to  identify  disease-associated  ncRNAs  that  may  improve  the  stratification  of  cancer  subtypes. 


Recently,  RNA-Seq  has  provided  a  method  to  delineate  the  entire 
set  of  transcriptional  aberrations  in  a  disease,  including  novel  tran¬ 
scripts  not  measured  by  conventional  analyses1-5.  To  facilitate  inter¬ 
pretation  of  sequence  read  data,  existing  computational  methods 
typically  process  individual  samples  using  either  short  read  gapped 
alignment  followed  by  ab  initio  reconstruction2,3  or  de  novo  assembly 
of  read  sequences  followed  by  sequence  alignment4,5.  These  meth¬ 
ods  provide  a  powerful  framework  to  uncover  uncharacterized  RNA 
species,  including  antisense  transcripts,  short  RNAs  <250  bp  or  long 
intergenic  ncRNAs  (lincRNAs)  >250  bp. 

Although  still  largely  unexplored,  ncRNAs,  particularly  lincRNAs, 
have  emerged  as  a  new  aspect  of  biology,  with  evidence  suggesting  that 
they  are  frequently  cell- type  specific,  contribute  important  functions 
to  numerous  systems6,7  and  may  interact  with  known  cancer  genes 
such  as  EZH2  (ref.  8).  Indeed,  several  well-described  examples,  such 
as  HOTAIR 8,9  and  ANRIL10,11,  indicate  that  ncRNAs  maybe  essential 
actors  in  cancer  biology  typically  facilitating  epigenetic  gene  repres¬ 
sion  through  chromatin -modifying  complexes12,13.  Moreover,  ncRNA 
expression  may  confer  clinical  information  about  disease  outcomes 
and  have  utility  as  diagnostic  tests9,14.  The  characterization  of  RNA 
species,  their  functions  and  their  clinical  applicability  is  therefore  a 
major  area  of  biological  and  clinical  importance. 

Here,  we  describe  a  comprehensive  analysis  of  lincRNAs  in  102 
prostate  cancer  tissue  samples  and  cell  lines  by  RNA-Seq.  We  apply 


ab  initio  computational  approaches  to  delineate  the  annotated  and 
unannotated  transcripts  in  this  disease,  and  we  find  121  ncRNAs, 
termed  PCATs,  whose  expression  patterns  distinguish  benign,  local¬ 
ized  cancer  and  metastatic  cancer  samples.  Notably,  we  discover 
PCAT-1,  a  previously  undescribed  prostate  cancer  ncRNA  that  demon¬ 
strates  either  repression  by  PRC2  or  an  active  role  in  promoting  cell 
proliferation  through  transcriptional  regulation  of  target  genes.  To 
our  knowledge,  our  findings  describe  the  first  comprehensive  study 
of  lincRNAs  in  prostate  cancer,  provide  a  computational  framework 
for  large-scale  RNA-Seq  analyses  and  describe  PCAT-1  as  a  prostate 
cancer  ncRNA  functionally  implicated  in  disease  progression. 

RESULTS 

RNA-Seq  analysis  of  the  prostate  cancer  transcriptome 

Over  two  decades  of  research  have  generated  a  genetic  model  of 
prostate  cancer  based  on  numerous  neoplastic  events,  such  as  loss 
of  the  PTEN 15  tumor  suppressor  gene  and  gain  of  oncogenic  ETS 
family  transcription  factor  gene  fusions16-18  in  large  subsets  of  pros¬ 
tate  cancer  patients.  As  some  patients  lack  these  genetic  aberrations, 
we  hypothesized  that  prostate  cancer  similarly  harbored  disease- 
associated  ncRNAs  that  characterized  specific  molecular  subtypes. 

To  pursue  this  hypothesis,  we  applied  transcriptome  sequencing 
on  a  cohort  of  102  prostate  tissues  and  cell  lines — 20  benign  adjacent 
prostates  (benign),  47  localized  prostate  cancers  (PCA),  14  metastatic 
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Figure  1  Analysis  of  transcriptome  data  for 
the  detection  of  unannotated  transcripts. 

(a)  Schematic  overview  of  the  methodology 
employed  in  this  study,  (b)  Graphical 
representation  of  the  bioinformatics  filters  used 
to  merge  individual  transcriptome  libraries  into 
a  single  consensus  transcriptome.  The  merged 
consensus  transcriptome  was  generated  by 
compiling  all  individual  transcriptome  libraries 
and  using  individual  decision  tree  classifiers 
for  each  chromosome  to  define  high-confidence 
‘expressed’  transcripts  and  low-confidence 
‘background’  transcripts,  which  were  discarded. 
The  example  decision  tree  on  the  left  was 
trained  on  transcripts  on  chromosome  1.  The 
graphics  on  the  right  illustrate  the  application 
of  the  informatics  filtration  pipeline  to  sample 
assembled  transcripts,  (c)  After  informatic 
processing  and  filtration  of  the  sequencing 
data,  transcripts  were  categorized  to  identify 
unannotated  ncRNAs.  Transcribed  pseudogenes 
were  isolated,  and  the  remaining  transcripts  were 
categorized  based  on  overlap  with  an  aggregated 
set  of  known  gene  annotations  into  annotated 
protein  coding,  noncoding  and  unannotated. 

Both  annotated  and  unannotated  ncRNA 
transcripts  were  then  separated  into  intronic, 
intergenic  and  antisense  categories  based  on 
their  relationship  to  protein-coding  genes. 


tumors  and  21  prostate  cell  lines.  From  a  total 
of  1.723  billion  sequence  fragments  from 
201  lanes  of  sequencing  (108  paired-end 
and  93  single  reads  on  the  Illumina  Genome 
Analyzer  and  Genome  Analyzer  II),  we  per¬ 
formed  short-read  gapped  alignment19  and 
recovered  1.41  billion  mapped  reads,  with 
a  median  of  14.7  million  mapped  reads  per 
sample  (Supplementary  Table  1).  We  used  the 
Cufflinks  ab  initio  assembly  approach3  to  pro¬ 


duce,  for  each  sample,  the  most  probable  set  of 

putative  transcripts  that  served  as  the  RNA  templates  for  the  sequence 
fragments  in  that  sample  (Fig.  la  and  Supplementary  Figs.  1  and  2). 

As  expected  from  a  large  tumor  tissue  cohort,  individual  transcript 
assemblies  may  have  sources  of  noise,  such  as  artifacts  of  the  sequence 
alignment  process,  unspliced  intronic  pre-mRNA  and  genomic  DNA 
contamination.  To  exclude  these  from  our  analyses,  we  trained  a  deci¬ 
sion  tree  to  classify  transcripts  as  expressed  versus  background  on 
the  basis  of  transcript  length,  number  of  exons,  recurrence  in  mul¬ 
tiple  samples  and  other  structural  characteristics  (Fig.  lb,  left,  and 
Supplementary  Methods).  The  classifier  demonstrated  a  sensitiv¬ 
ity  of  70.8%  and  specificity  of  88.3%  when  trained  using  transcripts 
that  overlapped  genes  in  the  AceView  database20,  including  11.7% 
of  unannotated  transcripts  that  were  classified  as  expressed  (Fig.  lb 
right).  We  then  clustered  the  expressed  transcripts  into  a  consen¬ 
sus  transcriptome  and  applied  additional  heuristic  filters  to  further 
refine  the  assembly  (Supplementary  Methods).  The  final  ab  initio 
transcriptome  assembly  yielded  35,415  distinct  transcriptional  loci 
(Supplementary  Table  2  and  Supplementary  Methods). 

Discovery  of  prostate  cancer  noncoding  RNAs 

We  compared  the  assembled  prostate  cancer  transcriptome  to  the 
UCSC,  Ensembl,  RefSeq,  Vega  and  ENCODE  gene  databases  to  iden¬ 
tify  and  categorize  transcripts  (Fig.  lc).  The  majority  of  the  transcripts 


(77.3%)  corresponded  to  annotated  protein  coding  genes  (72.1%) 
and  noncoding  RNAs  (5.2%),  but  a  substantial  percentage  (19.8%) 
lacked  any  overlap  and  were  designated  unannotated  (Fig.  2a). 
These  included  partially  intronic  antisense  (2.44%),  totally  intronic 
(12.1%)  and  intergenic  transcripts  (5.25%),  consistent  with  previ¬ 
ous  reports  of  unannotated  transcription21-23.  Because  of  the  added 
complexity  of  characterizing  antisense  or  partially  intronic  tran¬ 
scripts  without  strand-specific  RNA-Seq  libraries,  we  focused  on 
totally  intronic  and  intergenic  transcripts. 

Global  characterization  of  unannotated  intronic  and  intergenic 
transcripts  demonstrated  that  they  were  more  highly  expressed 
(Fig.  2b),  had  greater  overlap  with  expressed  sequence  tags  (ESTs) 
(Supplementary  Fig.  3)  and  displayed  a  clear  but  subtle  increase  in 
conservation  over  randomly  permuted  controls  (intergenic  transcripts 
P  =  2.7  x  10-4  ±  0.0002  for  0.4  <  co  <  0.8;  intronic  transcripts  P  =  2.6  x 
10-5  +  0.0017  for  0  <  co  <  0.4,  Fishers  exact  test,  Fig.  2c).  By  contrast, 
unannotated  transcripts  scored  lower  than  protein-coding  genes  for 
these  metrics,  which  corroborates  data  in  previous  reports2,24.  Notably, 
a  small  subset  of  unannotated  intronic  transcripts  showed  a  profound 
degree  of  conservation  (Fig.  2c,  inset).  Finally,  analysis  of  coding 
potential  revealed  that  only  5  of  6,144  transcripts  harbored  a  high- 
quality  open  reading  frame  (ORF),  indicating  that  the  vast  majority  of 
these  transcripts  represent  ncRNAs  (Supplementary  Fig.  4). 
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Figure  2  Prostate  cancer  transcriptome  sequencing  reveals  dysregulation  of  unannotated  transcripts,  (a)  Global  overview  of  transcription  in  prostate 
cancer.  The  pie  chart  on  the  left  displays  transcript  distribution  in  prostate  cancer.  The  pie  charts  on  the  right  display  unannotated  (upper)  or  annotated 
(lower)  ncRNAs  categorized  as  sense  transcripts  (intergenic  and  intronic)  and  antisense  transcripts,  respectively,  (b)  Line  graph  showing  that  unannotated 
transcripts  are  more  highly  expressed  (reads  per  kilobase  of  transcript  per  million  mapped  reads;  RPKM)  than  control  regions.  Negative  control  intervals 
were  generated  by  randomly  permuting  the  genomic  positions  of  the  transcripts,  (c)  Conservation  analysis  comparing  unannotated  transcripts  to  known 
genes  and  intronic  controls  shows  a  subtle  degree  of  purifying  selection  among  unannotated  transcripts.  The  inset  on  the  right  shows  an  enlarged  view, 
(d-g)  Intersection  plots  displaying  the  fraction  of  unannotated  transcripts  enriched  for  H3K4me2  (d),  H3K4me3  (e),  acetyl-H3  (f)  or  RNA  polymerase  II 
(g)  at  their  transcriptional  start  site  (TSS)  using  ChIP-Seq  and  RNA-Seq  data  for  the  VCaP  prostate  cancer  cell  line.  The  legend  applies  to  plots  in  b-g.  (h) 
A  pie  chart  displaying  the  distribution  of  differentially  expressed  transcripts  in  prostate  cancer  (FDR  <  0.01). 


To  determine  whether  our  unannotated  transcripts  were  sup¬ 
ported  by  histone  modifications  defining  active  transcriptional  units, 
we  used  published  prostate  cancer  chromatin  immunoprecipitation 
(ChlP)-Seq  data  for  two  prostate  cell  lines25,  VCaP  and  LNCaP 
(Supplementary  Table  3).  After  filtering  our  data  set  for  transcribed 
repetitive  elements  known  to  display  alternative  patterns  of  histone 
modifications26,  we  observed  a  strong  enrichment  for  histone  modi¬ 
fications  characterizing  transcriptional  start  sites  (TSSs)  and  active 
transcription,  including  H3K4me2,  H3K4me3,  acetyl-H3  and  RNA 
polymerase  II  (Fig.  2d-g),  but  not  H3K4mel,  which  characterizes 
Sjjk  enhancer  regions27  (Supplementary  Figs.  5  and  6).  Notably,  inter- 
genic  ncRNAs  showed  greater  enrichment  compared  to  intronic 
ncRNAs  in  these  analyses  (Fig.  2d-g). 

To  elucidate  global  changes  in  transcript  abundance  in  prostate 
cancer,  we  analyzed  differential  expression  for  all  transcripts.  We 
found  836  genes  differentially  expressed  between  benign  samples  and 
localized  tumors  (false-discovery  rate  (FDR)  <  0.01),  with  annotated 
protein-coding  and  ncRNA  genes  constituting  82.8%  and  7.4%  of 
differentially  expressed  genes,  respectively,  including  known  pros¬ 
tate  cancer  biomarkers  such  AMACR28,  HPN29  and  PC  A3  (ref.  14) 
(Fig.  2h,  Supplementary  Fig.  2  and  Supplementary  Table  4). 
Finally,  9.8%  of  differentially  expressed  genes  corresponded  to 
unannotated  ncRNAs,  including  3.2%  within  gene  introns  and  6.6%  in 
intergenic  regions. 

Characterization  of  PCATs 

As  ncRNAs  may  contribute  to  human  disease6-9,  we  identified  aber¬ 
rantly  expressed  uncharacterized  ncRNAs  in  prostate  cancer.  We 
found  a  total  of  1,859  unannotated  lincRNAs  throughout  the  human 
genome.  Overall,  these  intergenic  RNAs  resided  approximately  half¬ 
way  between  two  protein  coding  genes  (Supplementary  Fig.  7),  and 
over  one-third  (34.1%)  were  >10  kb  from  the  nearest  protein-coding 


gene,  which  is  consistent  with  previous  reports30  and  supports  the 
independence  of  intergenic  ncRNAs  genes.  For  example,  visualizing 
the  Chrl5q  arm  using  the  Circos  program  (http://circos.ca/)  illus¬ 
trated  genomic  positions  of  89  unannotated  intergenic  transcripts, 
including  one  differentially  expressed  gene  centromeric  to  TLE3 
(Supplementary  Fig.  8). 

A  focused  analysis  of  the  1,859  unannotated  intergenic  RNAs 
yielded  106  that  were  differentially  expressed  in  localized  tumors  (FDR 
<  0.05,  Fig.  3a).  A  cancer  outlier  expression  analysis  (Supplementary 
Methods)  similarly  nominated  numerous  unannotated  ncRNA  out¬ 
liers  (Fig.  3b)  as  well  as  known  prostate  cancer  outliers,  such  as 
ERG18,  ETV1  (refs.  17,18),  SPINK1  (ref.  31)  and  CRISP3  (ref.  32). 
Merging  these  results  produced  a  set  of  121  unannotated  transcripts 
that  accurately  discriminated  benign,  localized  tumor  and  metastatic 
prostate  samples  by  unsupervised  clustering  (Fig.  3a).  Indeed,  clus¬ 
tering  analyses  using  unannotated  ncRNA  outliers  also  suggested 
disease  subtypes  (Supplementary  Fig.  9).  These  121  unannotated 
transcripts  were  ranked  and  named  as  PCATs  according  to  their 
fold-change  in  localized  tumor  versus  benign  tissue  (Supplementary 
Tables  5-7). 

Validation  of  novel  ncRNAs 

To  gain  confidence  in  our  transcript  nominations,  we  validated  mul¬ 
tiple  unannotated  transcripts  in  vitro  by  reverse  transcription  PCR 
(RT-PCR)  and  quantitative  real-time  PCR  (qPCR)  (Supplementary 
Fig.  10).  qPCR  for  four  transcripts  ( PCAT-114 ,  PCAT-14,  PCAT-43 
and  PCAT-1)  on  two  independent  cohorts  of  prostate  tissues  con¬ 
firmed  predicted  cancer-specific  expression  patterns  (Fig.  3c-f  and 
Supplementary  Fig.  11).  Notably,  all  four  are  prostate-specific,  with 
minimal  expression  seen  by  qPCR  in  breast  ( n  =  14)  or  lung  cancer 
[n  -  16)  cell  lines  or  in  19  normal  tissue  types  (Supplementary  Table  8). 
This  is  further  supported  by  expression  analysis  of  these  transcripts  in 
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Figure  3  Unannotated  intergenic  transcripts 
differentiate  prostate  cancer  and  benign 
prostate  samples,  (a)  Unsupervised  clustering 
analyses  of  differentially  expressed  or  outlier 
unannotated  intergenic  transcripts  clusters 
benign  samples,  localized  tumors  and 
metastatic  cancers.  Expression  is  plotted 
as  log2  fold-change  relative  to  the  median 
of  the  benign  samples.  The  four  transcripts 
detailed  in  this  study  are  indicated  on  the 
side,  (b)  Cancer  outlier  expression  analysis 
for  the  prostate  cancer  transcriptome 
ranks  unannotated  transcripts  prominently. 

(c-f)  qPCR  on  an  independent  cohort  of 
prostate  and  nonprostate  samples  (benign 
(n  =  19),  PCA  (n  =  35),  metastatic  (MET) 

(/?  =  31),  prostate  cell  lines  (a?  =  7),  breast  cell 
lines  [n  =  14),  lung  cell  lines  (n  =  16),  other 
normal  samples  (/?  =  19);  Supplementary 
Table  8))  measures  expression  levels  of  four 
nominated  ncRNAs — PCAT-14  (c),  PCAT-43  (d), 
PCAT-1 14  (e),  PCAT-1  (f) — and  upregulated 
in  prostate  cancer.  Inset  tables  on  the  right 
quantify  ‘positive’  and  ‘negative’  expressing 
samples  using  the  cut-off  value  (shown  as  a 
black  dashed  lines).  Statistical  significance 
was  determined  using  a  Fisher’s  exact  test. 
qPCR  analysis  was  performed  by  normalizing 
to  GAPDH  and  the  median  expression  of  the 
benign  samples. 


our  RNA-Seq  compendium  of  13  tumor  types, 
representing  325  samples  (Supplementary 
Fig.  12).  This  tissue  specificity  was  not  neces¬ 
sarily  due  to  regulation  by  androgen  receptor 
signaling,  as  only  PCAT-14  expression  was 

induced  when  androgen  responsive  VCaP  and  LNCaP  cells  were 
treated  with  the  synthetic  androgen  R1881,  consistent  with  previous 
data  from  this  locus17  (Supplementary  Fig.  13).  PCAT-1  and  PCAT- 
14  also  showed  cancer- specific  upregulation  when  tested  on  a  panel  of 
matched  tumor-normal  pair  samples  (Supplementary  Fig.  14). 

Of  note,  PCAT-1 14,  which  ranks  as  the  fifth  best  outlier,  just  ahead 
of  ERG  (Fig.  3b  and  Supplementary  Table  7),  appears  as  part  of 
a  large,  >500  kb  locus  of  expression  in  a  gene  desert  in  Chr2q31. 
We  termed  this  region  second  chromosome  locus  associated  with 
prostate-1’  (SChLAPl)  (Supplementary  Fig.  15).  Careful  analysis  of 
the  SChLAPl  locus  revealed  both  discrete  transcripts  and  intronic 
transcription,  highlighting  this  region  as  an  intriguing  aspect  of  the 
prostate  cancer  transcriptome. 

PCAT-1 ,  an  unannotated  prostate  cancer  lincRNA 

To  explore  several  transcripts  more  closely,  we  carried  out  5'  and  3' 
rapid  amplification  of  cDNA  ends  (RACE)  for  PCAT-1  and  PCAT-14. 
Interestingly,  the  PCAT-14  locus  contained  components  of  viral  ORFs 
from  the  HERV-K  endogenous  retrovirus  family  (Supplementary 
Fig.  16),  whereas  PCAT-1  incorporates  portions  of  a  mariner  family 
transposase33,34,  an  Alu  and  a  viral  long  terminal  repeat  promoter 
region  (Fig.  4a  and  Supplementary  Fig.  17).  Whereas  PCAT-14  was 
upregulated  in  localized  prostate  cancer  but  largely  absent  in  metas- 
tases  (Fig.  3c),  PCAT-1  was  strikingly  upregulated  in  a  subset  of  meta¬ 
static  and  high-grade  localized  (Gleason  score  >7)  cancers  (Fig.  3f 
and  Supplementary  Fig.  11).  Because  of  this  notable  profile,  we 
hypothesized  that  PCAT-1  may  have  coordinated  expression  with  the 
oncoprotein  EZH2,  a  core  PRC2  protein  that  is  upregulated  in  solid 
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tumors  and  contributes  to  a  metastatic  phenotype35,36.  Surprisingly, 
we  found  that  PCAT-1  and  EZH2  expression  were  nearly  mutually 
exclusive  (Fig.  4b),  with  only  one  patient  showing  outlier  expression 
of  both.  This  suggests  that  outlier  PCAT-1  and  EZH2  expression  may 
define  two  subsets  of  high-grade  disease. 

PCAT-1  is  located  in  the  chromosome  8q24  gene  desert  -725  kb 
upstream  of  the  c-MYC  oncogene.  To  confirm  that  PCAT-1  is  a 
noncoding  gene,  we  cloned  the  full-length  PCAT-1  transcript  and 
performed  in  vitro  translational  assays,  which  were  negative  as  expected 
(Supplementary  Fig.  18).  Next,  because  Chr8q24  is  known  to  harbor 
prostate  cancer-associated  single  nucleotide  polymorphisms  (SNPs) 
and  to  exhibit  frequent  chromosomal  amplification37-42,  we  evaluated 
whether  the  relationship  between  EZH2  and  PCAT-1  was  specific  or 
generalized.  To  address  this,  we  measured  expression  levels  of  c-MYC 
and  NCOA2 ,  two  proposed  targets  of  Chr8q  amplification39,42,  by 
qPCR.  Neither  c-MYC  nor  NCOA2  levels  showed  striking  expression 
relationships  to  PCAT-1,  EZH2  or  each  other  (Supplementary  Fig.  19). 
Likewise,  PCAT-1  outlier  expression  was  not  dependent  on  Chr8q24 
amplification,  as  highly  expressing  localized  tumors  often  did  not  have 
8q24  amplification  and  high  copy  number  gain  of  8q24  was  not  suf¬ 
ficient  to  upregulate  PCAT-1  (Supplementary  Figs.  20  and  21). 

PCAT-1  function  and  regulation 

Despite  reports  showing  that  upregulation  of  the  ncRNA  HOTAIR 
participates  in  PRC2  function  in  breast  cancer9,  we  do  not  observe 
strong  expression  of  this  ncRNA  in  prostate  (Supplementary  Fig.  22), 
suggesting  that  other  ncRNAs  may  be  important  in  this  cancer.  To 
determine  the  mechanism  for  the  expression  profiles  of  PCAT-1  and 
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Figure  4  PCAT-1  is  a  marker  of  aggressive  cancer 
and  a  PRC2-repressed  ncRNA.  (a)  The  genomic 
location  of  PCAT-1  determined  by  5'  and  3'  RACE, 
with  DNA  sequence  features  indicated  by  the 
colored  boxes,  (b)  qPCR  for  PCAT-1  (yaxis)  and 
EZH2  (x axis)  on  a  cohort  of  benign  [n=  19), 
localized  tumor  (n  =  35)  and  metastatic  cancer 
(n  =  31)  samples.  The  inset  table  quantifies 
patient  subsets  demarcated  by  the  gray  dashed 
lines,  (c)  Knockdown  of  EZH2  in  VCaP  resulted 
in  upregulation  of  PCAT-1.  Data  were  normalized 
to  GAPDH  and  represented  as  fold-change. 

ERG  and  B-actin  serve  as  negative  controls.  The 
inset  western  blot  indicates  EZH2  knockdown. 

(d)  Treatment  of  VCaP  cel  Is  with  0. 1  jtM  of  the 
EZH2  inhibitor  DZNep  or  vehicle  control  (DMSO) 
shows  increased  expression  of  PCAT-1  transcript 
after  EZH2  inhibition,  (e)  PCAT-1  expression 
is  increased  upon  treatment  of  VCaP  cells  with 
the  demethylating  agent  5'azacytidine  (5'Aza), 
the  histone  deacetylase  inhibitor  SAHA  or  a 
combination  of  both.  qPCR  data  were  normalized  to 
the  average  of  (GAPDH  +  (3-actin)  and  represented 
as  fold-change.  GSTP1  and  FKBP5  are  positive  and 
negative  controls,  respectively,  (f)  ChIP  assays  for 
SUZ12  demonstrated  direct  binding  of  SUZ12  to 
the  PCAT-1  promoter.  Primer  locations  are  indicated 
(boxed  numbers)  in  the  PCAT-1  schematic. 
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EZH2 ,  we  inhibited  EZH2  activity  in  VCaP  cells,  which  express  low- 
to-moderate  levels  of  PCAT-1.  Knockdown  of  EZH2  by  short  hairpin 
(sh)RNA  or  pharmacologic  inhibition  of  EZH2  with  the  inhibitor 
3-deazaneplanocin  A  (DZNep)  caused  a  dramatic  upregulation  in 
PCAT-1  expression  levels  (Fig.  4c, d),  as  did  treatment  of  VCaP  cells 
with  the  demethylating  agent  5'deoxyazacytidine,  the  histone  deacety¬ 
lase  inhibitor  SAHA  or  both  (Fig.  4e).  ChIP  assays  also  demonstrated 
that  SUZ12,  a  core  PRC2  protein,  directly  binds  the  PCAT-1  promoter 
~1  kb  upstream  of  the  TSS  (Fig.  4f).  Notably,  RNA  immunoprecipita- 
tion  similarly  showed  binding  of  PCAT-1  to  SUZ12  protein  in  VCaP 
cells  (Supplementary  Fig.  23a).  RNA  immunoprecipitation  assays 
followed  by  RNase  A,  RNase  H  or  DNase  I  treatment  either  abolished, 
partially  preserved  or  totally  preserved  this  interaction,  respectively 
(Supplementary  Fig.  23b).  This  suggests  that  PCAT-1  exists  primarily 
as  a  single- stranded  RNA  and  secondarily  as  a  RNA/DNA  hybrid. 

To  explore  the  functional  role  of  PCAT-1  in  prostate  cancer,  we  stably 
overexpressed  full-length  PCAT-1  or  controls  in  RWPE  benign  immor¬ 
talized  prostate  cells.  We  observed  a  modest  but  consistent  increase 
in  cell  proliferation  when  PCAT-1  was  overexpressed  at  physiological 


levels  (Fig.  5a  and  Supplementary  Fig.  24).  Next,  we  designed  short 
interfering  (si)RNA  oligos  to  PCAT-1  and  performed  knockdown  exper¬ 
iments  in  LNCaP  cells,  which  express  higher  levels  of  PCAT-1  without 
PRC2-mediated  repression  (Supplementary  Fig.  25).  Supporting  our 
overexpression  data,  knockdown  of  PCAT-1  with  three  independent 
siRNA  oligos  resulted  in  a  25-50%  decrease  in  cell  proliferation  in 
LNCaP  cells  (Fig.  5b),  but  not  in  control  DU145  cells  lacking  PCAT-1 
expression  (Supplementary  Fig.  26)  or  VCaP  cells,  in  which  PCAT-1 
is  expressed  but  repressed  by  PRC2  (Supplementary  Fig.  27). 

Gene  expression  profiling  of  LNCaP  knockdown  samples  on 
cDNA  microarrays  indicated  that  PCAT-1  modulates  the  transcrip¬ 
tional  regulation  of  370  genes  (255  upregulated,  115  downregulated; 
FDR  <0.01)  (Supplementary  Fig.  28  and  Supplementary  Table  9). 
Gene  ontology  analysis  of  the  upregulated  genes  showed  preferen¬ 
tial  enrichment  for  gene  set  concepts  such  as  mitosis  and  cell  cycle, 
whereas  the  downregulated  genes  had  no  concepts  showing  statistical 
significance  (Fig.  5c  and  Supplementary  Table  10).  These  results 
suggest  that  the  function  of  PCAT-1  is  predominantly  repressive 
in  nature,  similar  to  other  lincRNAs.  We  next  validated  expression 
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Figure  5  PCAT-1  promotes  cell  proliferation,  (a)  Cell  proliferation  assays  for  RWPE  benign  immortalized  prostate  cells  stably  infected  with  PCAT-1 
lentivirus  or  RFP  and  LacZ  control  lentiviruses.  An  asterisk  (*)  indicates  P<  0.02  by  a  two-tailed  Student’s  Mest.  (b)  Cell  proliferation  assays  in  LNCaP 
using  PCAT-1  siRNAs.  An  asterisk  (*)  indicates  P<  0.005  by  a  two-tailed  Student’s  t- test,  (c)  Gene  ontology  analysis  of  PCAT-1  knockdown  microarray 
data  using  the  DAVID  program.  Blue  bars  represent  the  top  hits  for  upregulated  genes.  Red  bars  represent  the  top  hits  for  downregulated  genes. 

DAVID  enrichment  scores  are  represented  with  Benjamini-Hochberg-adjusted  P  values.  All  error  bars  in  this  figure  are  mean  ±  s.e.m. 
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Figure  6  Prostate  cancer  tissues  recapitulate  PCAT-1  signaling,  (a)  qPCR  expression  of  three  PCAT-1  target  genes  after  PCAT-1  knockdown  in  VCaP  and 
LNCaP  cells,  as  well  as  following  EZH2  knockdown  or  dual  EZH2  and  PCAT-1  knockdown  in  VCaP  cells.  qPCR  data  were  normalized  to  the  average  of 
{GAPDH  +  (3-actin)  and  represented  as  fold  change.  Error  bars  represent  mean  ±  s.e.m.  (b)  Standardized  log2-transformed  qPCR  expression  of  a  set 
of  tumors  and  metastases  with  outlier  expression  of  either  PCAT-1  or  EZH2.  The  shaded  squares  in  the  lower  left  show  Spearman  correlation  values 
between  the  indicated  genes  (*  indicates  P<  0.05).  Blue  and  red  indicate  negative  or  positive  correlation,  respectively.  The  upper  squares  show  the 
scatter  plot  matrix  and  fitted  trend  lines  for  the  same  comparisons,  (c)  A  heatmap  of  PCAT-1  target  genes  ( BRCA2 ,  CENPF,  CENPE)  in  EZf/Z-outlier 
and  PCAT-I-outlier  patient  samples  (see  Fig.  4b).  Expression  was  determined  by  qPCR  and  normalized  as  in  b.  (d)  A  predicted  network  generated  by 
the  HefaLMP  program  for  7  of  20  top  upregulated  genes  following  PCAT-1  knockdown  in  LNCaP  cells.  Gray  nodes  are  genes  found  following  PCAT-1 
knockdown.  Red  edges  indicate  co-expressed  genes;  black  edges  indicate  predicted  protein-protein  interactions;  and  purple  edges  indicate  verified 
protein-protein  interactions,  (e)  A  proposed  schematic  representing  PCAT-1  upregulation,  function  and  relationship  to  PRC2. 


changes  in  three  key  PCAT-1  target  genes  ( BRCA2 ,  CENPE  and 
CENPF)  whose  expression  is  upregulated  upon  PCAT-1  knockdown 
(Fig.  6a)  in  LNCaP  and  VCaP  cells,  the  latter  of  which  appear  less 
sensitive  to  PCAT-1  knockdown  likely  due  to  lower  overall  expression 
levels  of  this  transcript. 

PCAT-1  signatures  in  prostate  cancer 

Because  of  the  regulation  of  PCAT-1  by  PRC2  in  VCaP  cells,  we 
hypothesized  that  knockdown  of  EZH2  would  also  downregulate 
PCAT-1  targets  as  a  secondary  phenomenon  owing  to  the  subsequent 
upregulation  of  PCAT-1.  Simultaneous  knockdown  of  PCAT-1  and 
EZH2  would  thus  abrogate  expression  changes  in  PCAT-1  target 
genes.  Carrying  out  this  experiment  in  VCaP  cells  demonstrated  that 
PCAT-1  target  genes  were  indeed  downregulated  by  EZH2  knock¬ 
down,  and  that  this  change  was  either  partially  or  completely  reversed 
using  siRNA  oligos  to  PCAT-1  (Fig.  6a),  lending  support  to  the  role 
of  PCAT-1  as  a  transcriptional  repressor.  Taken  together,  these  results 
suggest  that  PCAT-1  biology  may  exhibit  two  distinct  modalities:  one 
in  which  PRC2  represses  PCAT-1  and  a  second  in  which  active  PCAT-1 
promotes  cell  proliferation.  PCAT-1  and  PRC2  may  therefore  charac¬ 
terize  distinct  subsets  of  prostate  cancer. 

To  examine  these  findings,  we  used  qPCR  to  measure  expres¬ 
sion  of  BRCA2 ,  CENPE  and  CENPF  in  our  cohort  of  tissue  samples. 
Consistent  with  our  model,  we  found  that  samples  expressing  PCAT-1 
tended  to  have  low  expression  of  PCAT-1  target  genes  (Fig.  6b). 
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Moreover,  comparing  EZH2- outlier  and  PCAT-1  -outlier  patients 
(Fig.  4b),  we  found  that  two  distinct  phenotypes  emerged.  Individuals 
with  high  EZH2  tended  to  have  high  levels  of  PCAT-1  target  genes, 
and  those  with  high  expression  of  PCAT-1  itself  displayed  the  opposite 
expression  pattern  of  target  genes  (Fig.  6c).  Network  analysis  of  the 
top  20  upregulated  genes  after  PCAT-1  knockdown  with  the  HefaLMP 
tool43  further  suggested  that  these  genes  form  a  coordinated  network 
(Fig.  6d),  corroborating  our  previous  observations.  Taken  together, 
these  results  provide  initial  data  into  the  composition  and  function 
of  the  prostate  cancer  ncRNA  transcriptome. 

DISCUSSION 

To  our  knowledge,  this  study  represents  the  largest  RNA-Seq  analysis 
to  date  and  the  first  to  comprehensively  analyze  a  common  epithelial 
cancer  from  a  large  cohort  of  human  tissue  samples.  As  such,  our  study 
has  adapted  existing  computational  tools  intended  for  small-scale  use3 
and  developed  new  methods  to  distill  large  numbers  of  transcrip¬ 
tome  data  sets  into  a  single  consensus  transcriptome  assembly  that 
accurately  represents  disease  biology  (Supplementary  Discussion). 

Among  the  numerous  uncharacterized  ncRNA  species  detected 
by  our  study,  we  have  focused  on  121  PCATs,  which  we  believe  rep¬ 
resent  a  set  of  uncharacterized  ncRNAs  that  may  have  important 
biological  functions  in  this  disease.  In  this  regard,  these  data  con¬ 
tribute  to  a  growing  body  of  literature  supporting  the  importance  of 
unannotated  ncRNA  species  in  cellular  biology  and  oncogenesis6-12, 
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and  broadly  our  study  confirms  the  utility  of  RNA-Seq  in  defining 
functionally  important  elements  of  the  genome2-4. 

Of  particular  interest  is  our  discovery  of  the  prostate-specific 
ncRNA  gene  PCAT-1 ,  which  is  markedly  overexpressed  in  a  subset 
of  prostate  cancers,  particularly  metastases,  and  may  contribute  to  cell 
proliferation  in  these  tumors.  It  is  also  notable  that  PCAT-1  resides  in 
the  8q24  gene  desert’  locus,  in  the  vicinity  of  well-studied  prostate 
cancer  risk  SNPs  and  the  c-MYC  oncogene,  suggesting  that  this  locus — 
and  its  frequent  amplification  in  cancer — may  be  linked  to  additional 
aspects  of  cancer  biology  (Supplementary  Discussion).  In  addition, 
the  interplay  between  PRC2  and  PCAT-1  further  suggests  that  this 
ncRNA  may  have  an  important  role  in  prostate  cancer  progression 
(Fig.  6e).  Other  ncRNAs  identified  by  this  analysis  may  similarly 
contribute  to  prostate  cancer  as  well.  Furthermore,  recent  preclinical 
efforts  to  detect  prostate  cancer  noninvasively  through  the  collection 
of  patient  urine  samples  have  shown  promise  for  several  urine-based 
prostate  cancer  biomarkers,  including  the  ncRNA  PC  A3  (refs.  44,45). 
Although  additional  studies  are  needed,  our  identification  of  ncRNA 
biomarkers  for  prostate  cancer  suggests  that  urine-based  assays  for 
these  ncRNAs  may  also  warrant  investigation,  particularly  for  those 
that  may  stratify  patient  molecular  subtypes. 

Our  findings  support  an  important  role  for  tissue-specific  ncRNAs 
in  prostate  cancer  and  suggest  that  cancer-specific  functions  of  these 
ncRNAs  may  help  to  drive  tumorigenesis.  We  further  speculate  that 
specific  ncRNA  signatures  may  occur  universally  in  all  disease  states 
and  that  applying  these  methodologies  to  other  diseases  may  reveal 
key  aspects  of  disease  biology  and  clinically  important  biomarkers. 

METHODS 

Methods  and  any  associated  references  are  available  in  the  online  version 
of  the  paper  at  http://www.nature.com/naturebiotechnology/. 

Accession  codes.  Data  from  RNA-Seq  experiments  are  deposited 
at  the  NCBI  Gene  Expression  Omnibus  as  GSE25183.  PCAT-1  and 
PC  AT- 14  nucleotide  sequences  are  deposited  at  GenBank  nucleotide 
database  (nuccore)  as  HQ605084  and  HQ605085,  respectively. 

Note:  Supplementary  information  is  available  on  the  Nature  Biotechnology  website. 
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ONLINE  METHODS 

Cell  lines,  treatments  and  tissues.  All  prostate  cell  lines  were  obtained  from 
the  American  Type  Culture  Collection,  except  for  PrEC  (benign  nonimmor- 
talized  prostate  epithelial  cells)  and  PrSMC  (prostate  smooth  muscle  cells), 
which  were  obtained  from  Lonza.  Cell  lines  were  maintained  using  standard 
media  and  conditions. 

For  androgen  treatment  experiments,  LNCaP  and  VCaP  cells  were  grown  in 
androgen- depleted  media  for  48  h  and  subsequently  treated  with  5nM  methyl- 
trienolone  (R1881,  NEN  Life  Science  Products)  or  an  equivalent  volume  of 
ethanol  for  48  h  before  harvesting  the  cells.  For  drug  treatments,  VCaP  cells 
were  treated  with  20  jlM  5'deoxyazacytidine  (Sigma),  500  nM  HDAC  inhibitor 
suberoylanilide  hydroxamic  acid  (SAHA)  (Biovision),  or  both  5'deoxyazacyti- 
dine  and  SAHA.  5'deoxyazacytidine  treatments  were  performed  for  6  d  with 
media  and  drug  reapplied  every  48  h.  SAHA  treatments  were  done  for  48  h. 
DMSO  treatments  were  done  for  6  d.  For  DZNep  treatments,  DZNep  was  dis¬ 
solved  in  DMSO  and  VCaP  cells  were  treated  with  either  0. 1  juM  of  DZNep  or 
vehicle  control;  RNA  was  harvested  at  72  h  and  144  h. 

Prostate  tissues  were  obtained  from  the  radical  prostatectomy  series  and 
Rapid  Autopsy  Program_ENREF_48  at  the  University  of  Michigan  tissue  core 
as  part  of  the  University  of  Michigan  Prostate  Cancer  Specialized  Program 
of  Research  Excellence  (S.P.O.R.E.).  All  tissue  samples  were  collected  with 
informed  consent  under  an  Institutional  Review  Board  (IRB)  approved  pro¬ 
tocol  at  the  University  of  Michigan. 

RNA  isolation,  cDNA  synthesis  and  PCR  experiments.  Total  RNA  was 
isolated  using  Trizol  and  an  RNeasy  Kit  (Invitrogen)  with  DNase  I  digestion 
according  to  the  manufacturer’s  instructions.  RNA  integrity  was  verified  on 
an  Agilent  Bioanalyzer  2100  (Agilent  Technologies).  cDNA  was  synthesized 
from  total  RNA  using  Superscript  III  (Invitrogen)  and  random  primers 
(Invitrogen).  Quantitative  Real-time  PCR  (qPCR)  was  done  using  Power 
SYBR  Green  Mastermix  (Applied  Biosystems)  on  an  Applied  Biosystems 
7900HT  Real-Time  PCR  System.  (RT-PCR  was  done  with  Platinum  Taq  High 
Fidelity  polymerase  (Invitrogen).  All  oligonucleotide  primers  are  listed  in 
Supplementary  Table  11.  For  PCR  product  sequencing,  PCR  products  were 
resolved  on  a  1.5%  agarose  gel,  and  either  sequenced  directly  or  extracted 
using  a  Gel  Extraction  kit  (Qiagen)  and  cloned  into  pcr4-TOPO  vectors 
(Invitrogen).  PCR  products  were  bidirectionally  sequenced  at  the  University 
of  Michigan  Sequencing  Core. 

RNA-ligase-mediated  rapid  amplification  of  cDNA  ends  (RACE).  5'  and 

3'  RACE  was  performed  using  the  GeneRacer  RLM-RACE  kit  (Invitrogen) 
according  to  the  manufacturer’s  instructions.  RACE  PCR  products  were 
obtained  using  Platinum  Taq  high-fidelity  polymerase  (Invitrogen),  the  sup¬ 
plied  GeneRacer  primers,  and  appropriate  gene- specific  primers  indicated  in 

Supplementary  Table  11. 

RNA-Seq  library  preparation.  2  jig  total  RNA  was  selected  for  polyA+  RNA 
using  Sera-Mag  oligo(dT)  beads  (Thermo  Scientific),  and  paired-end  next- 
generation  sequencing  libraries  were  prepared,  as  previously  described46, 
using  Illumina- supplied  universal  adaptor  oligos  and  PCR  primers  (Illumina). 
Samples  were  sequenced  in  a  single  lane  on  an  Illumina  Genome  Analyzer  I  or 
Genome  Analyzer  II  flow  cell  using  previously  described  protocols46.  36-45  mer 
paired-end  reads  were  done  according  to  the  protocol  provided  by  Illumina. 

Overexpression  studies.  PCAT-1  full-length  transcript  was  cloned  into  the 
pLenti6  vector  (Invitrogen)  along  with  RFP  and  LacZ  controls.  After  confir¬ 
mation  of  the  insert  sequence,  lentiviruses  were  generated  at  the  University  of 
Michigan  Vector  Core  and  transfected  into  the  benign  immortalized  prostate 
cell  line  RWPE.  RWPE  cells  stably  expressing  PCAT-1,  RFP  or  LacZ  were  gen¬ 
erated  by  selection  with  blasticidin  (Invitrogen),  and  10,000  cells  were  plated 
into  12-well  plates.  Cells  were  harvested  and  counted  at  day  2,  day  4  and  day 
6  post-plating  with  a  Coulter  counter. 


siRNA  knockdown  studies.  Cells  were  plated  and  transfected  with  20  JlM 
experimental  siRNA  oligos  or  nontargeting  controls  twice,  at  12  h  and  36  h 
post-plating.  Knockdowns  were  performed  with  Oligofectamine  in  OptiMEM 
media.  Knockdown  efficiency  was  determined  by  qPCR.  siRNA  sequences 
(in  sense  format)  for  PCAT-1  knockdown  were  as  follows:  siRNA  1  UU 
AAAGAGAUCCACAGUUAUU;  siRNA  2  GCAGAAACACCAAUGGAUA 
UU;  siRNA  3  AUACAUAAGACC AUGGAA AU ;  siRNA  4  GAACCUAACUGG 
ACUUUAAUU.  For  EZH2  siRNA,  the  following  sequence  was  used:  GAGG 
UUCAGACGAGCUGAUUU. 

shRNA  knockdown  and  western  blot  analysis.  Cells  were  seeded  at  50-60% 
confluency,  incubated  overnight,  and  transfected  with  EZH2  or  nontargeting 
shRNA  lentiviral  constructs  as  described  in  for  48  h.  GFP+  cells  were  drug- 
selected  using  1  Jlg/ml  puromycin.  RNA  and  protein  were  harvested  for  PCR 
and  western  blot  analysis  according  to  standard  protocols.  For  western  blot 
analysis,  PVDF  membranes  (GE  Healthcare)  were  incubated  overnight  at  4  °C 
with  either  EZH2  mouse  monoclonal  (1:1,000,  BD  Biosciences,  no.  612666), 
or  B-actin  (Abeam,  ab8226)  for  equal  loading. 

Gene  expression  profiling.  Agilent  Whole  Human  Genome  Oligo  Microarray 
was  used  for  cDNA  profiling  of  PCAT-1  siRNA  knockdown  samples  or  nontar¬ 
geting  control  according  to  standard  protocols_ENREF_50.  All  samples  were 
run  in  technical  triplicates  against  nontargeting  control  siRNA.  Expression 
array  data  was  processed  using  the  SAM  method47  with  an  FDR  <0.01.  Up- 
and  downregulated  probes  were  separated  and  analyzed  using  the  DAVID 
bioinformatics  platform48. 

ChIP.  Assays  were  done  as  previously  described25,  where  4-7  jig  of  the 
following  antibodies  were  used:  IgG  (Millipore,  PP64),  SUZ12  (Cell  Signaling, 
no.  3737)  and  SUZ12  (Abeam,  ab  12073).  ChIP-PCR  reactions  were  done  in 
triplicate  with  SYBRGreen  using  1:150th  of  the  ChIP  product  per  reaction. 

In  vitro  translation.  Full-length  PCAT-1,  Halo-tagged  ERG  or  GUS  positive 
control  were  cloned  into  the  PCR2.1  entry  vector  (Invitrogen)  and  in  vitro 
translational  assays  were  done  using  the  TnT  Quick  Coupled  Transcription/ 
Translation  System  (Promega)  with  1  mM  methionine  and  Transcend  Biotin- 
Lysyl-tRNA  (Promega)  according  to  the  manufacturer’s  instructions. 

Bioinformatic  analyses.  Sequencing  reads  were  aligned  with  TopHat19,  and 
ab  initio  assembly  was  performed  with  Cufflinks3.  Transcriptome  librar¬ 
ies  were  merged  and  statistical  classifiers  were  developed  and  employed  to 
filter  low-confidence  transcripts.  Nominated  transcripts  were  compared  to 
UCSC,  RefSeq,  Vega,  Ensembl  and  ENCODE  database,  and  coding  poten¬ 
tial  was  determined  with  the  txCdsPredict  program  from  UCSC.  Transcript 
conservation  was  determined  with  the  SiPhy  package.  Differential  expression 
analysis  was  performed  using  SAM  methodology,  and  outlier  analysis  using  a 
modified  COPA  method.  See  the  Supplementary  Methods  for  details  on  the 
bioinformatics  methods  used. 

Statistical  analyses  for  experimental  studies.  All  data  are  presented  as 
means  ±  s.e.m.  All  experimental  assays  were  performed  in  duplicate  or 
triplicate.  Statistical  analyses  shown  in  figures  represent  Fisher’s  exact  tests 
or  two-tailed  Student  f-tests,  as  indicated.  For  details  regarding  the  statis¬ 
tical  methods  employed  during  RNA-Seq  and  ChIP-Seq  data  analysis,  see 
Supplementary  Methods. 

46.  Maher,  C.A.  et  ai.  Chimeric  transcript  discovery  by  paired-end  transcriptome 
sequencing.  Proc.  Natl.  Acad.  Sci.  USA  106,  12353-12358  (2009). 

47. Tusher,  V.G.,  Tibshirani,  R.  &  Chu,  G.  Significance  analysis  of  microarrays  applied 
to  the  ionizing  radiation  response.  Proc.  Natl.  Acad.  Sci.  USA  98,  5116-5121 
(2001). 

48.  Dennis,  G.  Jr.  et  al.  DAVID:  Database  for  Annotation,  Visualization,  and  Integrated 
Discovery.  Genome  Biol.  4,  3  (2003). 
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Editor's  Summary 

Targeting  Outside  the  Box 

Out-of-the-box  thinking  is  highly  valued  in  all  creative  endeavors,  and  science  is  no  exception.  Similarly, 
out-of-the-cell,  or  extracellular,  drug  targets  have  many  advantages  over  intracellular  ones,  such  as  easy  access  by 
small-molecule  inhibitors  and  antibodies.  Because  one-third  of  all  cases  of  prostate  cancer — one  of  the  most 
prevalent  forms  of  the  disease  in  men  — are  aggressive  and  fast-growing,  and  traditional  treatments  are  often 
unhelpful  and  cause  troublesome  side  effects,  it  is  clear  that  some  out-of-the-box  thinking  is  required  to  address  this 
therapeutic  dilemma.  Now,  Ateeq  et  al.  have  identified  SPINK1  (serine  peptidase  inhibitor,  Kazal  type  1)  as  an 
extracellular  therapeutic  target  for  an  aggressive  subset  of  SPINK1  +  prostate  cancer.  SPINK1  is  highly  expressed  in 
~10%  of  prostate  cancers,  and  expression  has  been  correlated  with  aggressive  disease.  In  the  new  work,  the  authors 
showed  directly  that  SPINK1  contributes  to  the  aggressive  phenotype.  Forced  expression  of  recombinant  SPINK1 
increased  prostate  cancer  cell  proliferation  and  invasiveness,  whereas  knockdown  of  SPINK1  gene  expression  or 
treatment  with  a  SPINK1  -directed  monoclonal  antibody  resulted  in  decreased  cell  division,  invasiveness,  and  tumor 
growth.  Moreover,  SPINK1  mediated  its  neoplastic  effects  in  part  through  interactions  with  the  epidermal  growth  factor 
receptor  (EGFR).  Indeed,  antibodies  to  both  SPINK1  and  EGFR  blocked  the  growth  of  SPINK1+/ETS  -  tumors  more 
than  either  antibody  alone  and  did  not  affect  SPINK1  -  tumors.  Together,  these  findings  suggest  that  SPINK1 
represents  a  new,  specific,  and  — by  virtue  of  its  outside-of-the-box  location — druggable  target  for  a  potentially  lethal 
form  of  prostate  cancer. 
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CANCER 

Therapeutic  Targeting  of  SPINK1 -Positive 
Prostate  Cancer 

Bushra  Ateeq,1'2  Scott  A.  Tomlins,1'2  Bharathi  Laxman,1'2  Irfan  A.  Asangani,1'2  Qi  Cao,1'2 
Xuhong  Cao,1'3  Yong  Li,1'2  Xiaoju  Wang,1'2  Felix  Y.  Feng,1'4'5  Kenneth  J.  Pienta,1'5'6 
Sooryanarayana  Varambally,1'2  5  Arul  M.  Chinnaiyan1'2'3'5'7* 

Gene  fusions  involving  ETS  (erythroblastosis  virus  E26  transformation-specific)  family  transcription  factors  are 
found  in  -50%  of  prostate  cancers  and  as  such  can  be  used  as  a  basis  for  the  molecular  subclassification  of 
prostate  cancer.  Previously,  we  showed  that  marked  overexpression  of  SPINK1  ( serine  peptidase  inhibitor ,  Kazal 
type  7),  which  encodes  a  secreted  serine  protease  inhibitor,  defines  an  aggressive  molecular  subtype  of  ETS 
fusion-negative  prostate  cancers  ( SPINK1+/ETS~ ,  -10%  of  all  prostate  cancers).  Here,  we  examined  the  potential 
of  SPINK1  as  an  extracellular  therapeutic  target  in  prostate  cancer.  Recombinant  SPINK1  protein  (rSPINKI)  stim¬ 
ulated  cell  proliferation  in  benign  RWPE  as  well  as  cancerous  prostate  cells.  Indeed,  RWPE  cells  treated  with 
either  rSPINKI  or  conditioned  medium  from  22RV1  prostate  cancer  cells  ( SPINK1+/ETS~ )  significantly  increased 
cell  invasion  and  intravasation  when  compared  with  untreated  cells.  In  contrast,  knockdown  of  SPINK1  in  22RV1 
cells  inhibited  cell  proliferation,  cell  invasion,  and  tumor  growth  in  xenograft  assays.  22RV1  cell  proliferation, 
invasion,  and  intravasation  were  attenuated  by  a  monoclonal  antibody  (mAb)  to  SPINK1  as  well.  We  also  demon¬ 
strated  that  SPINK1  partially  mediated  its  neoplastic  effects  through  interaction  with  the  epidermal  growth  factor 
receptor  (EGFR).  Administration  of  antibodies  to  SPINK1  or  EGFR  (cetuximab)  in  mice  bearing  22RV1  xenografts 
attenuated  tumor  growth  by  more  than  60  and  40%,  respectively,  or  -75%  when  combined,  without  affecting 
PC3  xenograft  ( SPINK1~/ETS~ )  growth.  Thus,  this  study  suggests  that  SPINK1  may  be  a  therapeutic  target  in  a 
subset  of  patients  with  SPINK1+/ETS~  prostate  cancer.  Our  results  provide  a  rationale  for  both  the  development 
of  humanized  mAbs  to  SPINK1  and  evaluation  of  EGFR  inhibition  in  SPINK1+/ETS~  prostate  cancers. 


INTRODUCTION 

Therapies  targeted  against  specific  molecular  alterations  present  only 
in  cancer  cells  have  revolutionized  the  treatment  of  several  cancers. 
For  example,  targeting  ERBB2,  which  is  amplified  in  -20%  of  breast 
cancers,  with  the  humanized  monoclonal  antibody  (mAb)  trastuzumab 
(Herceptin)  has  resulted  in  improved  survival  for  breast  cancer  pa¬ 
tients.  Although  organ- confined  prostate  cancer  is  highly  curable, 
more  than  32,000  U.S.  men  are  expected  to  die  of  metastatic  prostate 
cancer  in  2010  (1).  Multiple  approved  therapies  (and  newer  agents  in 
late-stage  development)  target  the  androgen  signaling  axis  in  meta¬ 
static  disease;  however,  additional  targeted  therapies  are  lacking. 

We  previously  used  a  bioinformatics  approach,  cancer  outlier 
profile  analysis  (COP A),  to  systematically  prioritize  genes  with 
marked  over  expression  in  a  subset  of  cancers  (outlier  expression).  This 
strategy  identified  outlier  expression  of  the  ETS  (erythroblastosis  virus 
E26  transformation-specific)  family  members  ERG  and  ETV1  in  a 
subset  of  prostate  cancers  across  multiple  gene  expression  profiling 
studies.  It  also  led  to  the  discovery  of  recurrent  gene  fusions  involving 
the  5'  untranslated  region  of  the  androgen-regulated  gene  TMPRSS2 
with  ETS  transcription  factors  (ERG,  ETV1 ,  ETV4 ,  or  ETV5)  (2-5). 
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Subsequent  in  vitro  and  in  vivo  studies  have  demonstrated  a  driving 
role  for  ETS  fusions  in  prostate  oncogenesis  and  cancer  progression 
(6-9). 

Subsequently,  we  used  a  “meta- outlier  approach,”  which  used 
COP  A  to  prioritize  genes  that  consistently  showed  high-ranking 
outlier  expression  across  multiple  profiling  studies.  This  approach 
identified  SPINK1  (serine  peptidase  inhibitor,  Kazal  type  1)  as  a 
high-ranking  meta-outlier  in  a  subset  of  prostate  cancer  with  mutually 
exclusive  outlier  expression  of  ERG  and  ETV1  across  multiple  prostate 
cancer  profiling  studies  (10).  SPINK1,  also  known  as  pancreatic  secretory 
trypsin  inhibitor  (PSTI)  or  tumor-associated  trypsin  inhibitor  (TATI), 
encodes  a  56-amino  acid  peptide  thought  to  protect  the  pancreas 
from  autodigestion  by  preventing  premature  activation  of  pancreatic 
proteases  (11).  Apart  from  its  normal  expression  in  pancreatic  acinar 
cells,  SPINK1  mRNA  has  been  reported  to  be  expressed  in  various 
human  cancers  (12-18),  and  increased  serum  SPINK  1  concentration 
has  been  correlated  with  poor  prognosis  in  some  studies  (12,  13,  17). 
The  prostate  gland  also  secretes  a  variety  of  serine  proteases,  most 
notably  the  kallikrein  enzyme  PSA  (prostate-specific  antigen),  but  also 
trypsin  (19).  Thus,  SPINK1  may  have  a  role  in  modulating  the  activity 
of  cancer-related  proteases  in  other  tissues  besides  the  pancreas. 

We  confirmed  the  mutually  exclusive  overexpression  of  SPINK  1 
and  ETS  gene  fusions  using  a  combined  immunohistochemistry 
(for  SPINK1)  and  fluorescence  in  situ  hybridization  (FISH)  (for  ETS 
fusions)  approach  across  multiple  independent  cohorts,  and  demon¬ 
strated  that  SPINK1  outlier  expression  is  associated  with  an  aggressive 
subset  of  prostate  cancers  (10).  We  also  demonstrated  that  SPINK1 
outlier  expression  can  be  detected  noninvasively  in  urine  and  con- 
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tributes  to  a  multiplexed  panel  of  biomarkers,  which  outperforms 
serum  PSA  for  prostate  cancer  diagnosis  in  patients  presenting  for 
needle  biopsy  ( 10 ,  20).  Our  combined  analyses  of  more  than  1500 
prostate  cancer  cases  demonstrated  SPINK1 
outlier  expression  in  -10%  of  all  PSA- 
screened  prostate  cancers,  which  were  in¬ 
variably  negative  for  ETS  gene  fusions 
(SPINK1  +/ETS~ )  (10).  Furthermore,  SPINK1+ 
tumors  show  shorter  PSA  recurrence-free 
survival  in  prostatectomy-treated  patients 
(10)  and  shorter  progression-free  survival 
in  endocrine-treated  patients  (21). 

Unlike  ETS  gene  fusions  that  lead  to 
the  overexpression  of  a  transcription  fac¬ 
tor  (which  are  difficult  to  target  therapeu¬ 
tically),  SPINK  1  encodes  an  extracellular 
secreted  protein  and  thus  is  potentially  more 
amenable  to  therapeutic  targeting.  Here, 
we  qualify  SPINK  1  as  a  therapeutic  target 
in  SPINK1+/ETS~  prostate  cancer  and  dem¬ 
onstrate  the  therapeutic  potential  of  a  mAb 
to  SPINK1  in  preclinical  models.  Addition- 


Fig.  1.  SPINK!  has  oncogenic  effects  in  pros¬ 
tate  cells  in  vitro.  (A)  SPINK1  stimulated  cell 
proliferation  in  SPINK!~/ETS~  cell  lines.  Benign 
immortalized  prostate  cell  line  RWPE  and 
prostate  cancer  cell  lines  DU  145  and  PC3 
(all  SPINK!~/ETS~)  were  untreated  or  treated 
with  rSPINKI  (10  ng/ml).  Cell  proliferation 
was  measured  by  a  WST-1  colorimetric  assay 
at  the  indicated  time  points.  (B)  SPINK1  me¬ 
diates  invasion  of  RWPE  cells  as  measured 
by  Boyden  chamber  Matrigel  invasion  assay. 
RWPE  cells  were  treated  with  rSPINKI  (1 0  ng/ml) 
or  conditioned  media  (CM)  from  22RV1  cells 
(SPINK!+/ETS~).  (C)  As  in  (B),  except  using  22RV1 
cells  transfected  with  siRNA  against  SPINK!. 
SPINK! -sWenced  22RV1  cells  were  further  treated 
with  rSPINKI  (1 0  ng/ml)  or  CM  from  22RV1  cells. 
(D)  SPINK!  expression  in  SPINK!  knockdown 
22RV1  cells  (stable  pooled  shSPINK!  or  stable 
shSPINK!  clone  1 1 )  compared  to  nontargeting 
pooled  stable  control  (shA/5  vector)  cells  by 
qPCR  (transcript)  or  immunofluorescence  using 
an  antibody  against  SPINK1  (protein,  upper 
inset;  600x  magnification).  (E)  Invasion  assay 
using  shSPINK!  and  shA/5  cells.  Representative 
photomicrographs  (400x  magnification)  show¬ 
ing  cell  motility  assay  (top  inset)  are  shown. 
shA/5  vector  cells  exhibit  longer  cell  motility 
tracks  compared  to  shSPINK!  knockdown  cells. 
(F)  Cell  proliferation  assay  using  pooled 
shSPINK!,  shSPINK!  clone  1 1,  or  shA/5  cells  at 
the  indicated  time  points.  (G)  Soft  agar  colony 
assay  using  pooled  shSPINK!  and  shA/5  cells. 
All  experiments  were  independently  performed 
in  triplicate.  Data  represent  means  ±  SEM.  P 
values  from  significant  two-sided  Student's 
t  tests  are  given  (*P  <  0.05;  **P  <  0.001). 


ally,  we  demonstrate  that  SPINK1  mediates  its  oncogenic  effects  in  part 
through  epidermal  growth  factor  receptor  (EGFR)  and  that  a  mAb  to 
EGFR  shows  in  vitro  and  in  vivo  activity  in  SPINK1+  prostate  cancer. 
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RESULTS 

SPINK1  as  an  autocrine  factor  in  prostate  cancer 

To  further  investigate  the  role  of  SPINK  1  in  prostate  cancer,  we  deter¬ 
mined  the  effects  of  exogenous  SPINK  1  on  invasion  and  proliferation 
using  recombinant  hexahistidine  (6XHis) -tagged  SPINK1  protein 
(rSPINKl)  (fig.  SI  A)  or  conditioned  media  (CM)  collected  from  22RV1 
prostate  cancer  cells  (SPINK  1  +/ETS~)  (fig.  SIB)  (10).  We  treated  be¬ 
nign  immortalized  RWPE  prostate  epithelial  cells  and  DU145  and  PC3 
prostate  cancer  cells  (both  of  which  are  SPINK1~/ETS~)  with  rSPINKl 
(10  ng/ml),  which  resulted  in  a  significant  increase  in  cell  proliferation 
(Fig.  1  A).  We  next  characterized  the  effect  of  rSPINKl  or  22RV 1  CM  on 
cell  invasion  using  a  Boyden  chamber  Matrigel  invasion  assay.  As  shown 
in  Fig.  IB,  addition  of  rSPINKl  or  22RV 1  CM  to  RWPE  cells  significantly 
increased  invasion  (P  =  0.003  and  0.0009,  respectively).  Similar  effects 
were  observed  when  MCF7  breast  cancer  cells  were  treated  with  rSPINKl 
or  22RV1  CM  (fig.  SIC).  Multiple  recombinant  6XHis-tagged  control 
proteins  or  CM  collected  from  RWPE  or  LNCaP  prostate  cancer  cells 
did  not  induce  invasion  in  RWPE  cells  (figs.  SID  and  S2). 

We  previously  showed  that  transient  small  interfering  RNA 
(siRNA)-mediated  knockdown  of  SPINK1  in  22RV1  cells  decreased 
cell  invasion  (10).  Here,  we  extended  these  results  by  demonstrating 
that  the  addition  of  rSPINKl  or  22RV1  CM  rescued  the  invasive  phe¬ 
notype  of  22RV1  cells  in  which  SPINK1  was  knocked  down  (Fig.  1C; 
P  =  0.001  for  both  rSPINKl  and  22RV1  CM). 
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Fig.  2.  An  antibody  to  SPINK1  attenuates  in  vitro  proliferation  and  invasion  exclusively  in  SPINK1+/ETS 
prostate  cancer  cells.  (A)  Cell  proliferation  of  DU  145,  PC3,  and  22RV1  cells  was  assessed  in  the  presence  of 
SPINK1  mAb  or  IgG  mAb  (1  pg/ml).  (B)  As  in  (A),  except  using  22RV1  cells  and  SPINK1  mAb  or  IgG  mAb 
(0.5  to  1  |ig/ml).  (C)  Effect  of  SPINK1  mAb  or  IgG  mAb  on  invasion  of  SPINK1+/ETS~  cells  (22RV1  and 
CWR22Pc)  and  SPINK1+/ETS~  cells  (DU145,  PC3,  LNCaP,  and  VCaP).  All  experiments  were  independently 
performed  in  triplicates.  Data  represent  means  ±  SEM.  Rvalues  from  significant  two-sided  Student's  t  tests 
are  given  (*P  <  0.05;  **P  <  0.001). 


We  next  investigated  whether  the  exogenous  effect  of  SPINK  1  on 
cell  proliferation  and  invasion  was  dependent  on  protease  inhibitory 
activity  of  trypsin  [which  has  been  shown  to  be  simultaneously  ex¬ 
pressed  with  SPINK  1  in  different  tumor  types  (17,  22)]  or  PSA.  Initial 
experiments  demonstrated  that  PRSS1  (trypsinogen)  mRNA  expres¬ 
sion  in  22RV 1  cells  is  relatively  low  compared  with  the  CAPAN-1  pan¬ 
creatic  cancer  cell  line  (fig.  S3  A),  although  a  significant  increase  in 
PRSS1  transcript  was  observed  in  siRNA-mediated  SPINK1  knock¬ 
down  22RV1  cells  (fig.  S3B).  However,  as  shown  in  fig.  S3C,  stimu¬ 
lation  of  22RV1  cells  with  rSPINKl  or  EGF  did  not  affect  trypsin 
expression.  siRNA-mediated  knockdown  of  PRSS1  in  22RV1  cells 
also  had  no  effect  on  invasion  (fig.  S3,  D  and  E).  Similarly,  stimula¬ 
tion  of  22RV 1  cells  with  rSPINKl  or  EGF  did  not  significantly  affect 
PSA  expression  (fig.  S4A).  Finally,  blocking  PSA  with  a  mAb  did  not 
significantly  inhibit  22RV1  cell  invasion  (fig.  S4B).  Together,  these 
findings  demonstrate  that  extracellular  SPINK  1  induces  prostate 
cancer  cell  proliferation  and  invasion  independent  of  protease  inhib¬ 
itory  activity  of  trypsin  or  PSA.  Although  effects  on  other  proteases 
cannot  be  excluded,  our  results  suggest  that  SPINK  1  is  an  autocrine 
pro-proliferative  and  proinvasive  factor  with  effects  independent  of 
trypsin  and  PSA  activity. 


The  role  of  SPINK1  in  cell  proliferation  and  invasion 

To  further  investigate  the  role  of  SPINK1  in  cell  proliferation  and 
invasion,  we  generated  short  hairpin  RNA  (shRNA)  against  SPINK1 
and  established  stable  22RV 1  cells  where 
SPINK  1  was  silenced  (sh SPINK1).  Knock¬ 
down  of  SPINK1  in  both  pooled  and 
clonal  sh SPINK1  cells  compared  to  non¬ 
targeting  control  cells  (shNS  cells)  was 
confirmed  at  the  RNA  level  by  quantita¬ 
tive  polymerase  chain  reaction  (qPCR) 
(more  than  80%  in  both),  as  well  as  at 
the  protein  level  by  immunofluorescence 
staining  with  an  antibody  against  SPINK1 
(Fig.  ID).  Next,  we  investigated  the  role 
of  SPINK1  in  cell  invasion  and  motility 
using  sh SPINK1  cells.  As  anticipated, 
sh SPINK1  cells  showed  decreased  cell 
invasion  by  more  than  75%  in  a  Boyden 
chamber  Matrigel  assay  compared  to  non¬ 
specific  vector  control  (shNS)  cells  (Fig. 
IE;  P  =  0.002).  Reduction  of  cell  motility 
in  a  bead  motility  assay  was  also  ob¬ 
served  in  shSPINKl  cells  compared  to 
shNS  cells  (Fig.  IE,  top  panel). 

To  investigate  the  role  of  SPINK1  in 
cell  proliferation,  we  carried  out  assays 
using  pooled  shSPINKl ,  the  clone  with  the 
greatest  SPINK1  knockdown  (shSPINKl 
clone  11),  and  sh NS  cells.  Both  pooled 
(55%  reduction)  and  clonal  shSPINKl 
cells  (66%  reduction)  showed  significant¬ 
ly  decreased  proliferation  compared  to 
sh NS  cells  (Fig.  IF;  P  =  0.00002  in  both 
cases).  Further,  shSPINKl  cells  showed 
decreased  soft  agar  colony  formation  when 
compared  to  sh NS  cells  (Fig.  1G). 
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Fig.  3.  SPINK1  mediates  its  oncogenic  effects  in  part  through  EGFR.  (A)  Immunoprecipitation  using 
antibodies  to  IgG,  SPINK1,  or  GST  of  exogenous  SPINK1-GST,  GST,  or  GST-VEGFR  added  to  HEK  293 
cells  transfected  with  EGFR  and  immunoblotted  with  an  antibody  to  EGFR  (top  panel),  and  immuno¬ 
precipitation  using  antibodies  to  IgG  or  SPINK1  of  exogenous  SPINK1-GST  added  to  22RV1  cells  and 
immunoblotted  with  an  antibody  to  EGFR  (bottom  panel).  (B)  Western  blot  showing  EGFR  phospho¬ 
rylation  in  response  to  rSPINKl  (100  ng/ml)  or  EGF  (10  ng/ml)  stimulation.  (C)  Invasion  assay  showing 
siRNA-mediated  EGFR  knockdown  22RV1  cells  treated  with  rSPINKl  (10  ng/ml).  (D)  Same  as  in  (C),  ex¬ 
cept  with  RWPE  cells.  (E)  Invasion  assay  showing  rSPINKl  (10  ng/ml)-stimulated  RWPE  cells  in  the 
presence  or  absence  of  C225  [cetuximab  (50  pg/ml)]  or  IgG  mAb  (50  pg/ml).  (F)  Invasion  assay  show¬ 
ing  the  effect  of  IgG  or  C225  antibody  on  SPINK1+  and  SPINKr  cancer  cells.  (G)  As  in  (F),  except  22RV1 
cells  were  treated  with  a  combination  of  antibodies  to  SPINK1  (1  pg/ml)  and/or  C225  (50  pg/ml).  (H)  Cell 
proliferation  assay  using  the  indicated  cells  in  the  presence  of  IgG  mAb  or  C225.  All  experiments  were 
independently  performed  in  triplicates.  Data  represent  means  ±  SEM.  P  values  from  significant  two-sided 
Student's  t  tests  are  given  (*P  <  0.05;  **P  <  0.001). 


In  vitro  targeting  of  SPINK1 
using  a  mAb 

Because  our  results  above  demonstrate  a 
role  for  SPINK1  in  invasion  and  prolifer¬ 
ation,  and  SPINK  1  is  an  extracellular  se¬ 
creted  protein,  we  hypothesized  that  a 
mAb  against  SPINK  1  may  be  able  to  di¬ 
rectly  target  SPINK1+/ETS~  prostate  can¬ 
cer  cells.  Thus,  we  tested  the  effects  of 
an  antibody  to  SPINK1  on  22RV1  cell 
proliferation  and  invasion.  The  SPINK1 
mAb  (0.5  and  1  pg/ml)  significantly  in¬ 
hibited  22RV1  cell  proliferation  by  40 
and  50%,  respectively,  compared  to  a  con¬ 
trol  monoclonal  immunoglobulin  G  (IgG) 
antibody  (Fig.  2,  A  and  B;  P  =  0.0001  and 
P  =  0.0007,  respectively).  However,  the 
antibody  to  SPINK  1  had  no  effect  on 
DU145  and  PC3  cell  proliferation. 

In  addition  to  inhibiting  proliferation, 
the  mAb  to  SPINK1  (0.5  and  1  pg/ml)  sig¬ 
nificantly  attenuated  cell  invasion  by  69 
and  81%,  respectively,  compared  to  a  con¬ 
trol  IgG  mAb  in  22RV 1  cells  (Fig.  2C;  P  = 
0.002  and  P  =  0.007,  respectively).  Similar 
to  22RV1,  which  is  an  androgen  signaling- 
independent  derivative  of  primary  CWR22 
human  prostate  xenograft  tumors,  we  also 
investigated  CWR22Pc  cells,  an  androgen 
signaling-dependent  derivative  of  CWR22 
(23),  which  also  express  high  amounts  of 
SPINKl.  As  expected,  CWR22Pc  cell  in¬ 
vasion  was  blocked  by  47  and  54%  by  the 
mAb  to  SPINKl  at  0.5  and  1  pg/ml  of 
SPINKl  mAb  concentration  (Fig.  2C;  P  = 
0.003  and  P  =  0.002,  respectively).  The 
mAb  to  SPINKl  had  no  significant  effect 
on  invasion  of  SPINKF  prostate  cancer  cell 
lines  including  PC3,  DU145,  LNCaP,  or 
VCaP  (Fig.  2C).  Finally,  the  mAb  to 
SPINKl  attenuated  22RV1  cell  motility 
compared  to  IgG  control,  but  had  no  ef¬ 
fect  on  PC3  (SPINKF ZETST)  cell  motility 
(fig.  S5A). 

Oncogenic  effects  of  SPINKl 
in  part  through  interaction 
with  EGFR 

SPINKl  has  a  similar  structure  as  EGF, 
with  -50%  sequence  homology  and  three 
intrachain  disulfide  bridges  (24,  25).  To 
characterize  potential  SPINKl  and  EGFR 
interaction,  we  overexpressed  EGFR  in 
human  embryonic  kidney  (HEK)  293  cells 
and  incubated  the  lysates  with  SPINK  1- 
GST  (glutathione  S-transferase),  GST,  or 
GST-VEGF  (vascular  endothelial  growth 
factor)  receptor  2  (GST-VEGFR)  recom- 
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binant  proteins.  We  observed  a  strong  interaction  between  SPINKl  - 
GST  and  EGFR  but  not  with  GST  alone  or  GST-VEGFR  recombi¬ 
nant  protein  (Fig.  3 A,  top  panel).  Endogenous  SPINK  1  and  EGFR 
interaction  was  not  detected  by  immunoprecipitation  and  immuno- 
blotting  in  22RV 1  cells,  because  of  the  secretory  nature  of  the  SPINK1 
protein.  However,  addition  of  GST-SPINK1  to  22RV1  cells  followed 
by  immunoprecipitation  and  immunoblotting  confirmed  the  inter¬ 
action  of  SPINK1  and  endogenous  EGFR  in  22RV1  cells  (Fig.  3A, 
bottom  panel). 

To  further  delineate  the  role  of  EGFR  mediation  of  SPINK  1  in 
prostate  cancer,  we  next  assessed  whether  exogenous  SPINK1  was 
capable  of  inducing  EGFR  phosphorylation  (similar  to  the  cognate 
ligand  EGF).  Stimulating  22RV1  cells  with  rSPINKl  resulted  in 
EGFR  phosphorylation,  although  weaker  than  that  observed  with 
EGF  (Fig.  3B).  rSPINKl  stimulation  resulted  in  sustained  EGFR  phos¬ 
phorylation  over  a  90-min  time  course,  whereas  EGF  resulted  in  strong 
EGFR  phosphorylation,  which  diminished  after  only  10  min.  Simi¬ 
larly,  stable  sh SPINK1  knockdown  22RV 1  cells  (pooled  and  clonal) 
showed  decreased  phosphorylated  EGFR  (pEGFR),  with  slightly  de¬ 
creased  total  EGFR  (possibly  because  of  EGFR  degradation)  (fig. 
S6A).  Finally,  we  demonstrate  that  rSPINKl  is  able  to  induce  dimeriza¬ 
tion  of  EGFR,  although  more  weakly  than  EGF  (fig.  S6B). 

We  next  examined  the  functional  consequences  of  SPINK1-EGFR 
interaction  in  the  context  of  SPINK1+  prostate  cancer  using  22RV 1  cells. 
Transient  knockdown  of  EGFR  (fig.  S5B)  blocked  22RV 1  cell  invasion 
by  75%  (Fig.  3C;  P  =  0.004),  which  was  partially  rescued  by  addition 
of  exogenous  SPINK1.  A  similar  effect  of  EGFR  knockdown  was 
observed  in  RWPE  cells  treated  with  rSPINKl  (Fig.  3D;  P  =  0.014 
and  P  =  0.021,  respectively).  These  results  suggest  that  some  but 
not  all  of  SPINK  l’s  effects  are  mediated  by  EGFR. 

Because  mAbs  to  EGFR  are  Food  and  Drug  Administration  (FDA)- 
approved  for  certain  cancers,  we  sought  to  determine  whether  EGFR 
blockade  could  inhibit  the  oncogenic  effects  of  SPINK1.  We  first 
demonstrated  that  mAb  to  EGFR  (cetuximab,  C225)  blocked  the  cell- 
invasive  effects  of  rSPINKl  in  RWPE  cells  (Fig.  3E).  C225  also 
blocked  cell  invasion  of  SPINK1+  22RV1  cells  but  not  in  SPINKF  cell 
lines  DU145,  PC3,  LNCaP,  or  VCaP  (Fig.  3F).  Combining  mAbs  to 
SPINK  1  and  EGFR  had  an  additive  effect  in  the  inhibition  of  22RV1 
cell  invasion  (Fig.  3G;  P  =  0.001).  In  contrast  to  mAb  to  SPINK1  (Fig. 
2A),  C225  had  no  effect  on  22RV1  cell  proliferation  or  PC3  and 
DU  145  cell  proliferation  (Fig.  3H).  Together,  these  experiments  sug¬ 
gest  that  SPINK1  has  both  EGFR-dependent  and  EGFR-independent 
functions  in  prostate  cancer. 

As  a  preliminary  exploration  of  the  downstream  signaling  path¬ 
ways  involved  in  the  SPINK  1 -EGFR  axis,  we  studied  the  mitogen- 
activated  protein  kinase  (MAPK)  and  protein  kinase  B/AKT  pathways 
in  stable  SPINK1  knockdown  22RV1  cells  (sh SPINK1  clone  11).  We 
observed  decreased  pMEK  (phosphorylated  mitogen-activated  or  ex¬ 
tracellular  signal-regulated  protein  kinase  kinase),  pERK  (phosphoryl¬ 
ated  extracellular  signal-regulated  kinase),  and  pAKT  (phosphorylated 
AKT)  in  stable  sh SPINK1  cells  compared  to  control  shZVS  cells  (fig.  S5C). 
Likewise,  22RV 1  cells  treated  with  SPINK  1  mAb  antibody  showed  de¬ 
creased  pERK  (fig.  S5D).  These  observations  provide  the  foundation 
for  further  studies  of  the  SPINK  1 -EGFR  axis. 

The  role  of  SPINK1  in  vivo  and  as  a  therapeutic  target 

Our  in  vitro  studies  demonstrated  that  SPINK1  mediates  cell  prolifer¬ 
ation  and  invasion  in  SPINK1+  prostate  cancer  cells,  and  suggested  that 


a  mAb  can  target  extracellular  SPINK  1.  To  investigate  the  role  of 
SPINK1  in  intravasation,  a  key  step  involved  in  the  process  of  metas¬ 
tasis,  we  used  a  chick  chorioallantoic  membrane  (CAM)  model  sys¬ 
tem  (26)  and  demonstrate  that  rSPINKl  induced  intravasation  of 
benign  RWPE  cells  (Fig.  4A).  Similarly,  SPINK1  mAb  and  C225  sig¬ 
nificantly  inhibited  22RV1  cell  intravasation  (P  =  0.01  and  P  =  0.03, 
respectively),  but  did  not  significantly  inhibit  PC3  cell  intravasation 
(Fig.  4,  B  and  C). 

To  qualify  SPINK  1  as  a  potential  therapeutic  target  in  vivo,  we 
implanted  pooled  shSPINKl  -luciferase  (luc)  and  shNS-luc  22RV1 
cells  in  nude  male  mice.  At  both  4  and  5  weeks  after  implantation, 
22RV 1  -shSPINKl -luc  cells  formed  significantly  smaller  tumors  (55% 
reduction  at  week  4,  P  =  0.008,  and  63%  reduction  at  week  5,  P  = 
0.013)  compared  to  shMS-luc  cells  (Fig.  4,  D  and  H). 

To  demonstrate  preclinical  efficacy  of  the  mAb  to  SPINK1,  we 
treated  nude  mice  implanted  with  22RVl-luc  cells  with  either  the 
mAb  to  SPINK  1  or  an  isotype-matched  monoclonal  IgG  (10  mg/kg) 
twice  a  week.  As  shown  in  Fig.  4,  E  and  I,  administration  of  SPINK  1 
mAb  monotherapy  resulted  in  a  61%  reduction  of  tumor  burden  at 
week  4  (P  =  0.015)  and  58%  reduction  at  week  5  (P  =  0.015).  A  signif¬ 
icant  decrease  in  Ki-67-positive  immunostained  nuclei  was  observed 
in  the  SPINK  1  mAb-treated  group  compared  to  the  control  group 
(fig.  S7). 

Because  SPINK  1  mediates  its  oncogenic  effects  in  part  through 
EGFR,  we  similarly  assessed  the  mAb  to  EGFR  (C225)  using  the 
same  dosage  schedule.  C225  treatment  resulted  in  a  41%  reduction 
at  week  4  (P  =  0.04)  and  37%  reduction  at  week  5  (P  =  0.02)  (Fig.  4, 
E  and  I).  By  combining  mAbs  to  SPINK  1  and  EGFR,  we  observed 
an  additive  effect  in  vivo  showing  a  74  and  73%  reduction  in  the  growth 
of  22RV1  xenografts  at  weeks  4  (P  =  0.01)  and  5  (P  =  0.003),  respec¬ 
tively  (Fig.  4,  F  and  I). 

To  confirm  our  in  vitro  results,  which  suggested  no  effect  of 
SPINK  1  or  EGFR  inhibition  on  SPINKF  prostate  cancer,  we  per¬ 
formed  a  similar  xenograft  study  using  PC3  cells.  As  expected,  neither 
SPINK  1  mAb  nor  C225  significantly  inhibited  tumor  growth  in  PC3 
xenografted  mice  (Fig.  4,  G  and  I).  Finally,  to  investigate  the  potential 
toxicity  of  SPINK  1  mAb  therapy,  we  investigated  whether  the  mAb  to 
SPINK1  interacts  with  SPINK3,  the  murine  homolog  of  SPINKl.  The 
mAb  to  SPINKl  used  in  our  studies  does  not  recognize  murine 
SPINK3,  thus  explaining  the  lack  of  observed  toxicity  in  SPINKl 
mAb-treated  mice  (fig.  S8,  A  to  C). 


DISCUSSION 

Previous  studies  demonstrated  that  SPINKl  outlier  expression  identified 
a  subset  of  PTS-negative  prostate  cancers  (-10%  of  all  PSA-screened 
prostate  cancers),  although  the  mechanism  for  SPINKl  outlier  expres¬ 
sion  remains  unknown  (10).  SPINKl  defines  a  distinct  molecular  sub- 
type  of  prostate  cancer  characterized  by  lack  of  ETS  gene  fusions  as 
well  as  a  more  aggressive  phenotype  as  corroborated  by  independent 
groups  across  distinct  cohorts  of  prostate  cancer  patients  (16, 21).  Thus, 
our  working  hypothesis  is  that  SPINKF  prostate  cancer  represents  an 
aggressive  form  of  prostate  cancer  that  may  respond  to  different  ther¬ 
apies  than  ETS  gene  fusion-positive  prostate  cancers. 

Here,  we  show  that  SPINKl  promotes  prostate  cancer  proliferation 
and  invasion  through  autocrine  and  paracrine  signaling.  We  also  dem¬ 
onstrate  an  in  vivo  role  for  SPINKl  in  intravasation  and  tumor  xeno- 
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graft  growth.  At  present,  the  precise 
mechanism  and  signaling  pathways  re¬ 
sponsible  for  these  effects  in  SPINK1+ 
prostate  cancer  are  unclear.  A  recent 
study  showed  that  mutation  of  SPINK1 
at  leucine  18  (L18)  in  the  trypsin  inter¬ 
action  site  reduced  tumor  growth,  angi¬ 
ogenesis,  and  lung  metastases  in  HT-29 
5M21  human  colon  carcinoma  tumor 
xenografts,  suggesting  that  the  cancer- 
related  phenotypes  of  SPINK1  may  be 
related  to  its  anti-proteinase  activity  (27). 
Moreover,  the  invasive  behavior  of  these 
HT-29  5M21  colon  cancer  cells  was  abol¬ 
ished  with  an  antibody  to  SPINK  1  (27). 
However,  in  our  study,  we  did  not  observe 
any  effect  of  SPINK  1  on  trypsin  or  PSA, 
two  candidate  proteases  in  prostate  cancer. 

Recent  studies  also  indicate  that  SPINK1 
may  be  an  apoptosis  inhibitor  prevent  - 


Fig.  4.  SPINK1  is  a  therapeutic  target  in 
SPINK1+  prostate  cancer.  (A)  Chick  chorio¬ 
allantoic  membrane  (CAM)  assay  quantify¬ 
ing  intravasated  RWPE  cells  upon  stimulation 
with  rSPINKI  {n  =  6  in  each  group).  (B)  CAM 
assay  using  22RV1  cells  in  the  presence  of 
IgG  mAb,  SPINK1  mAb,  or  C225  (n  =  5  in 
each  group),  with  fold  change  of  intravasated 
cells  compared  to  IgG  mAb  plotted.  (C)  As  in 
(B),  except  using  PC3  cells.  (D)  Subcutaneous 
xenograft  growth  of  sh A/5-1  uciferase  (luc)  or 
sh5P/A//C7-luc  22RV1  cells  implanted  in  male 
BALB/c  nu/nu  mice  [n  =  10  in  each  group). 
(E)  As  in  (D),  except  using  22RV1-luc  cell  xeno¬ 
grafts  treated  with  control  IgG  mAb  {n  =  8), 
SPINK1  mAb  (n  =  6),  or  C225  (n  =  8)  (1 0  mg/kg) 
twice  a  week.  (F)  Same  as  in  (E),  except  mice 
{n  =  7  per  group)  were  treated  with  a  com¬ 
bination  of  SPINK1  and  C225  mAb  (10  mg/kg 
for  both).  (G)  As  in  (E)  and  (F),  except  using 
PC3-luc  xenografts  treated  with  control  IgG 
mAb,  SPINK1  mAb,  or  C225  (n  =  8  per  group) 
(10  mg/kg)  alone  or  in  combination  twice  a 
week.  (H)  Representative  bioluminescence 
images  from  mice  in  (D)  bearing  pooled 
shA/5-luc  or  shSPINKl-luc  xenografts  and  per¬ 
cent  reduction  in  tumor  volume  at  week  5.  (I) 
Same  as  (H),  except  bioluminescence  images 
from  mice  bearing  22RV1-luc  xenografts 
(red,  top  panel)  or  PC3-luc  (blue,  lower  panel) 
mice  treated  with  IgG  mAb,  SPINK1  mAb,  or 
C225  mAb  alone  or  in  combination,  with 
comparative  percent  reduction  plot  in  tumor 
volume  at  week  5.  Data  represent  means  ± 
SEM.  P  values  from  significant  two-sided 
Student's  t  tests  are  given  (*P  <  0.05;  **P  < 
0.001). 
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ing  serine  protease-dependent  cell  death  (28).  Here,  we  show  that 
SPINK1,  which  has  structural  similarities  with  EGF  (29),  binds  to 
EGFR,  and  inhibiting  SPINK  1  attenuates  key  downstream  mediators 
of  the  EGFR  pathway  including  MEK,  ERK,  and  AKT.  Furthermore, 
we  also  show  that  SPINK1  dimerizes  EGFR  and  induces  sustained 
phosphorylation  of  EGFR,  which  have  been  shown  to  be  critical  for 
downstream  signaling  activation  after  ligand  binding  (30).  However, 
in  contrast  to  SPINK  1  mAh,  EGFR  mAh  only  partially  inhibited  the 
cell-invasive  effects  of  22RV 1  cells  and  had  no  effect  on  cell  proliferation, 
suggesting  that  SPINK  1  engages  both  EGFR-dependent  and  EGFR- 
independent  pathways  to  mediate  its  oncogenic  effects.  SPINK  1  has  also 
been  shown  to  engage  the  EGFR/MAPK  cascade  in  NIH  3T3  fibroblasts 
and  pancreatic  cancer  cells  (31). 

This  study  provides  compelling  evidence  that  SPINK1  overexpres¬ 
sion  is  oncogenic  in  prostate  cancer  and  that  inhibition  of  SPINK1  via 
RNA  interference  or  blocking  antibodies  may  have  therapeutic  po¬ 
tential.  Our  preclinical  models  suggest  that  this  therapeutic  effect 
would  only  be  effective  in  patients  with  SPINK1+  prostate  cancer,  sug¬ 
gesting  that  such  therapies  would  need  to  be  evaluated  in  a  molec- 
ularly  guided  fashion.  Because  the  area  of  antibody-based  therapeutics 
for  extracellular  targets  is  well  developed,  based  on  examples  such  as 
trastuzumab  in  breast  cancers  with  ERBB2  overexpression,  we  postu¬ 
late  that  a  SPINK  1 -blocking  antibody  may  have  similar  efficacy  on  a 
molecularly  defined  subset  of  prostate  cancers.  We  have  previously 
demonstrated  that  patients  with  the  subset  of  SPINK 1+/ETS~  prostate 
cancers  can  be  reliably  identified  by  immunohistochemistry  ( 10 ,  20 ), 
as  would  be  required  for  a  molecularly  defined  clinical  trial.  Although 
humanized  SPINK  1  mAbs  are  not  yet  available  for  clinical  testing,  our 
studies  show  that  SPINK  1  partially  mediates  its  oncogenic  effects 
through  EGFR. 

This  finding  prompted  us  to  evaluate  the  utility  of  the  FDA-approved 
EGFR  mAh  cetuximab,  which  showed  in  vitro  and  in  vivo  activity  only 
against  SPINK1+  prostate  cancer  cells  (although  less  effective  than 
SPINK  1  mAh).  Phase  I/II  clinical  trials  of  cetuximab  (32)  and  EGFR 
small  molecules  have  been  largely  disappointing  in  metastatic  prostate 
cancer  (33,  34);  however,  a  small  subset  of  patients  have  had  responses, 
including  3  of  36  (8%)  patients  who  showed  >50%  PSA  decline  in  a 
Phase  Ib/IIa  clinical  trial  of  cetuximab  in  combination  with  doxorubicin 
in  castrate-resistant  metastatic  prostate  cancer  patients  (32).  Results 
from  our  study  provide  a  plausible  mechanism  for  why  only  the  lim¬ 
ited  subset  of  patients  with  positive  cancers  (-10%  of  all  cases)  may 
benefit  from  EGFR  inhibition.  This  hypothesis  can  be  assessed  ret¬ 
rospectively  and  in  biomarker-informed  clinical  trials  of  patients 
with  SPINK1+  prostate  cancer.  Because  the  mAh  to  SPINK1  used 
in  our  studies  did  not  interact  with  murine  SPINK3  (the  homolog 
of  SPINK1),  our  study  does  not  inform  on  the  potential  toxicity  of 
SPINK  1  mAh  therapy.  However,  an  FDA-approved  mAh  to  EGFR 
has  specific  in  vivo  activity  against  SPINK1+  prostate  cancer,  provid¬ 
ing  an  immediately  translatable  strategy  for  targeting  SPINK1+  cancers 
that  can  be  clinically  investigated  while  toxicity  of  humanized  SPINK  1 
antibody  therapy  is  explored. 

In  summary,  our  results  support  SPINK1  as  an  oncogene  in  a  sub¬ 
set  of  prostate  cancers  that  can  be  molecularly  identified,  and  provide 
the  rationale  to  develop  humanized  SPINK1  antibodies  for  human 
clinical  trials.  Our  work  also  reinforces  the  molecular  subclassification 
of  prostate  cancer  in  clinical  trials  (whether  through  SPINK/ETS  status 
or  other  relevant  biomarkers),  which  has  lagged  behind  other  com¬ 
mon  epithelial  cancers  (that  is,  breast,  lung,  and  colon). 


MATERIALS  AND  METHODS 

Cell  lines  and  SPINK1  knockdown 

The  benign  immortalized  prostate  cell  line  RWPE  as  well  as  pros¬ 
tate  cancer  cell  lines  DU145,  PC3,  and  22RV1  were  obtained  from 
the  American  Type  Culture  Collection  (ATCC)  and  were  grown  ac¬ 
cording  to  ATCC  guidelines.  For  stable  knockdown  of  SPINK1 , 
human  lentiviral  shRNAmir  individual  clone  (ID  V2LHS_153419) 
targeting  against  SPINK  1  or  nonsilencing  lentiviral  shRNAmir  in 
GIPZ  vectors  was  purchased  from  Open  Biosystems  (Thermo  Scientific 
Open  Biosystems).  Details  are  available  in  Supplementary  Materials 
and  Methods. 

Quantitative  PCR 

Total  RNA  was  isolated  with  a  miRNeasy  mini  kit  following  the 
manufacturer’s  instruction  (Qiagen).  Complementary  DNA  was 
synthesized  from  1  jug  of  total  RNA  with  Superscript  III  (Invitrogen) 
in  the  presence  of  random  primers.  qPCR  was  performed  with  the 
StepOne  Real-Time  PCR  system  (Applied  Biosystems).  Details  and 
primer  information  are  available  in  Supplementary  Materials  and 
Methods. 

Cell  proliferation  assay 

Proliferation  for  control  and  experimental  cells  was  measured  by  a 
colorimetric  assay  based  on  the  cleavage  of  the  tetrazolium  salt  WST-1 
by  mitochondrial  dehydrogenases  (cell  proliferation  reagent  WST-1; 
Roche  Diagnostics)  at  the  indicated  time  points  in  triplicate.  Cell  counts 
for  sh NS  vector  and  sh SPINK1  cells  were  estimated  by  trypsinizing  cells 
and  analysis  by  Coulter  counter  (Beckman  Coulter)  at  different  time 
points  in  triplicates. 

Basement  membrane  matrix  invasion  assay 

For  invasion  assays,  sh NS  vector-  or  shSPINKl  -transduced  cells,  as 
well  as  RWPE,  PC3,  and  22RV1  cells  were  used.  Equal  numbers  of 
the  indicated  cells  were  seeded  onto  the  basement  membrane  matrix 
(BD  Biosciences)  present  in  the  insert  of  a  24-well  culture  plate.  RPMI 
media  supplemented  with  10%  fetal  bovine  serum  were  added  to  the 
lower  chamber  as  a  chemoattractant.  After  48  hours,  noninvading  cells 
and  extracellular  matrix  were  removed  with  a  cotton  swab.  Invaded  cells 
were  stained  with  crystal  violet  and  photographed.  The  inserts  were 
treated  with  10%  acetic  acid,  and  absorbance  was  measured  at  560  nm. 

CAM  assay 

The  assay  was  performed  essentially  as  described  (26).  Two  million 
RWPE  cells  were  mixed  with  either  200  ng  of  multiple  tag  control 
protein  or  200  ng  of  rSPINKl  protein  and  applied  to  the  CAM  of 
11 -day-old  chicken  embryo.  Similarly,  2  million  22RV1  or  PC3 
cells  were  mixed  with  either  monoclonal  IgG  or  antibodies  to 
SPINK1  or  C225  (1  pg/ml)  and  applied  onto  the  upper  CAM  of  a 
fertilized  chicken  embryo.  Three  days  after  implantation,  the  relative 
number  of  cells  that  intravasate  into  the  vasculature  of  the  lower  CAM 
was  analyzed  by  extracting  genomic  DNA  with  the  Puregene  DNA 
purification  system.  Quantification  of  the  human  cells  in  the  extracted 
DNA  was  done  as  described  (35). 

22RV1  and  PC3  xenograft  models 

Four-week-old  male  BALB/c  nu/nu  mice  were  purchased  from  Charles 
River  Inc.  (Charles  River  Laboratory).  Stable  22RV1  shNS-luc  and  22RV1 
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shSPINKl-luc  cells  (5  x  105),  or  22RVl-luc  (2  x  105)  or  PC3-luc  (5  x 
105)  cells  were  resuspended  in  100  pi  of  saline  with  20%  Matrigel  (BD 
Biosciences)  and  were  implanted  subcutaneously  into  the  left  flank  re¬ 
gions  of  the  mice.  Details  are  available  in  Supplementary  Materials 
and  Methods. 

Statistical  analysis 

All  values  presented  in  the  study  were  expressed  as  means  ±  SEM.  The 
significant  differences  between  the  groups  were  analyzed  by  a  Student’s 
t  test,  and  a  P  value  of  <0.05  or  <0.001  was  considered  significant. 


SUPPLEMENTARY  MATERIAL 

www.sciencetranslationalmedicine.org/cgi/content/full/3/72/72ral7/DC1 
Materials  and  Methods 

Fig.  SI.  rSPINKI  or  CM  collected  from  22RV1  cells  induces  invasion  in  benign  or  cancer  cells. 
Fig.  S2.  CM  collected  from  22RV1  cells  induces  cell  invasion,  but  not  CM,  from  LNCaP  cells. 
Fig.  S3.  PRSS1  (trypsinl)  knockdown  in  22RV1  cells  has  no  effect  on  SPINK1 -mediated  cell  invasion. 
Fig.  S4.  Exogenous  rSPINKI  has  no  effect  on  PSA  in  22RV1  cells. 

Fig.  S5.  SPINK1  mAb  reduces  SPINK! +  cell  motility  and  SPINK1  knockdown  alters  MAPK  pathway. 
Fig.  S6.  Exogenous  SPINK1  induces  EGFR  dimerization  and  phosphorylation. 

Fig.  S7.  SPINK1  mAb  induces  decrease  in  tumor  proliferation  index. 

Fig.  S8.  Anti-SPINKI  mAb,  which  does  not  recognize  the  murine  homolog  of  SPINK1  (SPINK3), 
has  no  observed  toxic  effect  in  treated  mice. 
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which  miR-101  is  dysregulated  in  prostate  cancer. 

Completed  Grant  Support 

GlaxoSmithKline  (GSK)  (PI:  Varambally)  07/01/08  -  06/30/09 

The  Role  of  AVODART  in  Pre-Clinical  Models  of  Gene  Fusion  Positive  Prostate  Cancer 
The  major  goal  of  this  project  is  to  analyze  of  the  effect  of  AVODART  on  fusion  positive  prostate 
cancer  cell  proliferation,  invasion,  and  gene  expression  and  evaluate  the  effect  of  AVODART  on 
tumor  growth  in  a  TMPRSS2-ERG  xenograft  model  of  prostate  cancer. 

Role:  Co-investigator 


Honors  and  Awards 


1994-1995 

1994- 1996 

1995- 1996 
2010-2012 

2011 


Outstanding  Freshmen  Scholarship,  Peking  University,  Beijing,  China 
Fellowship  of  Honor  Science  Program,  Peking  University,  Beijing,  China 
Young  Eagle  Scholarship,  Peking  University,  Beijing,  China 
Department  of  Defense  Prostate  Cancer  Research  Program  Post-Doctoral 
Traineeship  Award 

American  Association  for  Cancer  Research  102nd  Annual  Meeting  AACR- Aflac 
Scholar-in-Training  Award 


Membership  in  Professional  Societies 

Associate  Member,  American  Association  for  Cancer  Research  (AACR) 


Teaching 

Mentor,  Undergraduate  Students 

Arun  Dahiya,  University  of  Michigan  College  of  LS&A,  2009-present 
Krystina  Auble,  University  of  Michigan  College  of  LS&A,  201 1-2012 
Joy  E.  Tsai,  University  of  Michigan  College  of  LS&A,  2008-2009 
Rupal  Shastri,  University  of  Michigan  College  of  LS&A,  2008-2009 
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Mithil  Pandhi,  University  of  Michigan  College  of  LS&A,  2006-2008 


Seminars  and  Extramural  Invited  Presentations 

1.  Multi-Institutional  SPORE  Prostate  Cancer  Program  Retreat,  “miR-101  loss  and  EZH2 
elevation  in  prostate  cancer”,  Baltimore,  MD,  March  25-27,  2009. 

2.  American  Association  for  Cancer  Research  (AACR)  102nd  Annual  Meeting,  “An  onco¬ 
protein  axis  linking  polycomb  repressive  complex  2  and  polycomb  repressive  complex  1 
through  miRNAs  in  cancer”,  Orlando,  FL,  April  2-6,  2011. 

3.  Department  of  Molecular  and  Cellular  Oncology,  The  University  of  Texas  MD  Anderson 
Cancer  Center,  “Coordinated  Regulation  of  Polycomb  Repressive  Complexes  Through 
miRNAs  in  Cancer”,  Houston,  TX,  September  1,  201 1. 

4.  Division  of  Life  Sciences,  Hong  Kong  University  of  Science  and  Technology, 
“Coordinated  Regulation  of  Polycomb  Repressive  Complexes  Through  miRNAs  in 
Cancer”,  Hong  Kong,  China,  September  27,  2011. 

5.  Center  for  Inflammation  and  Epigenetics,  The  Methodist  Hospital  Research  Institute, 
“Coordinated  Regulation  of  Polycomb  Repressive  Complexes  Through  miRNAs  in 
Cancer”,  Houston,  TX,  December  14,  201 1. 

6.  Institute  for  Integrative  Genome  Biology,  University  of  California  Riverside,  “Coordinated 
Regulation  of  Polycomb  Repressive  Complexes  Through  miRNAs  in  Cancer”,  Riverside, 
CA,  January  30,  2012. 

7.  Department  of  Pathology,  Mount  Sinai  School  of  Medicine,  “Coordinated  Regulation  of 
Polycomb  Repressive  Complexes  Through  miRNAs  in  Cancer”,  New  York,  NY,  February 
29,2012. 

8.  Fels  Institute  for  Cancer  Research  &  Molecular  Biology,  Temple  University,  “Coordinated 
Regulation  of  Polycomb  Repressive  Complexes  Through  miRNAs  in  Cancer”, 
Philadelphia,  PA,  March  21,  2012. 

9.  School  of  Biomedical  Sciences,  Chinese  University  of  Hong  Kong,  “Coordinated 
Regulation  of  Polycomb  Repressive  Complexes  Through  miRNAs  in  Cancer”,  Hong 
Kong,  China,  March  13,  2012. 

Internal  Presentations 


1.  University  of  Michigan  Comprehensive  Cancer  Center  Annual  Fall  Symposium,  “EZH2  is 
a  marker  of  aggressive  breast  cancer  and  promotes  neoplastic  transformation  of  breast 
epithelial  cells”,  Ann  Arbor,  MI,  Pctober  10,  2003. 
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Patents 


1.  UM-14513/US-1/ORD:  EZH2  Cancer  Markers.  Provisional  application  12/167,520.  Filed 
7/3/2008.  USPTO  Patent  Application  #:  2009001203 1 . 

2.  UM-30734/US-1/PRO:  MIR-101  Cancer  Markers.  Provisional  application  61/154,541. 
Filed  2/23/2009.  USPTO  Patent  Application  #:  20100311815. 

3.  UM-31379/US-2/PRO:  Compositions  and  Methods  for  Inhibiting  EZH2.  Provisional 
application  13/024,559.  Filed  2/10/2011.  USPTO  Application  #:  20100222420. 

4.  UM-31380/US-2/PRO:  Compositions  and  Methods  for  Inhibiting  MMSET.  Provisional 
application  13/025,290.  Filed  2/1 1/201 1 . 
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